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29 .  ABSTRACT  (Continued) 

"^are  developed  and  applied  at  representative  flight  conditions. 
The  stability  results  indicate  that  a  prime  determinant  of 
longitudinal  stability  in  the  transition  and  nozzle-borne 
flight  regimes  is  nozzle  angle;  this  result  has  an  effect  on 
the  choice  of  trajectories  to  minimize  aircraft  stability 
problems.  The  lateral-directional  stability  is  primarily  a 
function  of  forward  speed  and  angle  of  attack,  with  the  worst 
case  being  high  angle  of  attack,  low- speed  flight.  To  some 
extent,  flight  conditions  which  minimize  lateral-directional 
stability  problems  accentuate  longitudinal  difficulties.  A 
second  major  section  of  this  report  develops  the  theory  of 
command-generator  tracking  (CGT) ,  a  very  powerful  control 
system  design  method.,^  Since  the  CGT  allows  independent  solu¬ 
tions  to  the  stability  augmentation  and  command  augmentation 
problems,  the  designer \is  free  to  employ  the  most  efficient 
techniques  of  solution  to  these  two  problems.  The  application 
of  CGT  design  techniques  to  various  V/STOL  flight  condition 
regimes  includes  the  investigation  of  different  vehicle  com¬ 
mand  inputs.  For  air  combat  maneuvers,  a  normal  acceleration- 
roll  rate  command  input  is  desirable,  while  a  velocity  command 
vector  is  investigated  in  transition  and  thrust-borne  flight . 
Near  hover,  a  position-command  (4-D)  system  is  appropriate,  and 
an  attitude-command  system  is  compared  to  the  above  systems. 

Not  only  do  these  examples  illustrate  the  flexibility  of  the 
CGT  design  approach,  each  controller  provides  insights  into 
the  stability  and  control  characteristics  of  a  V/STOL  Type  B 
aircraft . 
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INTRODUCTION 


1 . 1  BACKGROUND 

Future  missions  of  the  U.S.  Navy  will  involve  in¬ 
creased  use  of  vertical/short- takeof f-and- landing  (V/STOL) 
aircraft  in  a  variety  of  roles.  It  is  important  that  these 
new  aircraft  be  designed  for  maximum  performance,  reliability, 
and  safety.  High-performance  V/STOL  aircraft,  such  as  the 
V/STOL  B  concept  now  receiving  study,  have  special  problems; 
they  must  not  only  hover  and  transition  to  forward  flight,  but 
fly  in  the  air  combat  environment  as  well.  Thus,  the  rela¬ 
tively  efficient  and  well-understood  mechanisms  currently  used 
to  produce  lift  at  hover  ( low-disc- loading  rotors)  can  not  be 
used  in  V/STOL  B  aircraft,  as  they  are  inappropriate  for  flight 
at  high  Mach  number.  Instead,  the  V/STOL  B  aircraft  is  likely 
to  rely  on  high-momentum  jet  thrust  for  hover.  This  introduces 
the  probability  of  severe  stability  and  control  problems  in 
low-speed  flight,  including  "suck-down"  and  ” fountain- flow" 
phenomenon  in  ground  effect. 

The  need  to  provide  lift  via  jet  thrust  in  hover  in¬ 
troduces  the  possibility  of  a  "fringe  benefit"  in  high-speed 
flight;  vectored  thrust  can  be  used  to  enhance  maneuverabi li ty 
in  air  combat.  Flight  tests  and  simulation  of  the  AV-8A 
Harrier  indicate  that  inflight  thrust  vectoring  can  provide 
substantial  gains  for  air-to-air  combat  and  also  may  be  effec¬ 
tive  for  air-to-ground  weapon  delivery  (Refs.  1  and  2).  The 
advantages  accrue  not  only  because  the  thrust  vector  can  be 
redirected,  but  also  because  the  high-performance  V/STOL  air¬ 
craft  has  reaction  control  effectors  for  hovering  control; 
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therefore,  the  aircraft  retains  angular  controllability  through¬ 
out  its  flight  envelope. 

Command  and  stability  augmentation  take  on  added 
significance  for  the  V/STOL  B  aircraft,  for  it  is  truly  a 
control-configured  vehicle  (CCV).  Without  command/stability 
augmentation,  the  V/STOL  B  aircraft  will  almost  certainly 


have  unacceptable  handling  qualities  at  one  or  more  flight 
conditions,  and  it  is  likely  that  maneuvering  will  compound 


these  difficulties  (Ref.  3). 


analysis,  both  of  which  are  well-suited  to  the  complex  nature 
of  V/STOL  flight,  and  they  have  been  carried  to  a  high  level 
of  refinement  in  our  earlier  work.  The  AV-8A  is  chosen  as  a 
baseline  aircraft  because  representative  data  for  that  air¬ 
craft  is  readily  available  (Ref.  6).  The  word  "representative 
should  be  emphasized;  our  proposed  goal  is  a  better  under¬ 
standing  of  generic  V/STOL  B  stability-and-control  character¬ 
istics  rather  than  a  detailed  study  of  a  particular  aircraft. 


1.2  REPORT  OUTLINE 

This  report  is  organized  according  to  the  following 
outline:  Chapter  2  addresses  the  open-loop  stability  of  the 
AV-8A’s  linear,  time- invariant  model  representation  and  the 
trim  conditions  that  occur  in  flight.  Chapter  3  presents 
the  design  and  simulation  results  for  the  command  generator 
tracker-proportional  integral  (CGT-PI)  controller.  Conclu¬ 
sions  and  recommendations  are  presented  in  Chapter  A.  Feed¬ 
forward  CGT  theory  is  developed  in  Appendix  A.  Pilot  command 
models  and  the  eigenvalue/eigenvector  based  "ideal"  aircraft 
model  are  shown  in  Appendix  B.  The  form  of  the  reference 
aircraft  mathematical  model  is  listed  in  Appendix  C.  A  list 
of  symbols  and  abbreviations  is  given  in  Appendix  D. 
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STABILITY  AND  CONTROL  ANALYSIS  IN  LOW-SPEED, 
TRANSITION,  AND  CONVENTIONAL  FLIGHT 


2.1  OVERVIEW 

The  extensive  range  of  flight  conditions  of  the  AV-8A 
Harrier  presents  a  challenge  to  the  aircraft  control  system 
designer  and  to  the  aircraft  pilot.  At  low  speed,  thrust  and 
gravity  forces  are  paramount,  and  the  vehicle  dynamically 
resembles  an  inverted  pendulum.  Aerodynamic  forces  become  im¬ 
portant  through  the  transition,  accompanied  by  drastic  changes 
in  the  flow  field  around  the  aircraft  as  the  nozzles  are 
swiveled  backwards.  At  the  same  time,  the  control  system 
transitions  from  reaction  control  jets  to  conventional  aero¬ 
dynamic  control  surfaces.  Appendix  C  contains  a  complete 
description  of  the  data  used  to  model  the  AV-8A,  and  an  ana¬ 
lysis  of  its  control  system.  Note  that  all  results  presented 
in  this  chapter  are  for  the  "bare  airframe"  (i.e.,  with  the 
stability  augmentation  sj^stem  switched  off). 

Section  2  of  this  chapter  examines  the  complete  range 
of  flight  conditions  from  low-speed  through  transition  to 
conventional  flight,  and  details  the  effect  of  these  velocity 
changes  on  vehicle  stability.  The  implications  of  these 
changes  on  vehicle  lateral  control  are  also  discussed,  as  are 
the  changes  in  vehicle  stability  due  to  ground  effect.  Sec¬ 
tion  3  details  the  manner  in  which  variations  in  angle  of 
attack  and  longitudinal  acceleration  affect  vehicle  stability 
at  an  intermediate  velocity,  while  Section  4  discusses  the 
vehicle  stability  characteristics  to  be  expected  during  a 
ski-jump  launch.  Section  5  summarizes  the  results  of  this 
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chapter,  and  presents  figures  which  give  the  relationship 
between  the  most  important  variables  and  their  effect  on  AV-8A 
stability . 


2.2  FLIGHT  CONDITIONS  AND  VELOCITY-DEPENDENT  STABILITY 
EFFECTS  IN  THE  TERMINAL  AREA 

2.2.1  Trim  Flight  Conditions 

The  flight  conditions  examined  in  this  section  span 
the  range  of  velocities  from  190  kt  (98  m/s)  to  10  kt  (5  m/s); 
included  are  conventional  wing-borne  flight  conditions,  transi 
tion  flight  conditions  and  low-speed  thrust-borne  flight 
conditions.  Figure  2  illustrates  the  angles  of  attack  (a) 
involved  in  these  flight  conditions.  For  conventional  flight 
(0j=O,  i.e.,  nozzles  fully  aft),  a  increases  with  decreasing 
speed . 
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Figure  2 


Angle  of  Attack 
5 


If  15  deg  a  is  the  upper  limit,  then  128  kt  (66  m/s)  is  the 
lower  speed  limit  in  conventional  flight  for  this  vehicle 
weight.  (High-lift  devices  are  not  modeled.)  The  Harrier  can 
fly  at  angles  of  attack  below  the  conventional  curve  by  lower¬ 
ing  the  nozzles  to  make  up  for  lost  wing  lift.  If  a  constant 
0  of  6  deg  is  assumed,  this  trade-off  can  be  effected  at  any 
speed  below  190  kt  (98  m/s).  Figure  3  presents  this  in  terms 
of  the  portions  of  vehicle  weight  supported  by  aerodynamic 
lift  and  propulsive  lift.  The  thrust-produced  vertical  com¬ 
ponent  of  the  supporting  force  is  directly  proportional  to 
throttle  and  nozzle  deflection  angle.  In  general,  the  thrust 
force  also  includes  ram  drag  forces,  although  these  do  not 
contribute  to  the  vertical  component  of  thrust  force  in  hori¬ 
zontal  flight.  The  wing-borne  lift  is  proportional  to  dynamic 
pressure,  and  takes  into  account  the  lift  decrement  caused  by 
the  thrust-induced  flow  field  changes.  Figure  3  also  shows 
zero  nozzle-angle  curves;  conceptually  there  are  a  series  of 
curves  with  varying  wing-borne/thrust-borne  ratios  in  the 
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Figure  3  Wing  Lift/Thrust  Lift  Trade-Off 


transition  speed  regime.  These  curves  represent  different 

combinations  of  a  and  0.. 

J 

The  variation  in  nozzle  angle  with  velocity  at  the 
nominal  angle  of  attack  is  shown  in  Fig.  4.  This  follows  the 
general  shape  of  the  thrust-borne-weight  curve  in  the  previous 
figure.  At  low  speeds,  the  thrust  vector  angle  with  respect 
to  vertical,  which  is  the  sum  of  the  nozzle  angle  relative  to 
its  datum  (0 j  )  ,  the  nozzle  datum  relative  to  the  body  (1.5 
deg)  and  the  body  pitch  angle  (0),  approaches  90  degrees. 

The  engine  speed  (and  hence  thrust)  necessary  to  trim 
at  these  flight  conditions  is  given  in  Fig.  5.  Th •.*  steady- 
state  relationship  typical  of  thrust-borne  flight  (increased 
engine  speed  is  necessary  at  lower  flight  speeds)  is  apparent 
in  this  figure.  The  engine  speeds  needed  for  fully  wing-borne 
flight  are  also  given  in  the  higher  flight  speed  range,  and 
the  difference  is  indicative  of  the  relative  efficiency  of 
wing-borne  flight. 
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Figure  4  Nozzle  Angle  Versus  Velocity 
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Figure  5  Engine  RPM  at  Various  Speeds 

The  stabilator  deflections  in  the  velocity  range  from 
10  kt  (5  m/s)  to  190  kt  (98  m/s)  are  given  in  Fig.  6,  illus¬ 
trating  the  mild  changes  in  pitch  trim  that  occur  between  the 
various  velocities  with  nozzles  down.  Also  of  note  is  the 
large  change  in  stabilator  position  between  the  "nozzles  up" 
position  and  the  "thrust  borne"  flight  condition.  The  aero 
data  tables  show  that  this  shift  is  due  to  the  large  change  in 
stabilator  deflection  for  zero  pitch  moment  caused  by  the  flow 
field  changes  as  the  nozzles  are  lowered.  For  example,  at  5 
deg  angle  of  attack,  zero  aerodynamic  pitch  moment  coefficient 
occurs  at  -6.95  deg  stabilator  and  zero  nozzle  angle,  but  at 
45  deg  nozzle  angle,  zero  aerodynamic  pitch  moment  coefficient 
occurs  at  6.05  deg  stabilator.  This  indicates  that  as  the 
nozzles  are  dropped  to  45  deg  a  significant  nose  up  pitch 
moment  will  arise.  This  requires  the  pilot's  attention  in 
order  to  transition  smoothly  from  wing-borne  to  nozzle-borne 
flight . 
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Figure  6  Stabilator  Deflection  Versus  Velocity 

The  pitch  jets  (which  thrust  downwards  from  either 
the  nose  or  tail  of  the  aircraft)  also  are  important,  and  Fig. 
7  details  the  operation  of  these  jets.  For  purposes  of  trim 
calculation,  the  jets  are  connected  to  the  stabilator  in  the 
same  manner  as  in  the  actual  aircraft  for  the  sas-off  case. 

The  only  exception  is  for  nozzle  angle  less  than  20  deg,  where 
the  trim  results  presented  here  assume  that  the  pitch  jets  are 
completely  off.  In  actuality,  they  are  phased  in  as  the  noz¬ 
zles  transition  from  zero  to  twenty  degrees.  For  this  reason, 
the  trim  control  results  in  this  region  are  represented  by 
dashed  lines. 

The  range  of  flight  conditions  described  here  is 
examined  in  the  following  section  as  to  the  vehicle  stability 
and  control  characteristics. 
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Figure  7  Pitch  Jet  Opening  Versus  Velocity 
2.2.2  Velocity-Dependent  Stability  Variations 

The  variations  in  longitudinal  and  lateral-directional 
stability  over  the  range  of  flight  conditions  from  10  to  190 
kt  (5  to  98  m/s)  are  discussed  in  this  section.  All  trim 
conditions  represent  symmetric  flight;  hence,  the  longitudinal 
linear  dynamics  are  uncoupled  from  the  lateral-directional 
linear  dynamics.  For  nozzles  completely  rearward  (0^  =  0  deg), 
Fig.  8  shows  that  the  longitudinal  modes  consist  of  a  fast, 
well-dampled  short  period  mode  and  a  slow,  poorly-damped 
phugoid  mode.  Both  of  these  are  typical  of  conventional 
aircraft . 


For  nozzles  downwards,  the  longitudinal  modes  are 
completely  different.  At  very  low  speed  (10  kt,  5  m/s), 
the  four  longitudinal  eigenvalues  consist  of  an  unstable 
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Figure  8  Longitudinal  Modes  from  Hover  to  190  kt 

oscillatory  pair,  the  third  mode  (so-called  because  it  is 
inherently  different  from  the  short-period  and  phugoid  os¬ 
cillations),  and  two  stable  real  roots.  The  slow  real  root 
exhibits  primarily  pitch-speed  characteristics,  while  the 
faster  one  involves  angle  of  attack  variation.  This  mode 
becomes  much  faster  as  the  aerodynamic  forces  grow  with  in¬ 
creasing  speed. 

The  unstable  third  mode  is  caused  by  a  static  pitch 
instability  in  much  the  same  way  as  a  eg  too  far  aft  can 
cause  difficulty  for  a  conventional  aircraft.  The  problem 
with  this  aircraft  is  caused  by  both  aerodynamic  and  engine 
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effects.  The  aerodynamic  pitch  moment  due  to  angle  of  attack 
(C  )  is  small  at  low  speeds  and  depends  on  nozzle  angle  so 
that  stabilizing  aerodynamic  moments  are  not  really  provided 
until  the  nozzle  angle  is  reduced  below  30  deg.  The  engine 
further  destabilizes  the  longitudinal  dynamics  because  of  the 
ram  drag  effect.  This  is  a  drag  force  caused  by  the  engine 
inlet  air  flow,  and  it  depends  on  airspeed  and  engine  rpm. 
Because  this  drag  force  is  applied  ahead  of  the  eg,  it  pro¬ 
duces  a  nose-up  moment  which  grows  with  increasing  angle  of 
attack.  The  net  effect  of  the  aerodynamic  and  engine- induced 
effects  is  that  the  vehicle  is  statically  unstable  below  about 
170  kt  (87  m/s). 

The  shapes  of  these  modes  are  shown  in  Fig.  9  for 
both  135  kt  (69  m/s)  and  180  kt  (93  m/s).  At  135  kt  (69  m/s), 
the  unstable  third  mode  will  cause  difficulty.  As  can  be 
seen,  this  mode  is  much  different  from  a  conventional  short 
period  mode.  The  amount  of  velocity  and  pitch  angle  variation 
is  much  larger,  and  the  phasing  between  pitch  rate  and  angle 
of  attack  (essentially  the  same  as  body  z-axis  velocity)  is 
much  different.  At  180  kt  (93  m/s),  the  aircraft  is  statical¬ 
ly  stable,  as  evidenced  by  the  presence  of  a  conventional 
short  period  mode.  The  short  period  is  formed  from  the  angle 
of  attack  convergence  and  a  root  resulting  from  the  decomposi¬ 
tion  of  the  third  mode  at  about  158  kt  (31  m/s).  This  decom¬ 
position  also  produces  a  very  slow  speed  divergence  mode,  and 
the  stable  pitch-speed  mode  remains  much  the  same. 

The  lateral  modes  of  motion,  shown  in  Fig.  10,  display 
conventional  Dutch  roll,  roll  and  spiral  modes  for  nozzles 
fully  aft  (6j=0),  with  a  poorly  damped  Dutch  roll  mode  worthy 
of  note.  Near  hover  (10  kt,  5  m/s),  the  aircraft  exhibits  two 
slow  oscillations,  one  stable  and  one  unstable.  The  Dutch 
roll  mode,  which  is  unstable  until  about  120  kt  (62  m/s), 
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Figure  9  Longitudinal  Modes  at  135  kt  (69  m/s) 
and  180  kt  (93  m/s) 
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Figure  10  Lateral  Modes  from  Hover  to  190  kt  (98  m/s) 

increases  in  frequency  with  increasing  V  while  only  slight 
forward  speed  increase  (to  35  kt ,  18  m/s)  splits  the  roll/ 
spiral  oscillation  into  conventional  roll  and  spiral  modes. 

In  the  transition  speed  range,  only  relatively  small  changes 
in  mode  speed  occur  as  the  transition  from  wing-borne  to 
nozzle-borne  flight  is  made.  There  are  some  significant 
changes  in  mode  shape,  however,  as  evidenced  by  Fig.  11.  The 
lateral  mode  shapes  at  160  kt  for  both  nozzle-borne  flight 
(8j=45.2  deg)  and  wing-borne  flight  (0^=0. 0  deg)  are  shown. 

The  value  of  Cn^  (yaw  moment  derivative  with  respect  to  side¬ 
slip)  is  about  half  as  large  at  this  wing-borne  flight  condi¬ 
tion  than  for  the  thrust-borne  case;  hence,  the  lateral  modes 
are  slower  and  exhibit  significantly  larger  amounts  of  sideslip 
relative  to  the  angular  rates.  In  general,  the  transition 
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LATERAL  MODE  SHAPES  AT  160  KT  (82  m/s) 
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Figure  11  Lateral  Mode  Shapes  at  160  kt  (82  m/s), 
a=6  deg,  0j=45.2  deg 


regime  lateral  modes  are  similar  to  those  of  a  conventional 
aircraft,  and  hence  should  not  present  a  major  piloting  problem 
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The  Harrier  exhibits  large  positive  dihedral  effect 
in  the  transition  flight  regime  (Ref.  8),  and  this  results  in 
a  limit  to  allowable  sideslip  excursions.  This  occurs  because 
the  roll  moment  induced  by  sideslip  (L^-Ap)  can  be  larger  than 
the  available  roll  control  power  (Lgrj*6rj  +  ^ga'6a)-  Little 
difficulty  is  exhibitied  at  low  velocities  because  the  aerody¬ 
namic  roil  moment  is  small,  and  no  difficulty  occurs  in  conven 
tional  flight  because  the  dihedral  effect  is  smaller  in 
wing-borne  flight  and  the  ailerons  are  more  effective.  For 
example.  Table  1  gives  values  for  in  both  wing-borne  and 
thrust-borne  flight.  Hence,  the  region  of  concern  is  the 
transition  region.  This  region  is  examined  here  by  the  use  of 
the  static  roll  moment  equation 

V  =  0  =  +  L6rj6rj  +  L6a5a  (1) 

The  amount  sideslip  that  can  be  balanced  by  a  given  amount  of 
control  can  be  calculated  by  Eq.  2 


L*  ■  L* 

ae  --  -  r12  6ri  -  r1 6 
p  rj  Lp  * 


TABLE  1 

DIHEDRAL  EFFECT  AT  160  kt  (82  m/s) 


FLIGHT  CONDITION 

EFFECTIVE  DIHEDRAL,  L„ 

P 

Wing-Borne  a=9.12  deg 

6  .=0.0  deg 

J 

Thrust-Borne  a=6  deg 

0j=45.2  deg 

-461  lb^-ft/deg,  -625  N'm/deg 

-1560  lbf- ft/deg ,  -2115  N-m/deg 

Hence,  the  maximum  allowable  sideslip  can  be  calculated  from 
Eq.  2  by  inserting  the  maximum  control  deflection  in  the 
equation.  Figure  12  shows  the  results  of  this  calculation  for 
the  roll  reaction  jets  and  the  aileron  separately,  and  for 
their  sum.  These  calculations  are  performed  at  6  deg  a;  the 
roll  control  deficiency  is  larger  at  higher  angles  of  attack. 
For  example  at  100  kt  (51  m/s),  the  maximum  sideslip  is  re¬ 
duced  2  deg,  or  13%,  by  increasing  the  angle  of  attack  to  12 
deg.  Also  shown  in  Fig.  12  is  the  sideslip  demonstration 
envelope  required  by  Section  3.3.11  of  Ref.  15.  The  area 
between  the  demonstration  envelope  and  the  sideslip  limited  by 
total  roll  control  application  is  a  region  where  an  uncon¬ 
trollable  roll  is  predicted  by  this  linear  analysis.  A  further 
comment  is  prompted  by  Section  3.3.11.3.1  of  Ref.  15  which  re¬ 
quires  that  for  Level  1  handling  qualities  no  more  than  50%  of 
the  available  roll  control  power  be  used  to  counter  "sideslip 
angles  which  might  be  experienced  in  service  employment". 
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Figure  12  Sideslip  Boundary  Due  to  Roll  Control  Limit 
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This  restriction  would  greatly  expand  the  speed  range  where 
the  maximum  controllable  sideslip  fell  below  the  flight  test 
demonstration  limit. 


This  problem  is  one  that  has  been  recognized  previous¬ 
ly,  as  stated  in  Ref.  8,  and  is  being  corrected  in  the  YAV-8B 
Harrier  by  increasing  the  roll  reaction  jet  effectiveness.  The 
YAV-8B  control  system,  described  in  Ref.  16,  includes  larger 
roll  reaction  jet  thrust  and  a  longer  moment  arm,  as  shown  in 
Table  2.  Also  given  is  the  total  roll  moment  due  to  full 
reaction  jet  thrust.  As  can  be  seen,  this  value  is  72%  larger 
for  the  YAV-8B. 

TABLE  2 

ROLL  REACTION  JET  PARAMETERS -AV-8A  VS  YAV-8B 


Roll  Jet 
Roll  Jet 


"Up"  Thrust 
"Down"  Thrust 


AV-8A _ 

612.  lbf,  2722  N 
237.5  lbf ,  1056  N 


YAV-8B 

750.  lbf,  3336  N 
350.  lbf,  1557  N 


Jet  Moment  Arm 


11.08  ft,  3.377  m 


14.75  ft,  4.496  m 


Total  Roll  Moment 


9412.5  lbf- f t , 
12760.  N-ih 


16225.  lbf- f t , 
22000.  N-tb 


Applying  a  factor  of  1.72  to  the  reaction  jet  curve 
in  Fig.  12  results  in  an  approximation  to  the  sideslip  boundary 
expected  of  the  YAV-8B.  These  results  are  shown  in  Fig.  13. 

The  sideslip  boundary  for  maximum  deflection  of  both  reaction 
jets  and  ailerons  is  well  outside  the  desired  demonstration 
boundary.  At  50  percent  total  control  power,  the  maximum 
allowable  sideslip  boundary,  although  lower  than  the  demon¬ 
stration  boundary  from  67  kt  (34  m/s)  to  about  190  kt  (98  m/s) 
may  correspond  to  sideslip  angles  which  might  be  expected  in 
service  employment.  It  should  be  noted  that  50%  aileron 
deflection  does  not  occur  at  50%  reaction  jet  thrust  in  the 
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Figure  13  Approximate  YAV-8B  Sideslip  Boundary 

YAV-8B,  but  rather  at  about  63%  raction  jet  thrust.  This  is 
due  to  the  reaction  jet  gearing,  which  results  in  full  roll 
reaction  jet  thrust  at  80%  lateral  stick  travel,  rather  than 
100%  as  in  the  AV-8A.  (See  Appendix  C.)  Hence,  the  roll 
control  derivative  (for  less  than  80%  lateral  stick)  is  about 
2.15  times  larger  in  the  YAV-8B  than  in  the  AV-8A. 

This  analysis  indicates  that  the  AV-8B  should  be  able 
to  demonstrate  the  entire  range  of  sideslip  desired  by  Ref.  15, 
and  may  satisfy  Level  1  restrictions  on  positive  effective 
dihedral . 

2.2.4  Influences  of  Ground  Effect  on  Stability 

The  model  for  ground  effect  that  is  used  here,  taken 
from  Ref.  17,  is  described  in  Appendix  C.  Vertical  force, 
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pitch  moment  and  roll  moment  are  all  affected  as  the  aircraft 
descends  from  above  20  feet  (3.1  m)  to  ground  level.  Figure 
14  shows  the  aircraft  mode  shapes  and  speeds  both  out  of 
ground  effect  (OGE)  and  in  ground  effect  (1GE).  As  shown  in 
Section  2.2.2,  there  is  an  unstable  longitudinal  third  mode  as 
well  as  an  unstable  Dutch  roll  mode  at  10  kt  (5  m/s)  OGE.  The 
modes  exhibit  periods  of  25  and  80  sec  respectively.  The  same 
flight  speed  in  ground  effect  results  in  drastically  different 
modes.  The  pitch  effects  stabilize  the  longitudinal  oscilla¬ 
tion,  which  takes  on  some  aspects  of  a  short  period  response, 
especially  in  the  pitch  rate-angle  of  attack  phasing.  The 
thrust  loss  due  to  ground  effect  results  in  the  destabiliza¬ 
tion  of  the  neutrally  stable  (zero-valued  eigenvalue)  altitude 
mode,  and  the  altitude  instability  exhibits  a  time  constant  of 
3.5  sec,  which  corresponds  to  a  time- to-double  of  2.4  sec. 
Relatively  little  pitching  (0.03  deg  pitch  up  per  foot  of  lost 
altitude)  and  speed  change  (-0.05  fps  per  foot  of  lost  altitude) 
occur  as  the  aircraft  drops.  The  rapid  divergence  of  this  mode 
and  the  difficulty  of  quickly  and  accurately  controlling  thrust 
suggest  that  steady  flight  in  ground  effect  should  not  be  at¬ 
tempted,  and  that  ground  effect  altitudes  should  only  be  en¬ 
tered  with  the  intention  of  landing. 

The  changes  in  lateral-directional  mode  shape  and 
speed  that  occur  due  to  ground  effects  are  also  very  severe. 

The  large  destabilizing  roll  moment  effect  splits  up  both  OGE 
lateral  oscillations.  Two  slow  yaw  and  sideslip  modes  are 
formed  (one  stable  and  one  unstable)  and  two  very  fast  roll 
modes  appear.  One  is  very  stable  ( time-to-half-amplitude  is 
0.28  sec)  while  the  other  is  equally  unstable  ( time- to-double- 
amplitude  is  0.29  sec).  This  very  fast  instability  causes 
greatly  increased  pilot  lateral  control  activity  at  low  alti¬ 
tudes,  as  shown  in  Ref.  7,  for  example. 


Changes  Due  to  Ground 


The  linear  model  of  the  V/STOL  aircraft  displays 
characteristics  in  ground  effect  which  agree  with  the  observed 
responses  of  the  aircraft  and  provide  information  about  the 
aircraft  instabilities  in  ground  effect.  An  examination  of 
the  effect  of  ground  effect  instabilities  on  control  system 
design  is  described  in  Chapter  3. 

2.3  STABILITY  VARIATIONS  DUE  TO  ANGLE  OF  ATTACK  AND 
CHANGING  VELOCITY 


2.3.1  Accelerated  Flight  Conditions  at  Various  Angles 
of  Attach 

An  airspeed  of  100  kt  (51. A  m/s)  is  chosen  for  the 
examination  of  AV-8A  stability  variations  with  a  and  acclera- 
tion  along  the  flight  path  (hereafter  referred  to  as  the  deri¬ 
vative  of  velocity  -  V).  At  this  airspeed,  trim  flight  condi¬ 
tions  are  found  at  angles  of  attack  from  0  to  16  deg.  At 
lower  a,  the  lack  of  wing  lift  must  be  made  up  by  vertically- 
oriented  thrust  force,  hence  a  nozzle  angle  to  angle  of  attack 
relationship  results  as  shown  in  Fig.  15.  Quasi-static  trim 
flight  conditions  (conditions  where  the  aircraft  is  trimmed 
about  a  steady  V)  are  also  shown  in  the  figure.  As  can  be 
seen,  the  nozzles  are  deflected  rearward  to  accelerate  and 
forward  to  decelerate.  The  general  reduction  in  nozzle  angle 
at  higher  angles  of  attack  is  due  primarily  to  the  increased 
body  pitch  angle.  The  actual  nozzle  angle  with  respect  to  the 
vertical  changes  only  slightly  as  angle  of  attack  is  varied. 

For  a  given  a,  increasing  or  decreasing  velocity  can 
be  done  essentially  at  constant  engine  throttle  setting. 

Table  3  illustrates  this  observation.  Indeed,  at  16  deg  a, 
these  results  indicate  that  slight  increases  in  throttle 
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Figure  15  No2zle  Angle  vs  Angle  of  Attack 


TABLE  3 


TRIM  ENGINE  RPM  (PERCENT)  AS  A  FUNCTION  OF 
ANGLE  OF  ATTACK  AND  LONGITUDINAL  ACCELERATION 
(at  100  kt  -  51  m/s) 


ACCELERATION 
ANGLE  OF  ATTACK 

-0.2  g ' s 

0.0  g '  s 

0.2  g '  s 

0  deg 

94.3 

94.1 

94.9 

4  deg 

- 

90.6 

; 

8  deg 

86.6 

86.7 

86.5 

12  deg 

- 

81.0 

- 

16  deg 

77.1 

75.8 

80.4 
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are  necessary  to  achieve  either  a  steady  accleration  or  de¬ 
celeration.  Also  noteworthy  in  Table  3  is  the  observation 
that  the  high-a  ,  high-drag  flight  conditions  require  less 
throttle  than  lower  drag  angles  of  attack.  This  is  due  to  the 
increase  in  lift  at  higher  a  which  reduces  the  amount  of 
thrust-borne  weight. 

The  relationship  between  flight  condition  (V,  a,  and 
6 j )  and  amount  of  wing- supported  weight  is  shown  in  Fig.  16. 
The  increase  in  wing- supported  weight  with  a  is  readily  seen, 
as  is  the  resulting  decrease  in  6 ^ .  (Since  8^  is  measured 
relative  to  the  fuselage,  most  of  this  decrease  is  due  to  the 
fuselage  pitch  increase.) 


V  «  +0.2g  V  =  O.Og 


100  KT  (51  m/s)  HORIZONTAL  FLIGHT 
V:  DERIVATIVE  OF  VELOCITY 


o  •  0  deg 


0*8  deg 


NOZZLE  ANGLE,  8j  (deg) 

Figure  16  Wing  Supported  Weight  Variation  with  Angle 
of  Attack 


Nozzle  angle  is  an  important  variable  in  explaining 
the  variation  in  wing-supported  lift  due  to  acceleration  along 


the  flight  path,  since  the  same  curve  seems  to  apply  to  all 
angles  of  attack  between  0  and  16  deg.  Above  about  65  deg 
nozzle  angle  and  below  about  45  deg  nozzle  angle,  moderate 
increases  in  wing-supported  weight  accompany  the  reduction  in 
that  cause  the  aircraft  to  accelerate.  Between  65  and  45 
deg  nozzle  angle,  however,  wing- supported  weight  rises  drama¬ 
tically  a  6j  is  reduced.  This  indicates  a  significant  change 
in  the  character  of  the  flow  field  around  the  aircraft  in  this 
region.  These  curves  also  explain  the  observation  made  con¬ 
cerning  Table  3:  No  additional  thrust  is  needed  to  accelerate 
because  extra  wing  lift  (at  constant  or)  is  available  when  the 
nozzles  are  closer  to  their  aft  position.  Efficiency  of 
flight  (decreased  thrust,  increased  wing  lift)  is  increased  by 
an  increase  in  or  and  a  decrease  in  6 j ,  especially  in  the  45  to 
65  deg  nozzle  angle  range. 

2.3.2  Stability  Variations  with  Angle  of  Attack 

The  basic  modes  of  motion  of  this  vehicle  include  a 
slow,  unstable  longitudinal  oscillation--the  third  mode,  and 
two  real  longitudinal  convergence  modes,  as  well  as  Dutch 
roll,  roll  and  spiral  modes.  The  loci  of  these  modes  as  a 
varies  are  shown  in  Fig.  17.  The  fast  convergences,  angle-of- 
attack  and  roll  modes,  change  somewhat  with  increasing  trim  a 
while  spiral  retains  its  slow,  stable  character. 

Between  12  and  16  deg  o.  the  longitudinal  third  mode 
decomposes  into  two  real  roots.  One  is  a  fast,  unstable  di¬ 
vergence  in  pitch  and  angle  of  attack,  while  the  other  almost 
immediately  combines  with  the  slow  pitch-speed  mode  to  form  a 
conventional  phugoid  mode. 

Dutch  roll  mode  is  unstable  throughout  most  of  this 
angle-of -at  tack  range.  The  actual  crossing  of  the  imaginary 
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Figure  17  Variations  in  Aircraft  Stability 
with  Angle  of  Attack 

axis  occurs  at  about  4  deg  a.  Both  frequency  and  speed  of  di¬ 
vergence  increase  with  a,  making  this  mode  the  most  trouble¬ 
some  to  the  pilot.  At  any  angle  of  attack  above  7  deg,  the 
Dutch  roll  amplitude  doubles  in  less  than  5  sec. 

2.3.3  Stability  Variations  with  Acceleration 

Acceleration  along  the  velocity  vector  (V)  is  varied 
from  +0.2  g  to  -0.2  g  and  its  effect  on  the  stability  of  the 
aircraft  is  described  in  this  section.  Figure  18  illustrates 
these  stability  variations  with  acceleration  at  three  different 
angles  of  attack  (0,  8,  and  16  deg). 
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Figure  18  Stability  Variations  with  Acceleration 

Acceleration  has  an  opposite  effect  on  the  two  longi 
tudinal  real  modes.  The  very  stable  angle-of-attack  mode  is 
made  even  more  stable  by  acceleration,  whereas  the  slow  pitch 
speed  mode  is  made  more  stable  by  deceleration.  A  value  of 
-0.2  g  is  sufficient  to  stabilize  this  mode  even  at  16  deg 
angle  of  attack.  On  the  other  hand,  0.2  g  acceleration  re¬ 
sults  in  an  interaction  with  the  third  mode  which  produces  a 
conventional  phugoid  oscillation  and  a  fast  angle  of  attack 
divergence . 

The  unstable  longitudinal  third  mode  diverges  faster 
under  acceleration  for  0  and  8  deg  angle  of  attack.  This  is 
the  opposite  of  the  effect  of  acceleration  on  the  pitch-speed 
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mode.  At  16  deg  angle  of  attack,  it  appears  that  any  value 
of  acceleration  or  deceleration  will  result  in  an  instabi 1 i ty . 
Deceleration  causes  an  unstable  third  mode,  with  significant 
forward  velocity  oscillation,  whereas  during  acceleration, 
a  fast  unstable  angle-of-attack  divergence  appears.  This 
change  in  appearance  of  the  instability  may  cause  control 
difficulties. 

In  the  lateral -directional  modes,  the  acceleration 
effect  is  much  less.  At  higher  angles  of  attack,  acceleration 
has  a  minor  stabilizing  effect  on  the  unstable  Dutch  roll 
mode,  whereas  the  roll  mode  is  slowed  by  acceleration .  Except 
at  low  angles  of  attack,  spiral  is  affected  little  by  the 
different  acceleration  trim  flight  conditions.  To  summarize, 
acceleration  has  a  secondary  effect  on  the  mode  speeds  rela¬ 
tive  to  the  primary  determinant,  angle  of  attack.  The  effect 
of  acceleration  varies  from  mode  to  mode,  and  often  seems  to 
take  the  form  of  damping  interchange,  where  one  mode  is  made 
more  stable  while  another  is  made  less  stable. 

2.  A  STABILITY  DURING  A  SKI- JUMP  LAUNCH 

2.A.1  Range  of  Flight  Conditions  Examined 

The  ski-jump  launch  of  a  fixed-wing  V/STOL  aircraft 
is  charac terized  by  an  upward  flight  path  angle  (y  =  10  to  20 
deg)  and  a  downward  flight  path  angular  rate  (i.e.,  the  tra¬ 
jectory  is  semi-ballistic).  The  performance  advantage  lies 
in  the  lower  launch  airspeeds  or  higher  launch  gross  weights 
that  are  achievable  (Ref.  9).  It  is  the  goal  of  this  section 
to  examine  the  stability  of  flight  conditions  typical  of 
ski-jump  launches. 
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The  airspeed  at  the  end  of  the  ski-jump  is  taken  as 
65  kt  (33  m/s),  which  is  in  the  lower  range  of  useable  values 
(Ref.  9  lists  65  to  80  kt  -  33  to  41  m/s  -  airspeeds  as  likely). 
At  this  speed,  little  wing  lift  is  gained  by  high  a  ,  so  a 
basic  angle  of  attack  of  6  deg  is  chosen.  Horizontal  flight 
is  examined  first.  Figure  19  details  the  engine  throttle 
setting  required  to  trim  the  aircraft  for  various  values  of  V 
from  0.0  to  0.6  g,  and  for  various  values  of  normal  accelera¬ 
tion.  For  this  velocity,  a  flight  path  rate  (y)  of  5  deg/sec 
implies  a  1.3  g  pull  up,  zero  y  is  1  g  flight,  while  -5  deg/ 
sec  y  gives  a  0.7  g  normal  specific  force  (i.e.,  -0.3  g  down¬ 
ward  acceleration).  The  important  aspect  of  this  plot  is  that 
vertical  specific  force  can  be  traded  for  acceleration  along 
the  flight  path  (V)  with  no  change  in  throttle  setting.  For 
example,  a  0.3  g  decrease  in  vertical  specific  force  allows  a 
0.68  g  increase  in  acceleration  along  the  flight  path  with  no 
increase  in  throttle  setting.  After  the  Harrier  leaves  the  ski- 
jump  with  a  large  upward  velocity  but  downward  acceleration, 
air  speed  increases  rapidly.  Eventually  horizontal  flight 
occurs  as  the  Harrier  reaches  the  proper  velocity  to  begin 
transition  to  wing  borne  flight. 

The  specific  force  polygons  shown  in  Fig.  20  illus¬ 
trate  this  effect.  On  the  left-hand  side  of  the  figure, 
constant  velocity,  1-g  flight  occurs  with  a  thrust  specific 
force  that  is  slightly  larger  than  the  vehicle  weight.  (At 
this  velocity  and  angle  of  attack,  the  aerodynamic  forces  are 
very  small,  and  include  a  slight  negative  lift  force.)  Along 
a  trajectory  that  exhibits  increasing  velocity  (V=0.6  g)  and 
less  than  1-g  vertical  net  specific  force  (Vy=-0.3  g),  less 
thrust  specific  force  (0.91  g)  is  needed  even  though  the  ve¬ 
hicle  is  increasing  speed  rapidly.  This  is  primarily  caused 
by  the  forward  rotation  of  the  thrust  vector  as  the  nozzles 
are  rotated  rearwards.  There  is  also  a  beneficial,  although 
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Figure  19  Thrust  vs  Acceleration  in  the  Ski-Jump  Launch 
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small,  aerodynamic  effect,  since  the  flowfield,  as  modified  by 
the  nozzles,  produces  positive  lift. 

To  provide  vertical  clearance  for  this  semi-ballistic 
trajectory,  the  ski- jump  launch  ramp  produces  a  large  upward 
initial  flight  path  angle,  y*  Increased  aerodynamic  lift, 
caused  by  increased  a  may  be  beneficial  also.  Figure  21  il¬ 
lustrates  the  throttle  changes  necessary  to  accommodate  a  and 
y  variations.  Increased  y  requires  increased  thrust  because 
of  the  work  necessary  to  raise  the  aircraft  mass  in  the  gravity 
field  to  a  higher  potential  energy.  Higher  a  allows  lower 
thrust  due  to  increased  aerodynamic  lift,  and  hence  a  reduced 
amount  of  thrust-supported  weight.  The  specific  force  polygons 
(Fig.  22)  illustrate  these  relationships.  The  increased  lift 
and  drag  at  12  deg  a  are  apparent,  as  is  the  increase  of  thrust 
necessary  to  trim  about  12  deg  y.  This  is  compensated,  to  some 
extent,  by  the  reduction  in  net  vertical  acceleration  from 
0.3  g  downwards  to  0.22  g  downwards. 
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THRUST  vs  ANGLES  OF  ATI  ACK  AND  FLIGHT  PATH  IN  THE 
SKI  JUMP  LAUNCH  (65  K 1  33m/s.  V'0  4g.  Y=  5  deg/sec) 
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Figure  21  Thrust  vs  Angles  of  Attack  and  Flight  Path  in 
the  Ski-Jump  Launch 
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Figure  22  Specific  Force  Polygons 
2.4.2  Stability  Variations  During  the  Ski-Jump  Launch 

The  effects  of  horizontal  and  vertical  acceleration 
on  AV-8A  stability  are  examined  at  constant  angles  of  attack 
and  flight  path,  a= 6  deg,  y=0  deg.  Figure  23  illustrates 
these  variations.  The  effects  on  lateral-directional  modes 
are  not  especially  large,  but  longitudinal  modes  vary  in  a 
pattern  closely  related  to  that  observed  in  Section  2.3. 
Acceleration  causes  the  very  stable  angle  of  attack  mode  to 
become  even  faster,  whereas  the  pitch-speed  mode  slows  and 
becomes  unstable.  The  unstable  third  mode  is  slowed  by  ac¬ 
celeration  and,  for  normal  specific  force  of  0.7  g  (Vy=-0.3 
g) ,  the  third  mode  decomposes  into  a  fast  angle  of  attack 
divergence  and  another  real  which  forms  a  conventional  phu- 
goid  oscillation  as  it  combines  with  the  pitch-speed  root. 
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Figure  23  Stability  Variations  Due  to  Acceleration  During  a 
Ski -Jump  Launch 


The  actual  ski-jump  launch  involves  a  significant 
positive  flight  path  angle,  y>  and  may  involve  higher  angle  of 
attack,  a.  Figure  24  examines  the  effects  of  these  variables 
on  vehicle  stability.  Flight  path  angle  has  negligible  effect 
on  the  lateral  directional  modes  and  on  the  phugoid  mode.  Both 
the  unstable  a  divergence  and  the  fast,  stable,  angle-of -attack 
modes  are  faster  at  y  of  12  deg  than  in  horizontal  flight,  but 
the  effect  is  not  especially  large.  Angle  of  attack,  as  shown 
in  Section  2.3,  does  have  a  very  significant  effect  on  all 
modes.  Its  effect  on  the  a  modes  is  similar  to  but  larger 
than  that  of  flight  path  angle  (on  a  normalized  basis),  and 
higher  a  stabilizes  the  phugoid  mode,  as  expected.  Angle  of 
attack  effects  on  roll  and  spiral  modes  are  small,  but  Dutch 
roll  becomes  significantly  more  unstable  as  a  increases. 
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Figure  24  Stability  Variations  Due  to  Angle  of  Attack 
and  Flight  Path  Angle  During  a  Ski-Jump 
Launch 

Since  angle  of  attack  increase  buys  relatively  little 
performance  increase  at  this  velocity,  as  compared  to  signi¬ 
ficantly  poorer  Dutch  roll  and  a  divergence  characteristics, 
this  analysis  indicates  that  it  would  be  desirable  to  delay 
the  a  increase  until  a  speed  is  reached  where  the  performance 
gain  outweighs  the  stability  difficulty. 


2.5  CHAPTER  SUMMARY 

The  range  of  AV-8A  flight  conditions  from  very  low 
speed  through  transition  to  conventional  flight  is  examined  in 
this  chapter.  The  effects  on  vehicle  stability  and  control  of 
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velocity,  angle  of  attack,  flight  path  angle  and  acceleration 
along  and  normal  to  the  trajectory  are  detailed. 

Longitudinal  Stability  depends  primarily  on  nozzle 
angle,  as  shown  in  Fig.  25,  which  is  a  composite  plot  showing 
the  trend  of  all  results  presented  in  this  chapter.  Aircraft 
velocity  only  has  a  significant  effect  on  the  fast,  stable  a 
mode.  All  other  effects  (angle  of  attack,  flight  path  angle, 
and  acceleration)  primarily  affect  longitudinal  stability 
because  they  require  a  change  in  nozzle  angle  for  trim.  Note 
that  this  figure  illustrates  the  general  trends  with  respect 
to  nozzle  angle;  not  all  the  mode  damping  values  fall  exactly 
on  the  arrows  indicated.  This  chart  suggests  that  those 
flight  conditions  which  require  nozzle  angles  less  than  45  deg 
should  be  avoided. 
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Figure  25  Longitudinal  Mode  Real  Part  vs  Nozzle  Angle 
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Lateral -directional  mode  stability  depends  primarily 
on  velocity  and  angle  of  attack,  with  Dutch  roll  mode  becoming 
unstable  at  higher  angles  of  attack  and  lower  velocity. 

Figure  26  illustrates  these  general  trends.  Other  flight 
condition  variables  (nozzle  angle,  flight  path  angle,  and 
accelerat ion )  have  only  secondary  effects  on  lateral-direc¬ 
tional  stability.  This  chart  suggests  that  low-speed,  high-a 
flight  conditions  be  avoided.  To  some  extent  this  suggestion 
is  in  conflict  with  the  one  given  above,  since  avoidance  of 
high-a  flight  conditions  argues  for  an  early  transition  to 
nozzle-borne  flight  as  speed  decreases,  but  this  will  result 
in  a  nozzle  angle  in  the  undesirable  region. 
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Figure  26  Lateral  Mode  Real  Part  vs 
Angle  of  Attack 
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3. 


DIGITAL  COMMAND  AUGMENTATION  FOR  THE  AV-8A  HARRIER 


A  flight  control  system,  shown  in  Fig.  27,  basically 
consists  of  feedback  control  from  the  aircraft  state  variables 
to  improve  stability  and  feedforward  control  from  the  pilot 
commands  to  achieve  desired  aircraft  trajectories  and  some 
means  of  adapting  this  control  system  to  changes  in  flight 
condition.  Many  state  feedback  and  control  gain  adaptation 
design  procedures  have  been  successfully  developed  and  ana¬ 
lyzed  (Refs.  18  to  21)  but  a  feedforward  control  design  pro¬ 
cedure  in  both  continuous  and  discrete- time  has  not  been  fully 
developed  for  aircraft  applications. 

A  common  procedure  used  in  pilot  feedforward  control 
is  to  mechanically,  and/or  electronically  connect  the  pilot's 
stick,  pedals,  and  throttle  controls  to  the  aerodynamic  surfaces 
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Figure  27  Basic  Control  Structure 
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and  reaction  control  jets  of  the  aircraft.  Gearing,  bob- 
weights,  and  force  producing  devices  are  used  to  provide  the 
proper  pilot  tactual  information  that  enhance  pilot  "feel". 
Pilot  control  of  the  surfaces  and  jets  means  that  the  pilot  in 
many  instances  must  coordinate  control  inputs  to  achieve  de¬ 
sired  aircraft  state  responses.  Often  control  interconnects 
are  employed  to  reduce  pilot  effort  in  turning  flight  and  to 
prevent  departures  at  high  angles  of  attack  due  to  adverse 
yaw. 

The  pilot's  inputs  do  not  have  to  be  directly  con¬ 
nected  to  the  aircraft's  control  inputs.  With  proper  feed¬ 
forward  control-theoretic  formulations,  pilot  command  inputs 
can  be  made  to  be  combinations  of  aircraft  states  and  controls 
Pilots  have  reacted  favorably  to  many  of  the  state  command 
systems  such  as  velocity  command  and  attitude  command  systems, 
Refs.  21  and  24,  and  the  F-16  command  system  described  in  Ref. 
23.  An  advanced  form  of  state  command  systems  is  model  follow 
ing.  Early  efforts  at  having  the  pilot  control  an  ideal  air¬ 
craft  model  then  having  the  actual  aircraft  states  follow  the 
ideal  model  states  are  discussed  in  Refs.  24  and  25.  Essen¬ 
tially  all  linear  state  command  systems  are  a  special  case  of 
the  general  model  following  problem.  Unfortunately,  under  the 
conditions  specified  in  Ref.  24,  it  would  appear  that  perfect 
model  following  is  rarely  possible.  It  is  the  purpose  of  this 
chapter  to  show  that,  under  a  more  suitable  set  of  conditions, 
perfect  model  output  following  i^s  almost  always  possible. 

Model  output  following  can  be  used  in  a  Digital  Flight  Con¬ 
trol  Systems  (DFCS),  and  it  can  be  designed  in  combination 
with  optimal  control  techniques  to  yield  a  desirable  feedback/ 
feedforward  DFCS  for  aircraft.  This  chapter  will  present  the 
results  of  applying  the  theoretical  model  following  tools 
developed  in  Appendix  A  to  the  AV-8A  Harrier. 
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General  CGT  theory  is  employed  to  determine  the  feed¬ 
forward  gains  and  control  structure  so  that  perfect  output 
tracking  is  accomplished.  Non-perfect  conditions  (parameter 
uncertainties,  etc)  occur  in  practice  so  a  proportional - 
integral  (PI)  structure  is  used  in  the  DFCS  design  to  compen¬ 
sate  for  condition  such  as  unmodeled  slowly  varying  distur¬ 
bances  and  modeling  errors. 

The  CGT  is  compatible  with  any  feedback  design  struc¬ 
ture,  such  as  optimal  control  or  eigenvalue/eigenvector  place¬ 
ment.  The  eigenvalue/eigenvector  placement  design,  as  the 
name  implies,  assigns  the  closed-loop  poles  of  the  aircraft  to 
certain  desirable  locations.  In  addition,  it  structurally 
changes  the  closed-loop  eigenvectors  so  that  the  commanded 
output  of  the  aircraft  emulates  the  output  of  the  command 
generator.  For  example,  if  the  command  generator  dynamics 
have  uncoupled  output  dynamics  then  the  eigenvalue/eigenvector 
placement  design  can  force  the  aircraft  to  have  uncoupled  out¬ 
put  dynamics.  An  uncoupled  velocity  command  design  will  be 
demonstrated  using  the  Harrier  and  compared  with  an  optimal 
control  velocity  command  design. 

Both  the  optimal  control  design  and  the  eigenvalue/ 
eigenvector  placement  designs  are  direct  digital  realizations 
which  retain  continuous- time  features.  In  both  cases  the  de¬ 
signer  may  use  continuous-time  criteria  which  is  automati cal ly 
transformed  to  discrete  time  and  solved  as  a  discrete  time 
problem.  The  discrete- time  solution  presented  here  is 
suitable  for  digital  flight  computer  implementation.  For 
evaluation  purposes,  the  discrete-time  closed-loop  system  is 
exactly  transformed  back  to  continuous  time  and  evaluated 
using  continuous-time  specifications  such  as  M1L-F-83300,  Ref. 
15.  It  is  shown  in  Ref.  13  that  the  transformed  discrete  time 
solution  for  the  optimal  control  approach  gives  essentially 
the  same  continuous- time  performance  as  if  the  problem  had 
been  solved  directly  in 


continuous  time.  The  same  is  true  for  the  eigenvalue/eigenvector 
placement  procedure.  Obtaining  equivalent  performance  means  that 
digital  designs  can  be  accomplished  at  lower  sample  rates,  which 
in  turn  reduces  computer  requirements,  and  allows  the  use  of 
more  complex  control  algorithms  in  flight.  Closed-loop  eigen¬ 
values  and  analysis  of  all  control  designs  for  the  Harrier  are 
presented  in  continuous- time  while  control  law  feedback  and 
feedforward  gains  are  presented  in  discrete  time.  The  stair¬ 
cased  output  of  the  control  law  is  smoothed  in  simulation  plots 
to  obtain  a  continuous-time  representation  of  control  movement. 

Once  the  theoretical  tools  for  digital  feedback/feed¬ 
forward  design  are  discussed,  the  next  task  is  to  construct 
the  command  generator  model.  The  AV-8A  Harrier  aircraft  re¬ 
quires  a  considerably  complex  model  since  the  aircraft  transi¬ 
tions  from  a  high  performance  aerodynamic  vehicle  at  high 
speeds  to  an  unstable  reaction  jet  controlled  platform  at  hover. 
One  command  generator  model  cannot  be  used  at  all  flight  condi¬ 
tions  to  always  provide  superior  performance.  The  AV-8A  has 
eight  available  controls  and  can  handle  a  wide  range  of  command 
generators.  Most  of  this  chapter  will  concentrate  on  present¬ 
ing  the  results  from  six  command  generators  used  with  the 
Harrier  at  different  flight  conditions.  Comparisons  will  then 
be  made  and  the  chapter  is  concluded  with  a  discussion  of  the 
command  generators  and  feedback  designs. 


3.1  COMMAND  GENERATOR  TRACKER  THEORY -DIGITAL  FLIGHT 
CONTROL  SYSTEM  DESIGN 

3.1.1  Explicit  Model  Following 

The  design  of  the  command  generator  tracker  propor¬ 
tional-integral  (CGT-PI)  control  law  is  based  upon  the  defini 
tion  of  a  coupled  linear  time- invariant  model  of  the  AV-8A 
described  as 


Ax ( t )  =  F  Ax ( t )  +  G  Au  ( t ) 


where  Au(t)  is  an  m-vector  representing  control  perturbations 
and  Ax ( t )  is  an  n-vector  representing  the  aircraft's  dynamic 
states.  The  purpose  of  the  CGT-P1  control  vector  is  to  sta¬ 
bilize  the  aircraft  about  a  trajectory  where  the  desired 
output  of  the  states  and  controls, 


Ay(t)  =  H  Ax ( t )  +  D  Au(t) 

perfectly  tracks  the  output  of  a  linear,  time- invariant  com¬ 
mand  generator, 


FA(t) 


G  Au  ( t ) 
m  — m 


=  VS»(t>  +  W1) 

The  m  -vector  Au  (t)  represents  the  pilot's  input  to  the 
m  — m 

command  generator  and  Ax  (t)  represents  the  n  -vector  of  the 
b  —m  K  m 

command  generators  dynamic  states.  The  aircraft  output, 

Av(t),  and  the  command  generator  output,  Aym( t ) ,  are  £-vector 
and  there  are  at  least  as  many  controls  as  there  are  commanded 
outputs.  The  trajectory  the  aircraft  states  and  controls  fol¬ 
low  when  perfect  tracking  occurs  is  the  "star"  trajectory  and 
has  the  property,  starting  from  the  initial  time,  tQ, 

i . 

HAx" ( t  )  +  DAu'  ( t )  =  HmAxm(t)  +  DmAum(t) 

The  first  simplfying  assumption  used  to  obtain  a  solu¬ 
tion  for  the  CGT-P1  is  that  Au  (t)  is  constant  before  and  after 

— m 

the  initial  time.  The  plant  is  assumed  to  be  tracking  the  com¬ 
mand  generator  previous  to  t^  and  at  t^,  Au^  steps  and  remains 
constant.  The  control  must  transfer  the  aircraft  from  the  old 
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to  the  new  star  trajectory  caused  by  stepping  Au^  while  penal¬ 
izing  the  deviations 

„  *  *  ^  .  .  * 

Ax  =  Ax  -  Ax  ;  Au  =  Au  -  Au  ;  Av  =  Au  -  Au 

in  the  following  manner 


The  control  rate  deviative  is  included  in  the  cost  function 
in  order  to  obtain  the  proportional - integral  property  in  the 
control  law  as  shown  in  Ref.  18. 

The  command  generator  tracker  for  arbitrary  Aum(t)  is 
presented  in  Ref.  12  for  the  continuous- time  case  and  in  Ref. 
11  for  the  discrete- time  case.  The  solution  requires  the  use 
of  compensators  whose  poles  are  the  transmission  zeroes  of  the 
plant.  In  addition,  up  to  n  Aum(t)  derivatives  may  be  needed 
for  the  continuous- time  case  and  up  to  n  future  values  of 
A-m  k  may  needed  in  the  discrete-time  case.  The  complexity 
of  arbitrary  Aum(t)  inputs  exceeds  the  scope  of  this  work.  On 
the  other  hand,  Ref.  1A  demonstrates  that  the  simplifying 
assumption  for  Auffl  presented  in  the  previous  paragraph  per¬ 
forms  very  well  if  Aum  is  slow’ly  varying  when  compared  to  the 
closed-loop  plant  time  constants. 

A  digital  control  law  which  minimizes  the  cost  func¬ 
tion  can  be  obtained  by  restricting  the  control,  Au ,  to  change 
only  at  equally  spaced  sampling  intervals,  At.  The  problem 
reduces  to  the  sampled-data  regulator  solved  in  Ref.  27.  The 
basic  equations  involved  are  listed  in  Table  A.  The  aircraft 
model  error  dynamics  are  obtained  by  subtracting  the  dynamic 
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TABLE  4 

TRANSFORMATION  OF  CONTINUOUS -TIME  OPTIMAL  CONTROL 
PROBLEM  TO  DISCRETE-TIME 


:0.\'T  I  .VIOL’S  -  TIME  AIRCRAFT  MODEL  ERROR  DYNAMICS  WITH  COMPENSATOR 
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equation  for  the  star  trajectory  shown  in  Eq .  A-l  (no  distur¬ 
bance)  from  the  aircraft  dynamics.  The  star  trajectory  is 
linearly  related  to  the  command  generator  states  and  controls 
as  shown  in  Table  5.  The  feedforward  matrix  linearly  relating 
the  command  generator  states  to  the  aircraft  states  satisfies 
the  algebraic  matrix  equation  shown  in  Table  5.  A  discussion 
and  techniques  for  solving  the  matrix  equation  are  presented 
in  Appendix  A.  Special  attention  is  given  to  the  importance 
of  transmission  zeroes  in  finding  a  solution.  The  discrete¬ 
time  control  law  solution  is  presented  in  Table  6  for  the  lin¬ 
earized  dynamics.  The  closed-loop  discrete- time  system  is  ob¬ 
tained  by  converting  the  PI  control  law  in  incremental  form 
back  to  the  position  form  as  shown  in  Ref.  26.  Eigenvalue 
evaluation  of  the  control  law  is  obtained  by  mapping  theclosed- 
loop  discrete- time  matrix  back  to  an  equivalent  continuous-time 
plant  using  the  natural  logarithm  of  a  matrix.  Only  eigen¬ 
values  in  the  primary  band  in  the  s-domain  are  computed.  High 
frequency  states  such  as  body  bending  modes  are  not  considered 
in  this  analysis,  hence  frequency  folding  problems  do  not  occur 
as  we  transform  back  and  forth  between  continuous  and  discrete 
time . 


The  CGT-PI  control  in  incremental  form  that  would  be 
implemented  onboard  the  aircraft  flight  computer  has  the  form 


-k  '  — k-1  ‘  Cl(-k'-k-l)  ‘  AtC2(^k-l‘^m,k-l) 


1 A21  +  C1A11]  (-m,k  '  -m,k-l} 


(3) 


and  uses  total  variables.  When  perfect  tracking  occurs  the 
command  error  y,  , -y  ,  is  zero,  the  control  contribution 

j.  j,  ^k  - 1  ■2-m ,  k 

•Cl(£k‘*k-1>  is  cancelled  clAll(Vk'S«,k-l)  and 

A~, (x  , -x  ,  , )  is  the  feedforward  term  to  add  to  u/  ,  to 

21  — m,k  — m,k-l  -k-1 

ensure  u^  maintains  the  tracking.  As  discussed  in  Appendix  A, 
the  command  generators  are  chosen  so  that  either  the  pilot's 
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TRANSFORMATION  OF  COMMAND  GENERATOR  TRACKER 
TO  DISCRETE-TIME 
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u  ,  inputs  eventually  must  return  to  zero  making  x  , -x  .  , 
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zero  (example:  acceleration  commands)  or  the  matrices  A22  and 
A|2  are  sparse  with  entries  occuring  because  of  geometrical 
considerations.  For  implementation  C^ ,  C2 ,  and  any  entries 
occuring  in  A22  and  A^^  which  cannot  be  reconciled  based  on 
the  physical  geometry  of  the  problem  would  be  gain  scheduled 
with  flight  conditions  as  discussed  in  Ref.  26.  It  is  possible 
to  eliminate  the  term  in  the  control  law  which  uses  A22  and 
A^2  using  a  discrete- time  version  of  Davison's  Servocompen- 
sator  shown  in  Ref.  23,  and  an  example  is  given  in  Appendix  A. 

The  prime  on  u^_2  indicates  that  the  control  command  should  be 
limited  to  the  natural  control  limits  of  the  aircraft  in  order  to 
prevent  windup  as  discussed  in  Ref.  1U  .  A  block  diagram  of  the 
control  law  is  shown  in  Fig.  28  for  linear  command  generators. 
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DISCRETE-TIME  RICCATI  EQUATION 
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TYPE  0  DISCRETE-TIME  OPTIMAL  CONTROL  LAW  GAINS 
(K:  K2J  =  (r  *  r^prD)‘1  <r^P*D  +  mt) 


TYPE  1  DISCRETE-TIME  OPTIMAL  CONTROL  LAW  GAINS 
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PI  CONTROL  LAW  IN  INCREMENTAL  FORM 


ERROR  DYNAMICS: 
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CLOSED-LOOP  DISCRETE-TIME  SYSTEM 
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EQUIVALENT  CLOSED- LOOP  CONTINUOUS-TIME  SYSTEM 
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3.1.2  Implicit  Model  Following 

In  this  section  an  alternate  procedure  is  presented 
for  determining  the  stabilization  gains  in  the  CGT-P1  control 
law.  The  procedure  provides  a  technique  for  choosing  gains 
which  cause  the  closed-loop  aircraft  dynamics  to  have  the  same 
eigenvalues  as  a  desirable  aircraft  model  and  restructures  the 
eigenvectors  of  the  closed-loop  aircraft  dynamics  so  that 
output  transients  of  the  aircraft  are  the  same  as  the  desir¬ 
able  aircraft  model.  The  procedure  is  the  solution  to  the 
implicit  model  following  problem  and  is  accomplished  using  the 
feedforward  matrices  in  the  CGT  under  special  conditions.  A 


discussion  and  derivation  of  the  eigenvalue/eigenvector  place¬ 
ment  procedure  is  given  in  Appendix  A. 4.  This  procedure  is 
used  in  one  of  the  AV-8A  control  designs  in  Section  3. A. 

The  equations  for  solving  the  eigenvalue/eigenvector 
placement  (EEP)  problem  for  the  CGT-PI  control  law  are  out¬ 
lined  in  Table  7.  The  same  aircraft  error  dynamics  as  in  the 
optimal  control  problem  are  used.  The  EEP  gains  are  to  sta¬ 
bilize  the  plant  about  the  star  trajectory.  It  is  important 
to  realize  that  the  command  generator  producing  the  star  tra¬ 
jectory  is  usually  not  the  same  as  the  ideal  aircraft  model 
used  in  the  EEP  problem.  The  ideal  aircraft  model  in  Table  7 
represents  the  designers  choice  of  how  the  closed-loop  air¬ 
craft  and  compensator  transient  dynamics  should  behave.  The 
designer  can  almost  always  only  match  m  Ay  transients  to  m  Ay 
desirable  transients,  where  m  is  the  number  of  aircraft  con¬ 
trols.  Insufficient  control  power  will  reduce  the  number  of 
matching  transients  even  further.  Once  the  EEP  gains  are 
determined,  the  control  law  is  the  same  as  Eq .  3  and  can  be 
evaluated  using  the  same  equivalent  closed-loop  continuous¬ 
time  system  shown  in  Table  6. 

Both  design  procedures  shown  in  Section  3.1.1  and  in 
this  section  are  direct  digital  designs.  In  optimal  control, 
performance  is  changed  by  adjusting  the  continuous- time  weight¬ 
ing  matrices  Q  and  R  independent  of  the  sampl ing- time .  In 
EEP,  performance  is  changed  by  adjusting  the  elements  in  , 

G^,  Tx  and  Tu  also  independent  of  the  sampl ing- time .  A  com¬ 
parison  between  an  optimal  control  design  and  the  EEP  design 
is  presented  in  Section  3. 3. A. 
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TABLE  7 


SOLUTION  OF  THE  DISCRETE-TIME 
EIGENVALUE/EIGENVECTOR  PROBLEM  t-2904 
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FEEDFORWARD  MATRIX  EQUATION  FOR  CONTROL  GAINS 
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3.2 


COMMAND  ANALYSIS  FOR  THE  AV-8A  AIRCRAFT 


The  AV-8A  aircraft  has  ten  control  surfaces,  reaction 
jets,  and  engine  controls  which  are  then  reduced  to  eight  con¬ 
trols  in  Appendix  C.  The  eight  controls  which  can  be  com¬ 
manded  by  the  DFCS  are:  RPM,  AN^  ;  stabilator,  A6g;  nozzle 
angle,  A0 ^ ;  pitch  jets,  A6pj  ;  for  the  longitudinal  axis  and 

aileron,  A6  ;  roll  jets,  A6  rudder,  A6  ;  and  yaw  jet,  A6 

a  r j  r  yj 

for  the  lateral-directional  axis.  Theoretically  it  should  be 
possible  to  command  eight  independent  combinations  of  states 
and  controls  with  eight  independent  controls.  Practically  for 
aircraft,  the  control  effect  matrix,  G,  has  at  most  rank  six, 
hence  sufficient  control  power  is  available  for  commanding  six 
combinations  of  states  and  control.  The  values  in  the  feed¬ 
forward  matrices  S2^  and  S22  indicate  how  much  travel  is 
needed  by  the  controls  to  produce  the  desired  response.  The 
more  control  power  available,  the  smaller  the  numbers  are  in 
S21  and  S22  relative  to  the  maximum  control  throw.  An  example 
of  values  in  S2^  and  S22  using  various  mixes  of  controls  and 
commands  is  presented  in  the  next  three  sections.  In  all 
cases  except  for  very  low  speed  flight,  the  command  systems 
for  the  AV-8A  designs  in  this  chapter  will  have  fewer  commands 
than  controls  because  of  insufficient  control  power. 

The  V/STOL  aircraft  undergoes  substantial  changes  in 

characteristics  as  the  vehicle  transitions  from  high  speed  to 

hover.  More  than  one  command  system  is  needed  to  optimize 

piloting  tasks  depending  on  the  flight  condition.  Examples  of 

using  more  than  one  command  system  onboard  an  aircraft  with 

the  pilot  either  choosing  the  appropriate  one  or  the  computer 

phasing  the  systems  in  and  out  are  shown  in  Refs.  23,  29,  and 

30.  In  this  study,  six  different  command  and  control  systems 

shown  in  Table  8  are  constructed,  designed,  and  simulated. 

Each  command  model  is  determined  by  the  desired  pilot  commands, 

u  ,  the  desired  model  output  vector,  y  ,  and  the  available 
— m  mii 


TABLE  8 

COMMANDS,  STATES,  AND  CONTROLS  FOR  THE  AV-8A 


TERMINAL  AREA  ATTITUDE 
COMMAND  RAMP/STEP  MODEL 


controls,  u.  All  command  models  have  five  pilot  inputs.  Two 

inputs  for  lateral  and  longtitudinal  stick,  one  input  for  the 

other  hand  (usually  consisting  of  a  speed  command),  one  input 

for  the  pedals  and  the  fifth  input  can  be  a  thumb  button  for  a 

rate  command  (or  a  thumb  wheel  for  a  position  command).  The 

rate  commands  are  model  inputs,  uffl ,  which  are  integrated  by 

the  model  to  generate  a  state,  x  ,  that  is  used  in  y  .  For 
°  —  m  ■Hu 

example,  a  step  in  the  vertical  velocity  command,  V  ,  causes 

z  ,  c 

the  position  command,  z  ,  in  the  model  to  be  a  ramp.  A  posi¬ 
tion  command  is  a  model  input  that  is  directly  used  in  such 
as  the  pitch  angle  command.  An  example  of  the  model  state 
equations  using  the  enroute  velocity  command  system  for  demon¬ 
stration  purposes  is  shown  in  Eqs.  B-l  and  B-2. 

In  the  first  four  command  models  in  Table  8  there  are 
more  controls  than  commands.  There  are  a  number  of  procedures 
for  designing  the  control  system  with  excess  controls.  A 
practical  but  conservative  approach  is  to  interconnect  two 
controls  to  form  one  control  and  continue  this  process  until 
the  commands  and  controls  are  equal,  as  is  done  in  the  current 
AV-8A  control  system.  The  interconnect  can  be  chosen  so  that 
the  controls  saturate  at  the  same  time.  An  alternate  approach 
is  to  use  a  weighted  pseudoinverse  to  find  the  feedforward 
matrices  in  Eq .  A-7  and  the  feedback  matrices  for  CGT-P1  in 
Eq .  A-28.  The  pseudoinverse  optimally  decides  how  to  use  the 
excess  controls  to  accommodate  the  commands.  A  pseudoinverse 
design  discussing  its  many  features  is  presented  in  Refs.  10 
and  14  for  a  fighter  aircraft.  More  implications  about  the 
pseudoinverse  design  which  is  the  method  used  for  the  AV-8A 
is  addressed  in  Section  3.3. 

The  V/STOL  flight  regimes  as  an  aircraft  approaches 
a  ship  to  land  can  be  considered  to  have  three  regions;  high 
speed  flight,  transition,  and  low  speed  flight  to  hover.  An 
example  of  these  regions  in  the  AV-8A  approach  profile,  taken 
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from  Ref.  7  is  shown  in  Fig.  29.  The  rest  of  this  section 
discusses  the  command  models  most  applicable  to  each  regime. 


3.2.1  High  Speed  Flight 

The  AV-8A  at  high  speeds  is  commanded  by  the  pilot 
using  the  air  combat  maneuvering  (ACM)  ramp/step  model  in 
Table  8.  The  high  speed  flight  condition  is  chosen  at  the 
extreme  of  the  available  AV-8A  flight  data,  87.6  m/s  (170  kt ) 
straight-and- level  flight.  The  pilot  command  generator  inputs 
are  vehicle  velocity  rate,  V,  normal  acceleration,  an ,  angle 
of  attack,  a,  roll  rate  about  the  velocity  vector,  pw,  and 
sideslip,  p.  The  command  generator  states  are  vehicle  ve¬ 
locity,  V,  flight  path  angle,  y,  and  wind-axis  roll  angle, 

0y.  The  ACM  ramp/step  model  is  the  most  complex  of  the  com¬ 
mand  generator  models  and  is  fully  described  in  Appendix  B.l. 
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Figure  29 


AV-8A  Approach  Profile 


The  ACM  ramp/step  model  is  a  generalization  using  the 
command  generator  tracker  of  the  ACM  command  vector  used  in 
Ref.  14.  The  design  in  Ref.  14  required  that  the  Euler  angles, 
0  and  6 ,  be  removed  from  the  aircraft  state  vector  for  design 
because  the  angles  do  not  approach  a  constant  value  when  a 
rate  is  commanded.  The  CGT  eliminates  this  difficulty  by 
using  command  generator  states,  y  and  ,  which  are  also 

not  constant  for  a  rate  command  and  are  related  linearly  to  6 
and  o.  The  ACM  ramp/step  design  greatly  improves  the  perfor¬ 
mance  of  outer-loop  tracking,  which  tended  to  be  inferior  in 
Ref.  14. 


Originally,  the  ACM  model  included  lateral  accelera 
tion,  Aa^7  m,  the  velocity  heading  angle,  A£m,  and  eight  con¬ 
trols.  The  resulting  feedforward  matrices,  S22  and  S22  (in 
English  units)  then  become 
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The  longi t ud i na 1  -  1  a t e ral -di rec t ional  coupling  evident  in  Eq .  A 

is  caused  by  the  coupling  in  the  control  effect  matrix,  f,  due 

to  the  reaction  jets.  For  minor  Aa  commands  all  four  lat- 

y  ,m 

eral  directional  controls  saturate  as  shown  in  the  last  column 
of  the  matrix  A99  in  Eq .  A.  Thus,  the  AV-8A  does  not  have 
adequate  side-force  control  power.  The  feedforward  matrices 
are  recomputed  by  eliminating  a^  m  and  from  the  model  and 
6  ^  j  and  6^  front  the  controls.  This  results  in 


0.141  0.782  -0.00338 
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The  values  in  Eq .  5  indicate  a  sufficient  operating  range 
before  control  saturation.  Futhermore,  the  control  operating 
range  is  not  significantly  decreased  by  removing  6  ^  and  Sr- 
as  active  controls.  The  yaw  jet.  6  .  is  included  because 

yj 

improved  directional  stability  at  high  angles  of  attack  could 
result  by  using  a  control  whose  control  power  does  not  depend 
on  the  fuselage  blocked  air  stream. 


3.2.2  Transition 


Three  command  generators  are  used  at  two  transition 
flight  conditions;  the  coordinated  turn  dynamic  model,  the 
coordinated  turn  ramp/step  model,  and  the  enroute  velocity 
ramp/step  model.  The  two  transition  flight  conditions  are  at 
51.5  m/s  (100  kt)  and  18  m/s  (35  kt ) .  The  dynamic  model  is 
included  in  the  command  generators  to  compare  dynamic  model 
following  with  ramp/step  model  following.  A  complete  de¬ 
scription  of  the  dynamic  model  is  given  in  Appendix  B.  Dy¬ 
namic  models  constitute  the  majority  of  modern  control  model 
following  results  as  evidenced  in  Refs.  25  and  31  to  33. 
Ramp/step  models  have  previously  been  ignored  in  modern  control 
analysis  (the  required  theory  has  only  recently  become  avail¬ 
able)  but  ramp/step  models  are  used  in  many  classical  designs 
(e.g.  ,  Ref.  22) . 


The  coordinated  turn  and  enroute  velocity  ramp/step 
models  differ  in  that  sideslip,  p,  is  tracked  in  the  former 
while  cross-range  velocity,  ,  is  tracked  in  the  later 
model.  The  downrange  and  crossrange  velocities  are  obtained 
by  rotating  the  inertial  velocities  of  the  vehicle  by  the  com¬ 
manded  body  heading  angle  .  In  a  steady-state  turn  without 
wind  and  commanded  to  zero,  sideslip  is  also  zero.  The 

enroute  velocity  ramp/step  system  is  similiar  to  the  flight 
tested  velocity  command  system  in  Ref.  26. 


The  coordinated  turn  ramp/step  and  dynamic  models 
have  the  same  command  vector,  y,  but  markedly  different  feed¬ 
forward  matrices.  Of  particular  interest  are  the  feedforward 
matrices  and  S 2 ^  which  must  be  gain  scheduled  for  CGT-PI 

implementation.  The  next  four  equations  show  and  S 2 2  f°r 

the  dynamic  and  ramp/step  models,  respectively,  at  51.5  m/s 
(100  kt).  The  dynamic  model  requires  at  least  65  nonzero 
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elements  in  and  $2^  while  the  ramp/step  model  requires 

only  18  nonzero  elements  to  perform  the  same  function.  Further 
comparisons  between  the  two  command  generators  are  presented 
in  Section  3.3.3. 


3.2.3  Low  Speed  Flight 

The  most  difficult  regime  to  control  the  AV-8A  is 
near  hover  where  precise  control  is  needed  to  land  the  air¬ 
craft.  Two  command  generators  are  used  for  low  speed  flight; 
a  velocity  ramp/step  model  and  an  attitude  ramp/step  model. 
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Both  command  systems  require  precise  position  information 
which  is  usually  available  from  a  radar  system  in  the  terminal 
area.  The  velocity  ramp/step  model  uses  the  pilot's  local- 
level  velocity  commands  to  formulate  a  desirable  position 
trajectory  in  inertial  space  which  the  aircraft  tracks.  Any 
lateral  drift  is  compensated  for  by  rolling  the  aircraft. 

Many  V/STOL  pilots  do  not  care  to  roll  the  aircraft  near  hover 
since  this  reduces  the  jet  thrust  supporting  the  aircraft. 
Further  difficulty  in  roll  control  is  caused  by  the  ground 
effect  (see  Appendix  C)  which  destabilizes  roll  as  the  air¬ 
craft  approaches  the  ground. 

The  attitude  command  replaces  the  cross-range  ve¬ 
locity  pilot  command  with  a  direct  roll  angle  command.  By 
fixing  the  pitch  and  roll  commands,  the  pilot  maintains  a 
level  platform  at  the  cost  of  lateral  drift.  The  attitude  and 
velocity  command  systems  are  designed  at  5.15  m/s  (10  kt )  using 
the  three  reaction  jets,  nozzle  angle,  and  RPM  to  uniquely 
accommodate  the  five  commands.  Analysis  of  AV-8A  models  at 
speeds  lower  than  10  kts  indicated  the  data  set  did  not  give 
results  which  agreed  with  low  speed  Harrier  behavior  (i.e.,  a 
highly  unstable  transmission  zero  developed  in  the  lateral 
dynamics ) . 


3.3  V/STOL  DIGITAL  FLIGHT  CONTROL  SYSTEM  RESULTS 
3.3.1  Control  Design  Procedures 

The  design  procedure  for  the  linear-optimal  CGT-P] 
design  involves  iterating  on  the  elements  in  the  continuous¬ 
time  cost  function  weighting  matrices  until  the  system  re¬ 
sponses  and  control  motions  have  desirable  transients.  In 
addition,  closed-loop  mapped  discrete-time  eigenvalues  should 
be  located  in  desirable  regions  in  the  left-half  complex 
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plane.  The  elements  in  Q  and  R  are  specified  as  the  inverse 
of  a  maximum  mean-square  value  as 


Q  is  positive  definite  with  cross  weighting  between  states  and 
controls,  caused  primarily  by  acceleration  weights,  while  R  is 
diagonal  and  also  positive  definite. 

Experience  has  shown  that  one  convenient  method  for 
choosing  Q  and  R  is  to  specify  Au  and  Au  as  fractions  of 
maximum  control  positions  and  rates,  then  perform  all  further 
iterations  on  Ax  ,  Ax  ,  and  Aym„. .  Each  command  response 
and  each  eigenvalue  is  primarily  affected  by  a  few  of  the 
remaining  state  weights  which  can  be  selectively  adjusted. 

The  closed-loop  eigenvectors  are  calculated  to  observe  which 
states  are  predominate  in  a  given  eigenvalue.  The  last  fea¬ 
ture  has  a  significant  effect  on  the  presentation  format  of 
the  closed-loop  eigenvalues  for  the  AV-SA.  Typical  aircraft 
have  readily  identifiable  modes  such  as  short  period.  Dutch 
roll,  spiral,  etc.  The  AV-8A  does  not  always  exhibit  such 
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behavior,  and  rather  than  create  new  names,  the  eigenvalues 
will  be  identified  by  the  state(s)  most  predominate  in  the 
associated  eigenvector. 

The  Q  and  R  matrices  at  the  four  design  flight  con¬ 
ditions  are  shown  in  Table  9.  The  weightings  remain  almost 
constant  except  at  low  speed  flight  where  the  most  aerodynamic 
variation  occurs.  M1L-F-83300,  the  flying  qualities  specifi¬ 
cation  for  piloted  V/STOL  aircraft,  Ref.  15,  is  used  as  an  aid 
in  choosing  Q  and  R.  The  CGT-PI  design  employs  full  state 
feedback  without  considering  the  effect  of  estimators,  hence 
the  closed-loop  eigenvalue  specifications  are  easily  met  for 
Level  1,  Category  A  and  C  flight'.  The  parts  of  MIL-F-83300 
that  deal  with  control  power  and  V/STOL  response  proved  to  be 
more  useful  in  the  eigenvalue/eigenvector  placement  design 
procedure  and  the  choice  of  the  coordinated  turn  dynamic  model 
shown  in  Appendix  B. 

The  sampling  rate  chosen  for  the  AV-8A  is  10  frames 
per  second.  A  number  of  investigations,  such  as  Refs.  19, 

26,  and  35,  have  shown  that  direct  digital  designs  can  use  low 
sample  rates  without  degrading  system  performance.  Reference 
35,  in  particular,  shows  that  for  aircraft  systems,  10  frames 
per  second  seems  to  be  the  break  point  between  adequate  and 
poor  performance.  The  digital  design  for  a  helicopter  in 
Ref.  26  was  successfully  flight  tested  using  10  frames  per 
second  with  a  control  law  very  similiar  to  the  CGT-PI .  The 


Control  design  without  considering  estimator  effects  is 
partially  justified  by  the  separation  theorem  but  recent 
results  in  Ref.  34  indicates  that  there  may  not  be  key 
guaranteed  gain  and  phase  margins  using  optimal  control 
when  estimators  are  employed.  A  useful  area  of  future 
research  would  be  to  compute  CGT-PI  gain  and  phase  margins 
with  estimator  effects. 
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complex  P1F  control  law  in  Ref.  26  required  60%  of  a  computa¬ 
tion  frame  to  perform  all  control  calculations  indicating  that 
onboard  CGT-P1  computations  can  be  done  within  a  computation 
frame  if  it  operates  at  10  frames  per  second. 


3.3.2  High  Speed  Flight 


Vectored  thrust,  A0  .  ,  and  reaction  control,  A6  .,  are 

j  yj 

the  two  unusual  features  used  with  the  ACM  command  model  in 


high  speed  flight.  The  reaction  control  is  expected  to  im¬ 
prove  sideslip  response  while  vectored  thrust  will  improve 
normal  acceleration  response  with  minimal  angle  of  attack 
changes.  A  complication  occurs  when  using  vectored  thrust 
because  0j  is  at  the  position  limit  (0^=0°)  for  straight  and 
level  flight.  This  section  discusses  how  CGT-PI  effectively 
accommodates  the  high  speed  design  complexities. 


The  open- loop  and  mapped  closed- loop  eigenvalues  for 
the  AV-8A  at  high  speed  is  shown  in  Table  10.  Spiral  and  a  Aq 
mode  are  open-loop  unstable  and  interact  with  the  integrator 
compensator  states  to  form  stable  closed-loop  complex  pairs. 
The  primary  effect  of  using  the  yaw  nozzle  is  to  further  sta- 
blize  the  sideslip  mode  (Av)  as  expected.  Figure  30  shows  the 
effect  of  the  reduced  Av  time  constant  for  a  Ap  step  com- 
mand.  Comparing  the  control  responses  in  Figs.  30c  and  d,  the 
CGT-PI  control  law  hardly  changes  aileron  authority  while 
transferring  rudder  authority  to  the  yaw  nozzle.  Note  that  in 
Fig.  30c,  where  the  controls  and  commands  are  equal,  the  con¬ 
trols  have  an  initial  transient  which  optimally  transfers  the 
system  states  to  the  new  *  trajectory  by  locking  on  to  the 
control  *  trajectories.  In  Fig.  30d  there  is  extra  freedom 
caused  by  the  extra  reaction  jet  control  and  the  controls  do 
not  necessarily  lock  on  to  the  control  *  trajectory  indicated 
by  the  pseudoinverse  method  for  finding  the  feedforward  matri¬ 
ces.  The  states,  however,  always  lock  on  to  the  *  trajectory. 
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TABLE  10 

EFFECTS  OF  CGT-P1  ON  THE  DYNAMIC  MODES  AT  87.5  m/s  (170  kt ) 


T-290? 


DYNAMIC 

MODE 

OPEN -LOOP 

CHARACTERISTICS 

CLOSED- LOOP 
(46  s  =  2 

WITHOUT 

.0,  A6ei  = 

YAW  NOZZLE 
3.0> 

CLOSED- 

(A6S 

LOOP  WITH  YAW 

=3.0,  46ej=l. 

NOZZLE 

0) 

NATURAL 
FREQUENCY , 
rad/sec 

DAMPING 
RATI  0 

TIME 

CONSTANT , 
sec 

NATURAL 
FREQUENCY . 
rad/sec 

damping 

ratio 

TIME 

CONSTANT , 
sec 

NATURAL 

FREQUENCY 

rad/sec 

DAMPING 

1  RATIO 

TIMK 

CONSTANT , 
sec 

(ap-.lAPw 

- 

- 

- 

2.33 

0 . 596 

- 

2.55 

0.637 

- 

AO  -  Aq 

- 

- 

- 

1.44 

0.597 

- 

0.99 

0.764 

- 

Ar-A<t> 

2.47 

0.252 

- 

2.76 

0.363 

- 

3.08 

0.370 

- 

Au- JAV 

- 

- 

- 

0.333 

0.837 

- 

0.333 

0.839 

- 

Aw 

- 

- 

0.348 

- 

- 

0.239 

- 

- 

0.366 

ROLL  MODE 

- 

- 

0.970 

- 

- 

0.408 

- 

- 

0.387 

fa.. 

- 

- 

_ 

_ 

0.602 

1.81 

0.91 

- 

* 

* 

- 

0.854 

Av 

- 

- 

- 

- 

- 

2.42 

- 

0.779 

SPIRAL 

- 

- 

-0.372 

- 

- 

- 

- 

- 

- 

FIHICOID 

0  199 

0.481 

- 

- 

- 

- 

- 

- 

- 

Aq 

- 

* 

-0.767 

- 

- 

- 

- 

- 

- 

Any  mismatch  in  one  control  such  as  A6  in  Fig.  30d  must  be 
compensated  for  by  a  (barely  perceptable)  mismatch  in  another 
control,  Adiyj  in  Fig.  30d.  The  control  trajectory  freedom  in 
the  pseudoinverse  design  is  a  characteristic  of  CGT-PI  (but 
not  a  characteristic  of  the  Type  0  CGT,  as  shown  in  Ref.  10) 
and  no  disadvantage  with  the  method  has  been  encountered  even 
in  the  full  nonlinear  evaluation,  performed  in  Ref.  14,  for  a 
fighter  aircraft.  In  fact,  many  advantages  are  available 
especially  when  the  extra  controls  favorably  saturate  for 
maximum  performance  when  large  commands  are  requested,  as 
discussed  in  Ref.  13. 

Recall  from  Fig.  28  that  CGT-PI  handles  saturation  bv 
leaving  the  control  at  the  boundary  and  stabilizes  the  vehicle 
using  the  remaining  unsaturated  controls.  Figure  31  shows  the 
state  and  control  responses  for  the  longitudinal  commands  with 
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Figure  30  Comparison  of  Directional  Response  With  and 
Without  Yaw  Nozzle  (A6  .)  at  87.3  m/s  (170 


Figure  31  Longitudinal  Command  Response  at  87.5  m/s  (170  kt) 
With  Nozzle  Angle  (A6n  .  )  Unlimited 


unlimited  A8^  movement  which  are  adequate  except  that  A6^ 
moves  beyond  its  limits  when  a  positive  Act  is  commanded.  The 
state  feedback  gain  for  the  CGT-P1  which  corresponds  to  the 
closed-loop  eigenvalues  in  the  middle  column  of  Table  10  is  as 
follows 


~0.434 

0.586 

0.075 

0.0032 

-0 

-0 

~0 

~0 

-0.611 

0.124 

-0.439 

-0.106 

~0 

~0 

~0 

~0 

0.006 

0.061 

-0.235 

-0.211 

~0 

~0 

~0 

~0 

-0 

-0 

~0 

~0 

0.454 

-0.317 

0.343 

0.961 

-0 

-0 

~0 

-0 

0.209 

-0.725 

0.0732 

-0.024 

The  ordering  of  the  states  and  controls  is  the  same  as  in 
Table  8.  The  control,  AS ^ ,  plays  an  important  roll  in  sta¬ 
bilizing  the  vehicle  with  sizable  pitch  rate  feedback,  -0.235, 
and  the  largest  Aw  feedback  gain,  -0.211.  The  importance  of 
A0 j  in  stability,  but  not  in  command  response,  can  be  down 
played  by  increasing  the  control’s  rate  weighting.  The  state 

feedback  gain,  caused  by  decreasing  A6 .  ,  which  corresponds 

J  *  max 

to  the  closed-loop  eigenvalues  in  the  last  column  of  Table  10 
is 


0.418 

0.584 

0 . 06^4 

0.00125 

-0 

-0.0333 

-0 

-0 

-0.819 

0.167 

-0.586 

-0.190 

-0 

-0 

-0 

-0 

0.209 

0.00446 

-0.0227 

-0.0690 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

-0 

0.253 

-0.0512 

0.277 

0.860 

0.0392 

0.0563 

~0 

-0 

0.0761 

-0.108 

0.0506 

0.108 

-0.0442 

-0.0596 

-0 

-0 

0.418 

-0.783 

0.161 

0.266 

(6) 

Without  substantially  changing  the  closed-loop  longitudinal 
eigenvalues  the  Aq-A8j  and  Aw-A6j  gains  have  been  reduced  by 
almost  an  order  of  magnitude.  When  the  longitudinal  ACM  com¬ 
mands  are  repeated  in  Fig.  32  with  A8j  position  limited  using 
the  feedback  gain  in  Eq .  6,  the  aircraft  remains  stable,  AV 
and  Ay  commands  are  accommodated,  but  the  aircraft  cannot 
significantly  maintain  an  increased  command  in  angle  of  attack 
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while  holding  AV  and  A0j  constant.  The  importance  of  vectored 

thrust  is  to  reduce  angle  of  attack  variations  when  large 

normal  acceleration  commands  are  requested  and  CGT-PI  performs 

this  feature  very  well  for  a  Aa  /V  command  of  1  sec  ^  which 

produces  the  A-y  ramp  in  Fig.  32a  and  the  control  responses  in 

Fig.  32b.  The  maximum  angle  of  attack  variation  before  Aa 

returns  to  its  commanded  value  of  0.0  is  0.1  deg  when  Aa 

6  n  ,m 

is  commanded.  In  flight,  the  angle  of  attack  ACM  command 
could  be  preprogrammed  depending  on  forward  velocity  to  avoid 
0j  commands  beyond  saturation  (with  pilot  override  available 
but  usually  not  used). 

3.3.3  Transition 


This  section  presents  design  results  for  the  two 
transition  flight  conditions,  51.5  m/s  (100  kt )  and  18.0  m/s 
(35  kt ) .  Transition  is  characterized  as  a  shift  from  the 
aerodynamic  control  surfaces  to  the  reaction  jets  for  sta¬ 
bility  and  control.  Two  features  are  expected  to  occur  during 
transition;  the  feedforward  gain  elements  and  the  feedback 
gains  for  A6  ,  A6  ,  and  A6  should  decrease  with  forward  ve- 

b  d  i 

loci ty . 


The  feedforward  matrix  S ^  is  shown  in  the  following 
equations  using  the  enroute  velocity  command  model  for  51.5 
m/s  (100  kt )  and  18.0  m/s  (35  kt),  respectively. 


ANj 

-0.338 

-0 . 838 

0.243 

0.0141 

-0.00506 

-0.0662 

-0.A02 

-0.368 

0.0136 

-0.00665 

AO" 

J 

-0.171 

-1.53 

-1 .8A 

0  0134 

-0.00766 

"h 

AA  . 

P  1 

r  S21 A-m  + 

-0.128 

~0 

-0.529 

-0.0222 

-0.275 

-0.0392 

-0.0127 

0.102 

0.00639 

-0.0478 

-0 

-0  00838 

-0.0133 

-0.256 

0.138 

A\ 

-0.00701 

-0.0171 

0.00626 

0.0470 

-0.00891 

~0 

0.00126 

-0.00107 

0.0713 

-0.0283 

Aa 


DR  ,m 


AV 


t.  ,m 


AP. 


a,\'R.m 

A<i' 
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-0.121 

-0.135 

-0.0196 

~0 

-0.00263 

~0 

-0.00156 

-0.00356 

-0 

-0 

-0.0252 

-1.02 

-1.77 

-0 

-0 

-0.0600 

-0.122 

-0.0969 

-0.00102 

0.00222 

-0 

-0.00113 

-0.00231 

0.00195 

~0 

~0 

-0 

~0 

-0.0161 

0.0209 

~0 

-0 

~0 

0.00277 

-0 

-0 

-0 

~0 

-0.00198 

-0.0131 

The  elements  for  AS  ,  AS  ,  and  AS  decrease  by  an  order  of 

s  a  r 

magnitude  as  expected.  Similarly,  the  feedback  gains  decrease 
as  shown  in  Fig.  33  using  the  A<!>  to  AS  gain  and  the  Aq  to  A6 
gain  as  typical  examples.  Although  only  four  points  are 
plotted  for  each  gain  shown  in  Fig.  33,  indications  are  that 
the  AV-8A  gains  have  easily  scheduled  smooth  variations  with 
flight  condition.  The  CGT-P1  designs  should  readily  convert 
to  gain- scheduled  nonlinear  control  laws  using  the  same  tech¬ 
nique  as  in  Ref.  26. 

The  closed-loop  mapped  eigenvalues  for  CGT-PI  in 
transition  flight  are  shown  in  Tables  11  and  12.  The  open-loop 
characteristics  have  unstable  complex  pairs  in  both  the  longi¬ 
tudinal  and  lateral  modes.  The  feedback  gains  stabilize  all 
modes,  cause  the  complex  pairs  to  have  damping  ratios  greater 
than  0.5,  and  rearranges  the  eigenvectors.  It  should  be  kept 
in  mind  that  the  identification  of  states  with  eigenvalues  is 
an  eigenvector  simplification  useful  in  choosing  Q  and  R.  For 
example.  A®  has  contributions  from  almost  all  the  lateral- 
directional  modes;  identifying  Ad>  with  Ar  in  a  complex  pair  is 
obtained  by  a  process  of  elimination  rather  than  a  recognized 
predominate  involvement.  As  the  forward  velocity  decreases 
(Table  11  to  Table  12),  the  decreased  aerodynamic  effects 
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Figure  33  Feedback  Gain  Variations  With  Forward  Velocity 

TABLE  11 

EFFECTS  OF  CGT-P1  IN  TRANSITION  FLIGHT 


(51 

.  5  m/s 

(100 

kt ) 

i  •  290 8 

OPEN- LOOP 

CHARACTERISTICS 

ALL 

CLOSED- LOOP 

CONTROLS 

CHARACTERISTICS 

NO  REACTION  JETS 

CLOSED- LOOP  CHARACTERISTICS 

DYNAMIC 

MODE 

NATURAL 
FREQUENCY , 
rad/sec 

DAMPING 

RATIO 

TIME 

CONSTANT , 
sec 

NATURAL 
FREQUENCY . 
rad/*ee 

DAMPING 

RATIO 

TIME 

CONSTANT . 
aec 

NATURAL  nAMPinir  ?Ifl£ 

FREQUENCY.  WATir  CONSTANT . 

rad/sec  tec 

TABLE  12 


EFFECTS  OF  CGT-P1  IN  TRANSITION  FLIGHT 


(18.0 

m/s  (35 

kt)  ) 

T-2909 

.  . 

OPEN- LOOP  CHARACTERISTICS 

ALL 

CLOSED- LOOP 

CONTROLS 

CHARACTERISTICS 

DYNAMIC 

MODE 

NATURAL 

FREQUENCY 

rad/sec 

DAMPING 

RATIO 

TIME 

CONSTANT . 

sec 

NATURAL  n.up i 

FREQUENCY  . 
rad/ sec  Ratio 

TIME 

CONSTANT , 
sec 

/A*-/Vcs 

- 

- 

- 

?.69 

0.666 

- 

Aq-Av 

- 

- 

- 

1.91 

0.713 

- 

Ar-As 

- 

- 

- 

1.73 

0.633 

- 

**-'4VDR 

- 

* 

- 

0.706 

0.668 

- 

Au-jAr 

- 

- 

- 

0.437 

0.761 

- 

ROLL  MODE 

- 

- 

1.29 

0.330 

Aw 

- 

- 

- 

0.571 

A8 

- 

- 

- 

0.609 

/A0 

• 

- 

- 

1.33 

AV 

- 

- 

- 

7  .  «*A 

SPIRAL 

- 

- 

15.67 

- 

A&-AU 

0.351 

-0.589 

- 

- 

Aq 

• 

• 

1.47 

- 

Aw 

- 

- 

9.32 

- 

Av-Ar 

0.672 

-0.310 

* 

* 

■ 

cause  the  closed-loop  poles  to  move  closer  to  the  imaginary 
axis  but  with  approximately  the  same  damping  ratio  in  the 
complex  pairs.  The  pseudoinverse  CGT-PI  design  has  the  same 
number  of  closed-loop  eigenvalues  with  and  without  the  reac¬ 
tion  jets  in  transition  flight.  Removing  the  reaction  jets  as 
shown  n  Table  11  causes  a  marked  degradation  in  system  closed 
loop  performance  for  fixed  Q  and  R. 

The  most  interesting  aspect  of  the  transition  designs 
is  in  comparison  between  the  dynamic  and  ramp/step  models  for 
coordinated  turns.  Figures  34  and  35  show  state  and  control 
transients  for  step  pitch  angle  commands  from  each  type  of 
command  generator,  while  Figs.  36  and  37  show  state  and  con¬ 
trol  transients  for  ramp  yaw  angle  commands.  The  dynamic 
models  used  for  model  following  have  eigenvalues  near  the  same 
regions  in  the  complex  plane  as  the  closed-loop  eigenvalues 
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for  CGT-PI.  If  this  is  a  reasonable  design  practice  (i.e.,  a 
poorly  damped,  slow  dynamic  model  is  not  a  desirable  command 
generator)  then  dynamic  models  offer  little  advantage  over 
ramp/step  command  generators  and,  in  fact,  have  a  clear  dis¬ 
advantage  from  an  implementation  point  of  view  (see  Section 
3.2.2).  The  dynamic  model  causes  the  aircraft  to  have  de- 

JL 

sirable  transients  as  the  sta  e  Ax  tracks  Ax"  as  expected,  but 
the  Ax  transient  for  the  ramp  and  step  models  are  just  as 
good.  Stable  dynamic  models  for  command  generator  tracking 
have  been  researched  starting  from  Ref.  17  but  it  is  onl}7 
recently,  Refs.  11  and  12,  that  theory  has  become  available 
for  designing  ramp  and  unstable  command  generator  trackers 
using  optimal  control . 

The  last  two  simulation  results  for  transition  flight 

show  the  five  command  responses  for  18.0  m/s  (35  kt)  in  Fig. 

38  using  the  enroute  velocity  command  model  and  the  effect  of 

putting  a  square  pulse  in  the  vertical  velocity  command  at 

51.5  m/s  (100  kt)  in  Fig.  39.  The  pitch  angle  and  yaw  angle 

response  at  18.0  m/s  (35  kt)  are  almost  the  same  as  the  51.5 

m/s  (100  kt)  design  responses  in  Figs.  34  and  35  indicating 

desirable  closed-loop  uniformity  during  transition.  The 

square  pulse  in  Fig.  39  demonstrates  an  important  property  of 

the  CGT-PI.  Even  though  the  design  requires  the  command  model 

input  to  be  constant;  the  input  in  implementation  can  be  very 

time-varying  without  causing  poor  transients.  If  inputs  are 

slowly  varying  or  are  changed  after  the  system  error  is  small, 

the  resulting  plant  control  inputs  are  essentially  optimal. 

In  Fig.  39,  the  pilot  control  command  of  0.305  m/s  (1  fps )  for 

AV  begins  the  descent,  then  just  as  Ax  approaches  Ax  ,  the 

AV  is  returned  to  zero.  The  AV-8A  control  transients  to 
z  ,m 

begin  and  stop  the  descent  are  essentially  optimal  mirror 
images  of  each  other. 
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The  AV-8A  in  low  altitude,  low  speed  flight,  as  with 
all  V/STOL,  is  very  difficult  for  a  pilot  to  control  and  is 
known  to  have  excessive  pilot  workload  (Ref.  3).  All  of  the 
design  tools  developed  in  Appendices  A  and  B  are  applied  to¬ 
ward  relieving  pilot  workload  in  low  speed  flight.  This  sec¬ 
tion  will  compare  the  attitude  and  velocity  command  systems, 
study  the  ground  effect  (Appendix  C),  and  present  simulations 
from  optimal  control  and  eigenvalue/eigenvector  placement 
designs . 


Table  13  shows  the  eigenvalue  characteristics  for  the 
AV - 8 A  with  velocity  or  attitude  commands.  The  only  signifi¬ 
cant  difference  is  that  the  Av-Ar  (or  Dutch  roll  like)  mode  in 
the  velocity  system  splits  in  the  attitude  command  system  to 
form  an  almost  neutrally  stable  lateral  speed  mode  associated 
with  Av.  The  mode  associated  with  Ar  combines  with  an  inte¬ 
grator  compensator  state  to  form  a  complex  pair.  Any  attempt 
to  stablize  the  Av  mode  in  the  attitude  command  system  de¬ 
grades  the  steady  state  behavior  of  the  A<t>m  command.  The  dif¬ 
ference  between  velocity  and  attitude  closed-loop  eigenvalues 
is  typical  behavior  and  is  also  reported  in  Ref.  26  for  a 
helicopter  study.  Figure  40  compares  state  responses  for 
velocity  and  attitude  commands.  When  yaw  angle  is  commanded 
to  be  a  ramp  with  the  three  positions  command  to  zero  in  the 
velocity  command  system,  the  AV-8A  turns  over  a  spot  on  the 
ground  with  a  small  ramp  increase  in  roll  angle  as  shown  in 
Fig.  40c.  When  a  cross-range  velocity  command  of  0.305  m/s 
(1  fps )  with  yaw  angle  commanded  to  zero  is  requested  for  the 
velocity  command  system  the  most  noticeable  effect  in  Fig.  40d 
is  the  roll  angle  peak  transient  of  approximately  1.3  degrees. 
A  (not  unreasonable)  step  crosswind  gust  disturbance  of  5.15 
m/s  (10  kt)  could  cause  the  velocity  command  system  to  bank 
the  aircraft  near  20  deg  in  transient  to  maintain  position. 
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The  bank  angle  transient  in  the  velocity  command  system  can  be 
reduced  at  the  expense  of  increased  lateral  drift. 

An  alternative  is  to  give  the  pilot  full  authority  in 
correcting  lateral  drift  using  the  attitude  hold  command  sys¬ 
tem  with  an  explicit  roll  command.  Any  crosswind  causes  the 
aircraft  to  drift  slowly  with  constant  attitude.  V/STOL 
pilots,  reported  in  Refs.  7  and  17,  prefer  a  level  attitude  in 
the  descent  to  touchdown  and  do  not  actively  correct  for  small 
lateral  drift. 

As  the  V/STOL  descends,  the  vehicle  enters  the  ground 
effect  with  subsequent  changes  in  characteristics,  particularly 
in  roll.  Table  14  shows  the  closed-loop  eigenvalues  for  the 
optimal  in  ground  effect  (IGE)  gains  and  indicates  what  happens 
if  the  out  of  ground  effect  (OGE)  gains  used  in  Table  13  are 
used  with  the  IGE  AV-SA  model.  The  lateral-directional  mode 
Ap-jAb  becomes  very  lightly  damped.  A  check  of  the  feedback 
gains  in  Fig.  33  shows  that  the  A<t>  to  A6rj  gain  should  be  con¬ 
siderably  reduced  when  IGE  occurs.  Both  of  these  observations 
are  in  agreement  with  Ref.  9  where  it  is  demonstrated  than 
pilots  should  not  have  tight  attitude  control  IGE  or  P10 
(pilot  induced  oscillation)  occurs  in  roll.  The  CGT-PI  con¬ 
trol  law  would  require  special  gain  adaptation  procedures  to 
avoid  poor  performance  as  the  vehicle  descends  to  touchdown . 

The  pilot  in  low  speed  flight  for  the  AV-8A  wi thout 
command  augmentation,  has  three  longitudinal  controls;  pitch 
attitude,  thrust  magnitude,  and  thrust  direction  to  control 
vehicle  velocity  and  orientation  in  the  longitudinal  axis.  The 
controls  require  rapid,  simultaneous,  coordinated  activity  if 
the  pilot  wishes  to  change  altitude  without  changing  pitch 
angle  and  speed,  and  likewise  for  changing  speed  without 
changing  altitude  and  pitch  angle.  It  is  recognized  (Ref.  3) 
that  a  suitably  decoupled  command  augmentation  fly-by-wire 


TABLE  14 


EFFECTS  OF  GAIN  ADAPTATION  IN  THE  GROUND  EFFECT  AT 
5.15  m/s  (10  kt)  FOR  THE  VELOCITY  COMMAND  SYSTEM 


l-.'Oll 


CLOSED- LOOP  CHARACTER I ST! CS 
ADAPTED  IGE  GAINS 

CLOSED-LOO]'  CHARACTER  1ST  ICS 
UNADAPTED  OGE  GAINS 

DYNAMIC 

mode: 

NATURAL 

FREQUENCY 

rad/sec 

DAMPING 
’  RATIO 

TIME 

CONSTANT , 
sec 

NATURAL 

FREQUENCY 

rad/sec 

DAMPING 
’  RATIO 

TIME 

CONSTANT. 

sec 

ad-  /  Ai|' 

2.05 

0.89 

- 

1.89 

0.0396 

- 

Aq-Aw 

1.69 

0 . 563 

- 

1.37 

0.621 

- 

Av-Ar 

1.52 

0.577 

- 

1  .54 

0.571 

- 

AX- fAX 

0.582 

0.612 

- 

0.588 

0.621 

- 

Az-  fAZ 

0.671 

0.728 

- 

0.610 

0 . 604 

- 

Ay- f  Av 

0.593 

0.682 

- 

0.482 

0.731 

- 

A0 

- 

- 

0.561 

- 

- 

0.371 

Ait 

- 

- 

0.354 

- 

- 

0 . 200 

At) i 

- 

- 

0.678 

- 

- 

0.685 

/A0 

- 

- 

1  .25 

- 

- 

1.05 

Au 

- 

- 

2.02 

- 

- 

2.02 

system  is  required  for  VSTOL  aircraft  to  meet  control  objec¬ 
tives  and  that  display  improvements  alone  will  not  overcome 
deficiencies.  The  direct  digital  designed  eigenvalue/ 
eigenvector  placement  CGT-PI  is  an  ideal  candidate  for  meet¬ 
ing  these  control  requirements. 

The  eigenvalue/eigenvector  placement  (EEP)  model 
presented  in  Appendix  B.3  decouples  forward  velocity  and 
altitude  response  while  placing  the  closed-loop  eigenvalues  at 
approximately  the  same  places  as  the  optimal  control  designed 
velocity  command  system.  The  PI  gains  for  the  EEP  design  are 
(the  ordering  of  the  longitudinal  states  and  controls  are 
shown  in  Table  B-l  ) 
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0.0005 

-2.16 

2.93 

0.298 

3.09 

-2.19 

-2.  IS 

0.112 

0.737 

-3.31 

-0.186 

-0.308 

-0.106 

-0 . 581 

1.73 

-0.0696 

1 .  39 

-0.59A 

AtC0 


0.00001 

-0.0A56 

-0.00221 


-0.039A  -0.0656 

-0.00355  0.003A1 

0.0312  -0.0176 


while  the  optimal  control  gains  are 


0.0911 

-0.910 

0.0571 

0.292 

-0.101 

-1  ,2A 

-0.979 

0.0325 

0.17A 

-1.71 

0.C198 

-0.181 

-0.A37 

-0.256 

1.30 

-0.539 

0.681 

-0.396 

A  tC2 


0.002A  -0.0210  -0.0266 

-0.02A9  -0.017A  0.000891 

-0.0115  0.0702  -0.000127 


Optimal  control,  as  expected,  has  smaller  gains  (i.e.,  lower 
control  variances  when  noise  is  present)  and  considerably 
lower  gains  for  the  thrust  magnitude  control;  all  desirable 
features.  A  comparison  of  the  response  of  the  tw7o  designs  is 
shown  in  Figs.  A1  to  A3  for  the  velocity  command  system.  The 
two  low  frequency  complex  pairs  associated  with  position  in 
the  optimal  control  design  are  replaced  by  stable  real  eigen¬ 
values  in  the  EEP  system  and  this  reflects  into  less  oscilla¬ 
tory  control  movements  particularly  when  down-range  velocity 
is  commanded  in  Fig.  A3.  The  comparison  of  the  two  designs 
indicates  that  the  direct  digital  EEP  design  has  the  potential 
for  competing  with  and  complementing  optimal  control  and  is  an 
easier  method  for  evolving  time  domain  design  objectives  into 
actual  control  performance.  The  EEP  design  simulated  in  the 
figures  is  a  first  time,  one  iteration  effort  that  did  not 
take  advantage  of  the  considerable  design  freedom  available 
in  the  choice  of  the  ideal  model. 
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EIGENVALUE/EIGENVECTOR  PLACEMENT 


Comparison  of  Optimal  Control  and 
Eigenvalue/Eigenvector  Placement  Pitch 
Angle  Step  Responses  at  5.15  m/s  (10  k 


genvalue/Eigenvector  Pla 
sit  ion  Ramp  Responses  at 


3. A 


CHAPTER  SUMMARY 


This  chapter  uses  the  recently  developed  theory  of 
discrete- time  feedforward  control  coupled  with  two  powerful 
direct  digital  multivariable  stabilization  design  procedures 
to  obtain  command  augmentation  systems  for  the  AV-8A  Harrier. 

The  feedforward  control  allows  the  design  of  output  model  fol¬ 
lowing  for  all  types  of  command  generator  models;  these  models 
are  constructed  such  that  the  aircraft  follows  their  output 
trajectory.  Six  different  command  generator  models,  especially 
formulated  for  V/STOL  control,  are  constructed,  analyzed,  and 
simulated.  The  two  stabilization  design  procedures  can  interact 
equally  with  each  of  the  command  generator  models.  Each  stabi¬ 
lization  procedure  has  the  desirable  Type  1  property,  and  each 
procedure  can  be  implemented  using  current  fly-by-wire  all 
digital  technology. 

The  results  with  the  command  generator  models  de¬ 
monstrate  that  rate  commands,  with  attitude  or  velocity  hold, 
(the  ramp/step  designs)  commonly  used  in  current  classical 
analog  designs  can  also  be  accomplished  with  optimal  control. 

In  a  comparison  between  an  optimal  control  ramp/step  design 
and  an  optimal  control  dynamic  model  design,  (the  dynamic 
model  being  the  common  optimal  control  multivariable  approach 
in  the  literature),  the  ramp/step  approach  yields  comparable 
performance  with  less  implementation  complexity.  Five  of  the 
six  command  generators  investigated  are  ramp/step  designs. 

Each  of  the  command  generators  successfully  uses  five  pilot 
inputs  to  command  five  aircraft  outputs  using  all  available 
controls  independently.  Thrust  vectoring  and  reaction  control 
at  high  forward  velocity  is  shown  to  be  difficult,  but  with 
proper  considerations  is  feasible  to  design  and  certainly  im¬ 
proves  AV-8A  performance  in  high  speed  flight.  The  transition 
and  terminal  area  commands  are  velocity  oriented,  which  is  in 
agreement  with  pilot  opinions.  The  results  in  this  chapter 
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indicate,  however,  that  in  low  speed  flight,  a  roll  command 
rather  than  a  cross-range  velocity  command  is  more  appropriate 
because  pilots  are  reported  to  prefer  a  level  platform  rather 
than  automated  lateral  drift  correction  as  the  vehicle  de¬ 
scends  to  touchdown. 


The  stability  augmentation  designs  for  the  AV-8A 
exhibit  uniform  handling  qualities  over  the  flight  conditions 
investigated .  Figure  44  shows  two  of  the  closed-loop  mapped 
high  frequency  eigenvalues  overlaid  on  a  specifications  plot, 
taken  from  Ref.  32.  Results  from  Ref.  36  indicate  desirable 


positions  of  V/STOL  complex  pair  eigenvalue  locations  for 
Dutch  roll  (lat  in  Fig.  44)  and  short  period  (Ion  in  Fig.  44). 
The  three  V/STOL  flight  condition  results  from  the  design  in 


this  chapter  all  fall  within  the  optimum  frequency  and  dampin 


window.  When  the  V/STOL  enters  high  speed  flight,  the  aero 
dynamic  controls  take  over  at  low  dynamic  pressure  dropping 


UNDAMPED  NATURAL  FREQUENCY  w„  (RAD/SECI 


Figure  44  CGT-PI  High  Frequency  Complex 
Pair  Performance 
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the  natural  frequency  of  short  period  and  the  damping  of  Dutch 
roll  as  shown  in  Fig.  44. 


A  new  direct  digital  design  technique  is  presented 
and  successfully  applied  to  the  AV-8A  at  low  speeds.  The 
technique  is  shown  to  be  comparable  to  optimal  control  results 
and  uniquely  specifies  a  stability  augmentation  system  by 
eigenvalue/eigenvector  placement  of  the  AV-8A  dynamics.  The 
new  technique  readily  transforms  designer  time  domain  objec¬ 
tives  such  as  decoupled  response  for  the  AV-8A  at  low  speeds 
and  pole  placement  into  practice,  but  frequency  domain  perfor¬ 
mance  would  have  to  be  checked  for  acceptability. 

In  summary,  feedback  and  feedforward  control  theory 
is  available  for  designing  algorithms  which  greatly  relieve 
pilot  workload  in  commanding  and  stabilizing  the  complex, 
unstable,  nonlinear  dynamics  of  the  V/STOL  aircraft.  The 
discrete- time  feedforward/linear  optimal  control  design  com¬ 
bination  can  be  a  desirable  candidate  for  improving  and  ex¬ 
panding  the  stability  and  control  characteristics  of  current 
and  future  aircraft. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  has  presented  a  number  of  advanced  tech¬ 
niques  for  designing  the  stability  and  control  systems  of  high 
performance  aircraft  and  related  performance  results.  In  par¬ 
ticular,  the  digital  synthesis  procedure  was  performed  on  the 
AV-8A  Harrier  in  all  flight  regimes;  ACM,  high  speed  flight, 
transition  and  hover.  The  topics  which  have  been  addressed 
include  the  stability  characteristics  of  the  Harrier,  new 
command  and  stabilization  control  theory,  and  it's  application 
to  the  Harrier. 


h . 1  RESULTS 


The  results  obtained  in  this  investigation  fall  into 
eight  categories: 


•  Characterization  of  hovering,  low-speed 
and  transition  flight  stability  relation¬ 
ships,  including  the  effects  of  velocity, 
angle  of  attack,  flight  path  angle,  and 
acceleration 

•  Characterization  of  stability  during 
a  ski-jump  launch 

•  Characterization  of  low-speed  stability 
in  "ground  effect"  conditions 

•  Quantitative  description  of  vectored 
thrust  and  reaction  jet  effects  on  the 
ability  to  stabilize  and  control  the 
aircraft  in  maneuvering  flight 

•  Development  of  an  eigenvalue/eigenvector 
placement  control  law  for  an  alternative 
synthesis  technique  to  optimal  control 


theory  and  to  obtain  multivariable  de¬ 
coupling 

•  Development  of  a  feedforward  command 
theory  (the  command  generator  tracker) 
to  ensure  asymptotic  output  tracking  of 
commands.  A  "Type  1"  structure  is  added 
to  account  for  aircraft  parameter  varia¬ 
tions 

•  Design  of  six  pilot  command  models  for 
precise  handling  and  tracking  of  appro¬ 
priate  aircraft  variables  during  all 
flight  regimes 

•  Preliminary  design  and  analysis  of  a 
digital  command  augmentation  system 
for  the  subject  aircraft. 

Actual  results  under  each  of  these  categories  are  discussed  in 
the  following  paragraphs. 

Stability  of  the  aircraft  is  examined  using  linear, 
time-invariant  dynamic  models.  Stability  derivatives  are 
based  on  the  AV-8A  aerodynamic  data  contained  in  Ref.  6. 
Vehicle  stability  evaluation  proceeds  in  a  manner  similar  to 
the  fighter  aircraft  analyses  reported  in  Refs.  4  and  5.  The 
linear  dynamic  models  include  all  the  effects  of  inertial  and 
aerodynamic  cross-couplings,  and  also  model  the  effects  of 
vehicle  acceleration  on  stability. 

Although  the  general  V/STOL  model  used  here  includes 
the  independent  effects  of  each  aerodynamic  and  reaction  con¬ 
trol,  the  trim  conditions  for  stability  analysis  are  found 
assuming  the  same  SAS-off  control  interconnects  as  on  the 
actual  Harrier.  Trim  conditions  are  found  in  accelerated  or 
unaccelerated  flight  using  the  nonlinear  aircraft  dynamics. 

For  the  work  which  examines  vectored  thrust  and  control  de¬ 
sign,  the  reaction  jets  and  aerodynamic  surface  controls  are 
treated  separately  to  maximize  the  payoff  of  the  extra  control 
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freedom.  Stability  relationships  that  are  observed  include 
those  discussed  in  the  following  paragraphs: 

Longitudinal  stability  is  essentially  a  function  of 
nozzle  angle,  0^.  Above  50  deg  0  ^  ,  the  static  instability  due 
to  ram  drag  expresses  itself  as  a  relatively  slow,  unstable 
oscillation  called  the  third  mode,  so-called  because  it  is 
inherently  different  from  the  conventional  short  period  and 
phugoid  modes.  For  0^  below  45  deg,  the  static  instability 
appears  as  an  angle-of-attack  divergence  that  is  even  more 
unstable  at  lower  0 j .  Finally,  for  the  conventional  flight 
configuration  (0^  =  0  deg)  at  higher  flight  speeds,  the  in¬ 
stability  completely  disappears  and  conventional,  stable, 
phugoid  and  short-period  modes  appear.  The  shipboard  ex¬ 
perience  with  the  AV-8A,  Ref.  7,  agrees  with  these  results. 

Lateral  stability  is  primarily  a  function  of  velocity 
and  angle  of  attack,  with  the  difficulty  arising  due  to  Dutch 
roll  instability  in  the  higher-a ,  lower-speed  regime.  In  the 
lateral  vehicle  response,  the  classic  roll  and  spiral  modes 
and  the  Dutch  roll  oscillation  describe  the  vehicle  response 
throughout  the  speed  regime . 

Ground  effect  causes  a  large  and  very  undesirable 
change  in  vehicle  stability.  Fast,  unstable  roll  and  altitude 
divergences  appear  that  essentially  prohibit  any  steady  flight 
in  ground  effect. 

Sideslip  envelope  restrictions  in  the  transition 
speed  range  are  significant  due  to  a  large  positive  dihedral 
effect  and  limited  roll  control  power,  as  has  been  discussed 
in  Ref.  8.  A  significant  increase  in  roll  control  power  on 
the  YAV-8B  lifts  the  sideslip  envelope  restrictions. 
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The  AV-8A  stability  during  a  ski-jump  launch  is 
studied.  This  launch  technique  (Ref.  9)  allows  a  V/STOL  air¬ 
craft  to  take  off  at  higher  gross  weight  or  lower  airspeed 
than  is  possible  with  a  conventional  short  take  off.  The 
stability  study  concludes  that  high  angles  of  attack  should  be 
avoided  until  the  speed  is  high  enough  to  justify  the  Dutch 
roll  instability  that  worsens  as  a  increases.  The  highest 
value  of  acceleration  along  the  flight  path  also  requires  the 
smallest  value  of  nozzle  angle,  and  therefore  produces  the 
fastest  longitudinal  instability.  An  optimum  (least  unstable) 
trade-off  between  normal  and  longitudinal  acceleration  may 
exists . 


Vectored  thrust  and  reaction  control  in  forward 
flight  for  the  AV-8A  are  examined  using  a  series  of  command 
models  and  the  new  theory  of  feedforward  control  reported  in 
Refs.  10  to  14.  Feedforward  control  solves  the  general  output¬ 
following  problem:  given  any  desirable  model  to  be  followed, 
either  stable  or  unstable,  how  should  the  aircraft  controls 
move  and  what  happens  to  the  aircraft  states  when  an  output 
vector  of  the  aircraft  states  perfectly  tracks  an  output 
vector  of  the  model  states?  When  the  dimension  of  the  output 
of  the  aircraft  states  is  equal  to  or  less  than  the  number  of 
controls,  the  feedforward  control  problem  almost  always  has  a 
solution . 


Feedforward  control  solutions  are  performed  for  the 
AV-8A  using  six  different  command  models.  Each  command  model 
is  adapted  for  a  particular  aspect  of  flying  the  V/STOL,  such 
as  for  air  combat  maneuvering  in  high  speed  flight,  or  for 
3-dimensional  position  control  near  hover.  Even  though  the 
generalized  AV-8A  can  have  eight  separately  adjustable  con¬ 
trols  only  five  outputs  of  the  aircraft  can  be  commanded. 
There  is  insufficient  lateral -directional  control  power  to 
independently  command  three  lateral-directional  outputs  for 
the  AV-8A. 
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Vectored  thrust  at  high  velocities  with  proper  con¬ 
trol  phasing  allows  for  large  normal  acceleration  commands 
with  minimal  increases  in  angle  of  attack  as  demonstrated 
using  simulations.  Near  hover,  proper  movement  of  thrust 
magnitude  and  nozzle  angle  along  with  independent  pitch  jet 
control  allows  for  complete  decoupling  of  axial  and  vertical 
motion.  Independent  use  of  reaction  jets  in  the  lateral- 
directional  axes  at  high  velocities  is  shown  to  decrease 
sideslip  excursions  in  rolling. 

A  highlight  of  the  AV-8A  digital  control  design  study 
is  that  a  significant  contribution  to  linear  control  system 
theory  has  been  made  which  has  a  direct  and  powerful  impact  on 
practical  digital  command  augmentation  system  design.  Feed¬ 
forward  control  theory,  linear  optimal  regulator  theory,  and 
proportional-integral  control,  all  in  discrete- time ,  are 
combined  to  form  the  command  generator  tracker-proportional 
integral  (CGT-PI )  full  authority  controller.  Using  CGT-PI , 
states  with  non-constant  steady  state  values  are  weighted  in 
the  inf inite- time  discrete  cost  function  and  yield  optimal 
results.  Inner-loop  stability  and  outer-loop  automatic  guid¬ 
ance  are  simultaneously  designed  and  simulated  using  CGT-PI. 
For  the  first  time,  a  comparison  is  made  between  a  model  to  be 
followed  which  has  idealized  dynamics  but  the  same  dimension 
as  the  aircraft  model  and  models  to  be  followed  which  are 
simply  integrators  with  a  lower  dimension  than  the  aircraft 
model.  The  idealized  high  dimensional  dynamics  model  is  shown 
not  to  offer  appreciable  advantages  over  the  simpler  models. 
For  the  first  time,  a  comparison  is  made  between  discrete- time 
optimal  control  and  discrete- time  eigenvalue/eigenvector  place 
ment.  Both  stabilization  procedures  have  the  desirable  closed 
loop  eigenvalue  mapping  property  reported  in  Ref.  13  which 
enables  the  digital  design  to  be  evaluated  using  continuous¬ 
time  specifications.  The  results  show  that  each  approach  has 
its  own  advantages.  The  best  procedure  is  to  use  the  concepts 
together  for  maximum  benefits. 
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The  preceding  paragraphs  have  described  the  contribu 
tions  made  in  understanding  the  handling  qualities  of  a  jet 
V/STOL  in  hover,  transition,  and  high  speed  flight.  New 
practical  control  results  which  are  better  able  to  handle  the 
demanding  multivariable  control  complexity  of  the  jet  V/STOL 
are  now  available  as  a  result  of  this  study. 


4.2  CONCLUSIONS 


This  investigation  provided  the  following  signigicant 
results  and  conclusions: 


•  Stability  Characteristics 

The  linear,  time-invariant  model,  generated 
by  linearizing  the  vehicle  about  trim  points 
(flight  conditions),  exhibited  several  un¬ 
stable  characteristics.  Longitudinal  sta¬ 
bility  is  greatly  influenced  by  nozzle  angle; 
angles  less  than  45°  should  be  avoided  in  for¬ 
ward  flight,  if  consistent  with  performance 
bounds.  Lateral  stability  is  influenced  by 
velocity  and  angle  of  attack;  low  speed,  high 
angles-of-attack  destabilize  the  Dutch  roll 
mode.  Thus,  sufficient  velocities  should  be 
obtained  in  the  ski-jump  launch  before  angle- 
of-attack  is  increased.  Ground  effect  causes 
large  destabilizing  divergence  in  roll  and 
altitude;  steady  flight  is  essentially  pro¬ 
hibited  . 

•  Departure-Prevention  Stability  Augmentation 

Control  laws  which  prevent  vehicle  instability 
and  departure  from  controlled  flight  can  be  de¬ 
signed  readily.  The  optimal  linear  regulator 
handles  coupled  dynamics  by  crossfeed  paths 
which  guarantee  stability.  An  alternative 
method  was  developed  here  to  accomplish  this 
and  decouple  aircraft  modes  as  well.  The 
desired  closed-loop  eigenvalues  (poles)  and 
eigenvectors  (modes)  are  chosen  in  the  design 
process  and  the  augmented  vehicle  will  emulate 


these  characteristics.  These  techniques 
enabled  the  augmented  Harrier  to  perform 
according  to  MIL  specifications. 

Departure-Prevention  Command  Augmentation  Systems 


Precision  tracking  response  to  pilot  commands  is 
afforded  by  employing  the  command  generator 
tracker  theory  developed  in  this  report.  Feed¬ 
forward  gains  are  generated  using  this  technique 
which  command  the  aircraft  to  track  the  pilot’s 
commands.  A  "Type  1"  proportional  integral 
control  was  included  to  ensure  perfect  tracking 
for  wide  variations  in  the  aircraft's  parameters. 
This  technique  enabled  the  Harrier  to  perform 
properly  in  all  flight  regimes  to  commanded 
inputs . 

Pilot  Command  Models 

Adaptive  command  models  were  employed  to  obtain 
precise  handling  and  tracking  qualities  at  all 
flight  regimes.  Six  different  command  models 
were  constructed  for  a  particular  aspect  of 
V/STOL  flight  and  the  gains  can  be  scheduled 
as  a  function  of  flight  condition.  Hover 
requires  three-dimensional  position  control 
and  combat  maneuvering  requires  accelera¬ 
tion  commands  --  up  to  five  commands  may  be 
given  simultaneously  and  independently. 

This  scheme  was  necessary  for  a  pilot  to 
precisely  handle  the  Harrier  under  such 
varying  flight  regimes. 

CCV  Design 

The  above  techniques  were  employed  in  con¬ 
structing  the  AV-8A  control  laws  since  eight 
separate  controls  are  available  to  obtain 
six  degrees  of  vehicle  freedom.  The  con¬ 
troller  was  generated  for  ACM,  enroute, 
transition,  terminal  velocity  and  attitude 
flight  regimes.  Performance  was  superb 
in  all  regimes  --  except  attitude  control 
at  Hover.  Side  face  thrust  must  be  im¬ 
proved  . 


4.3 


RECOMMENDATIONS 


It  is  recommended  that  the  following  studies  be 
undertaken  to  extend  and  demonstrate  the  utility  of  the  work 
described  in  this  report. 


Evaluate  the  CGT  Using  Nonlinear  Aircraft 


Simulation 

Using  the  digitally  implemented  controller 
and  gain  adaptation  logic ,  evaluate  the 
Type  1  CGT  and  pilot  command  models  by 
numerical  simulation. 

Pilot-Aircraft  Stability 

Determine  pilot-aircraft  stability  boundaries 
as  a  function  of  the  pilots'  adaptation  to 
changing  flight  conditions.  Construct  pilot 
strategies  for  controlling  the  "bare"  air¬ 
frame  and  suggest  ARI  configurations  to 
alleviate  difficulties. 

Evaluate  the  Sensitivity  of  Controller  Gain 


Schedules  to  the  Aircraft  Model 

Evaluate  the  robustness  of  the  control 
law  with  respect  to  aircraft  parameter 
or  trajectory  variations. 

Effects  of  Partial  State  Feedback 

The  aircraft  sensors  employed  to  measure 
the  aircraft  degrees  of  freedom  and  landing 
coordinates  require  state  estimators  to 
filter  noise  properly  and  reconstruct  any 
non-measured  states.  The  inclusion  of  this 
effort  must  be  considered  before  implemen¬ 
tation  in  the  actual  aircraft.  A  final 
design  can  be  issued  at  this  time. 

Investigate  the  Computational  Advantages  of 


Multi-Rate  Control  and  Estimation 


The  size  of  the  onboard  computer  is  largely  a 
function  of  algorithm  complexity  and  update 
rate.  The  advantages  of  sampling  channels  at 
various  rates  to  reduce  computational  require 
ments  should  be  evaluated. 


APPENDIX  A 

DISCRETE-TIME  COMMAND  GENERATOR  TRACKER 


A  roost  basic  control  problem  is  to  require  a  plant, 
such  as  the  V/STOL,  to  respond  to  command  inputs,  in  the 
presence  of  disturbances,  such  that  the  output  of  the  plant 
maintains  a  desired  trajectory  in  real  time.  A  method  for 
solving  this  tracking  problem  is  to  formulate  the  desired  tra 
jectory  and  disturbances  as  the  output  of  a  system  of  linear 
differential  equations.  The  tracking  problem  then  reduces  to 
finding  the  appropriate  feedforward  gains  from  the  command 
generator  linear  system  states  and  disturbance  states  to  the 
plant  control  inputs  so  that  command  generator  tracking  is 
obtained.  It  is  the  purpose  of  this  appendix  to  solve  the 
command  generator  tracking  (CGT)  problem  for  linear,  time- 
invariant,  discrete- time  plants. 

The  solutions  of  the  command  generator  tracking 
problem  for  discrete- time ,  linear,  time- invariant  plants  and 
models  is  a  straightforward  analogy  of  the  continuous- time 
solution  presented  in  Ref.  26.  The  only  new  concept  needed, 
in  both  the  continuous  and  discrete- time  cases,  is  the  intro¬ 
duction  of  the  star  trajectory.  The  star  trajectories  are 
curves  in  time  that  the  plant  states  and  controls  must  follow 
so  that  the  plant  output  perfectly  matches  the  model  output 
when  modeled  disturbances  affect  the  plant.  The  star  tra¬ 
jectory  primarily  serves  as  a  notational  convenience  and  an 
aid  in  constructing  the  command  generator  tracking  control 
laws . 


The  star  trajectory  is  a  generalization  of  steady- 
state  and  has  appeared  for  continuous- time  plants  in  Refs.  37 
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and  38.  The  star  trajectory  is  determined  by  solving  a  set  of 
'  algebraic  matrix  equations.  Two  new  numerical  procedures  are 

shown  in  this  appendix  to  solve  the  algebraic  equations.  An 
important  by-product  of  the  numerical  procedures  is  the  indi¬ 
cation  of  sufficient  conditions  for  a  solution  of  the  alge¬ 
braic  equations.  The  transmission  zeroes  of  the  discrete 
^  plant  play  an  important  role  in  understanding  the  sufficient 

conditions  for  finding  the  star  trajectories.  Transmission 
zeroes  are  the  concept  of  zeroes  of  a  transfer  function  gen¬ 
eralized  to  multivariable  systems  and  have  become  an  important 
concept  in  the  control  literature  (Ref.  39).  The  first  part 
of  this  appendix  defines  discrete- time  transmission  zeroes  and 
solves  the  CGT  problem  for  linear,  time-invariant  plants. 

1 

1  An  important  feature  of  the  command  generator  tracker 

is  that  it  can  be  used  in  combination  with  classical  or  modern 
feedback  control  design  procedures.  Proper  feedback  control, 
especially  proportional  plus  integral  (PI)  control,  causes  any 
errors  between  the  actual  plant  trajectories  and  the  ideal 
plant  trajectories  to  decay  to  zero.  PI  control  laws  are 
useful  when  used  with  the  CGT  because  the  integral  compensator 
estimates  the  model  control's  feedforward  contribution  to  the 
plant  control,  forces  the  difference  between  plant  output  and 
model  output  to  zero  for  certain  types  of  plant  modeling 
errors  and  can  accommodate  unmodeled  bias  disturbances  in  the 
plant.  The  discrete- time  CGT-PI  combination  is  the  structure 
used  with  optimal  control  theory  to  design  the  control  laws 
for  the  V/STOL  in  Chapter  4.  A  derivation  of  the  CGT-PI 
control  law  is  shown  in  the  second  part  of  this  appendix. 

This  appendix  concludes  with  a  demonstration  of 
eigenvalue/eigenvector  assignment  for  discrete- time  plants. 
Eigenvalue/eigenvector  assignment  of  the  plant  modes  can  be 
accomplished  using  a  special  version  of  the  CGT  and  is  a  form 
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of  implicit  model  following.  An  example  using  CGT  eigenvalue/ 
eigenvector  placement  is  applied  to  the  V/STOL  in  Chapter  A 
and  shows  how  feedback  can  be  used  to  decouple  the  modes  of 
the  plant. 

A . 1  THE  STAR  TRAJECTORY 

A  discrete-time ,  linear,  time-invariant  plant  may  be 
represented  as 


“k+l  =  *“k  ♦  riIk  - 


Ay^  =  HAx^  +  DAu^  +  BAw^ 


AHk+i  *  ♦„4«k 

where  Ax^  is  an  n-state  vector,  Au^  is  an  m-control  vector, 

Ay^  is  an  £-observation  vector,  Aw^  is  a  p-disturbance  vector 

and  <t> ,  T,  E,  H,  D,  B,  4>  are  constant  matrices  of  suitable 

w 

dimension.  It  is  assumed  that  the  pair  (4>,r)  is  stabilizable 
and  the  pair  (4>,H)  is  detectable.  A  command  generator  (or 
model  to  be  followed)  can  similarly  be  expressed  as 


Ax  ,  ,  =  4>  Ax  ,  +  T  Au  , 

— m,k+l  m  — m , k  m  — m , k 


Ay  ,  =  H  Ax  ,  +  D  Au  . 
*-m ,  k  m  - ro ,  k  m  — m ,  k 


Au  .  -  Au  .  ,  =  0 
— m ,  k  — m ,  k  - 1 


where  Ax  ,  is  an  n  -state  vector,  Aum  .  is  an  m  -constant 
— m ,  k  m  — m ,  k  m 

control  vector,  Ay^  ^  is  an  £-observation  vector  and  4>m,  r^, 
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H  ,  D  are  constant  matrices  of  suitable  dimension.  The 
m  m 

objective  of  the  command  generator  tracker  is  to  force  the 
error  state 


4£k 

AX  , 

— m ,  k 

At k  =  [H  D  B] 

4— k 

'  "V  Dm‘ 

-a“k. 

Au  , 

— m ,  k 

to  go  to  zero.  When  the  command  error  is  zero,  the  plant 
states  and  controls  are  said  to  be  tracking  the  star  tra¬ 
jectory  . 


The  star  trajectory  is  defined  to  be  the  state  and 
control  trajectories,  Ax^,  and  Au^  that  satisfy  the  plant 
dynamics 

^  JU 

Ax^+i  «  4>Axk  +  rAHk  +  (A- 1 ) 

and  causes  the  command  error  to  be  zero,  i.e., 


r  *1 

“  * 

A— k 

Ax  . 

— m ,  k 

D  B] 

Jc 

■  IHm  V 

iuk 

4*k_ 

Au  , 

— m ,  k 

(A-2) 


for  all  k>0.  An  excellent  analogy  is  that  the  star  trajectory 
for  linear  systems  is  the  counterpart  of  the  nominal  trajectory 
in  nonlinear  systems.  Since  the  star  trajectory  is  linear,  it 
would  be  expected  that  the  star  trajectory  satisfies  a  linear 
relationship  with  Axm,  Ai^,  and  Aw.  When  the  command  A^  is 
constant  and  the  disturbance  is  constant,  (i.e.,  4>  ,  ♦  ,  Dm 

w  in  in 

are  identity  matrices,  and  B,  Tm,  Hm  are  zero  matrices)  Eqs. 

A-l  and  A-2  can  be  combined  to  form 
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(♦-i)  r 


H 


r  *1 

Ax 

-EAw 

: k 

Au 

Au 

-m 

J 

a.  m 

_  _ 

If  the  quad  partition,  ^  ^  ,  is  invertible,  Ax  and  Au 

can  be  found  using  the  inverse 


: k  ~ 

Ax 

k 

“ 

Au 

<*-I)  r 
H  D 


*  -1 

-EAw 

Au 

-m 

(A-3) 


For  constant  commands  and  disturbances,  the  degenerate  star 
trajectories  are  points  in  Rn  and  Rm  and  are  linearly  related 
to  Aw  and  Au^.  Using  Eq.  A-3  as  a  guide,  the  star  trajec¬ 
tories  when  A^  is  not  a  constant  is  assumed  to  satisfy  the 
linear  relationship 


k 

A*k 

A11  A12  A13 

Ax  . 

— m ,  K 

= 

Au 

-m 

it 

4uk 

A21  A22  A23 

4«k 

The  star  trajectory  can  be  determined  if  solutions  for  to 

A23  exists  in  terms  of  plant  and  model  parameters.  The  feed¬ 
forward  matrices  A^  to  A23  are  assumed  to  be  constant  and 
have  the  appropriate  dimensions  indicated  in  Eq.  A-A. 


A  matrix  equation  for  the  feedforward  matrices  can  be 
found  by  using  the  plant  equations  and  the  command  error  equa¬ 
tion.  Augmenting  the  plant  with  the  command  equation,  Eq.  A-2, 


10A 


t 

[ 


subtracting  [Ax,  0]  from  both  sides  and  substituting  Eq. 

K  •A'T  'A'T  T 

A-4  for  the  vector  (Ax^  Au^  ]  produces. 


A*k+1  "A*k 


(<t>-i)  r 


H  D 


A11  A12  A13 


A21  A22  A23 


Ax  i, 
— m ,  k 

Au  . 
— m ,  k 


(A-5) 


Subtracting  Axj^  from  Axj^+1  in  the  top  row  of  Eq.  A-4  and  aug¬ 
menting  the  alternate  version  of  the  command  equation  produces 


JL  “ 

A-k+l  "A^k 


All^<t'm”1  >  Alirm  A13<VX> 


AX  , 
— m ,  k 


( A-6 ) 


Equating  the  right-hand  expressions  in  Eqs.  A-5  and  A-6  with 
A-m  k*  A^n’  anci  A-k  arbitrary  vectors  produces  the  matrix 
equation  for  the  feedforward  matrices, 


H  D  I  I  A21  A22  A23 


(4>-l)r  An  A12  A13  _  Alirm  A13(<t>w‘I)"E 


(A-7 ) 


The  next  section  presents  two  numerical  procedures  for  cal¬ 
culating  A^  to  A22  in  Ed-  A-7  and  indicates  sufficient  con¬ 
ditions  for  the  existance  of  a  solution. 


NUMERICAL  METHODS 


A. 2 


Two  numerical  procedures  for  solving  the  feedforward 
matrix  equation  are  presented  in  this  section  for  the  case 
when  the  number  of  controls,  m,  is  greater  than  or  equal  to 
the  number  of  commands,  £.  The  first  procedure  is  numerically 
inefficient  but  intuitively  instructive  about  the  conditions 
for  the  existence  of  a  feedforward  matrix  solution.  The 
second  procedure  is  numerically  efficient  and  is  used  to  gen¬ 
erate  the  results  in  Chapter  4. 

The  feedforward  matrix  solution  will  depend  on  a 
relationship  between  the  eigenvalues  of  the  model  and  the 
discrete-time  transmission  zeroes  of  the  discrete  plant. 

Using  a  discrete- time  analogy  to  the  results  in  Ref.  40,  a 
discrete-time  transmission  zero  is  any  complex  number,  A, 
which  satisfies 


rank 


(4>-A  I )  r 


H 


<  n+B. 


(A-8) 


The  transmission  zeroes  can  be  determined  using  the  eigen 
value/eigenvector  relationship  in  Ref.  41  given  by 


4> 

r 

=  A 

I 

o' 

Ji 

i _ 

H 

D 

4 

0 

0 

• 

• 

• 

» 

. 

» 

where  is  the  transmission  zeroes  eigenvector, 

the  right-weighted-pseudoinverse  , 


(A-9) 


Using 
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m  uwj 


<t>  r 

nn 

ni2 

ii 

■  “ 

I  0 

H  D 

n21 

°22 

. 

0  I 

Equation  A-9  can  be  rewritten  as 

"n  £x  *  X  4 


(A-10) 


(A-ll) 


'21 


=  f  t 


The  inverse  of  all  eigenvalues  o  f  nn  which  are  invariant 
under  all  pseudoinverses  in  Eq.  A-10  are  the  transmission 
zeroes.  An  invariant  zero  eigenvalue,  ^ ,  of  means  the 
transmission  zero  (1/A.n  .)  is  at  infinity  as  in  the  classical 
sense  for  scalar  transfer  functions.  When  £  and  m  are  equal, 
the  pseudoinverse  becomes  an  inverse  and  all  eigenvalues  of 
nil  are  transmissi°n  zeroes.  An  instructive  but  numerically 
poor  way  to  compute  the  pseudoinverse  is 


nil  ni2 

(♦-i)  r 

T 

1 

(«f-D  r 

(o-i)  r 

T 

=  Q'1 
P 

Q’1 

T 

n21  n22 

H  D 

{ 

H  D 

H  D 

• 

- 

\ 

“  • 

-  ^ 

/ 

where  Qp  is  the  pseudoinverse  weighting  matrix.  For  the  nu¬ 
merical  results  in  the  Chapter  4,  the  weighted  pseudoinverse 
is  actually  computed  using  singular  value  decomposition 
(Ref.  42). 
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For  the  first  numerical  procedure  one  can  use  to 

compute  the  feedforward  matrices,  assume  that  Tm  is  the  diago 

nalizing  matrix  for  4>  and  T  is  the  diagonalizing  matrix  for 

m  w 

4>w •  Treating  each  matrix  column  partition  in  Eq.  A-7  sepa¬ 
rately,  each  partition  can  be  rewritten  as  follows 


<t>  T 


H  D 


An  T 
11  m 


A01  T 
21  m 


A,  ,T  T_14>  T 
11  m  m  mm 


HT 
m  m 


(A-12) 


»  - 

(♦-i)  r 

A12 

II 

" 

A11  rm 

H  D 

A22 

Dm 

(A-13) 


<t*  r 


H  D 


A13  Tw 


A23  Tw 


-  et» 


■BT 

w 


(A-14) 


If  A  .  denotes  an  eigenvalue  of  <t>  and  A  .  denotes  an  eigen- 
m ,  i  b  m  w ,  l 

value  of  4>w,  then  Eqs.  A-12  and  A-14  can  be  rewritten  as 


(*-x  .i)  r 

m ,  i 

(A11 

Tm>  • 
m  i 

ii 

t - 

o 

H  D 

(A21 

1 

H 

(H  T  )  . 
m  mi 

>  • 

>n 


m 
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(4> -A  .  I)  r  (A,  o  T  )  .  (-ET  )  . 

w,i  13  w'l  w'i 


D  (A23  "Vi 


(-BT  ) . 
w  1 


i  =  l , .  .  .  ,n 


where  (A)^  denotes  the  l  column  of  the  matrix  A.  Comparing 
the  above  two  equations  and  Eq.  A-12,  with  the  definition  of  a 
transmission  zero  in  Eq.  A-8,  it  follows  that  the  feedforward 
matrices  have  a  solution  provided  no  eigenvalue  of  either 
or  4>w  is  equal  to  a  transmission  zero  of  the  plant  and  no 
discrete- time  transmission  zero  is  equal  to  one.  Under  these 
conditions,  the  quad  partition  matrix  on  the  left-hand  side  of 
the  above  feedforward  matrix  equations  have  full  rank  and 
solutions  for  A^  to  A22  are  determined  by  solving  each  equa¬ 
tion  using  singular  value  decomposi tion .  Numerical  procedures 
for  handling  complex  matrices  are  found  in  Ref.  43. 

In  the  second  numerical  procedure,  the  right  pseudo¬ 
inverse  of  the  quad  partition  matrix  is  used  directly  in  Eq . 
A-7  to  produce, 


A11  ‘  nll  A11  “m  +  nl2  Hm 


(A-15 ) 


A13  '  nil  A13  "w  nil  E  '  n 12  B 


(A- 16 ) 


A21  =  n21  A11  *111  +  n22  Hm 


A23  =  n21  A13  ^w  n21  E  ‘  n22  B 


A12  =  X11  A11  rm  +  x12  Dm 


A22  =  X21  A11  rm  +  x22  Dm 

where 


(♦-I)  r 

H  D 


X11 

x  12 

ii 

I 

0 

X21 

x22 

1 - 

o 

- , 

( A- 17 ) 


The  pseudoinverse,  x»  shown  in  Eq .  A- 17  has  the  same  weighting 
matrix  used  in  constructing  n  and  is  similarly  computed  using 
singular  value  decomposition.  The  non- trivial  equations  are 
Eqs.  A-15  and  A-16  which  have  the  form 


X  =  A  X  B  +  C 


where  X  is  unknown.  A  numerical  algorithm  to  solve  Eqs.  A-15 
and  16  is  available  by  generalizing  the  results  in  Ref.  44  and 
using  iterative  residulization .  The  algorithm  has  a  solution 
provided 


f  1 


A 


n,i 


A 


w,  j 


t  1 


where  An  i  are  the  eigenvalues  of  n^.  Using  Eq .  A-ll  and  the 
above  two  eigenvalue  conditions,  Eqs.  A-15,  A-16,  and  A-17 
have  a  solution  provided  no  eigenvalue  of  either  4>m  or  4>w  is 
equal  to  a  transmission  zero  of  the  plant  and  no  discrete- time 
transmission  zeroes  are  equal  to  one.  The  conditions  needed 
to  determine  the  feedforward  matrices  using  either  numerical 
procedure  presented  is  the  same.  The  second  numerical  pro¬ 
cedure  has  the  advantage  of  using  only  real  numbers  (i.e.,  no 
complex  arithmetic  is  needed)  and  the  algorithms  employed  are 
noted  for  their  efficiency. 
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A. 3  PROPORTIONAL- INTEGRAL  CONTROL  ABOUT  THE  STAR  TRAJECTORY 

If  the  plant  shown  in  Eq.  A-l  is  stable,  then  the 
open-loop  control  obtained  by  commanding  Au  with  Au  will 
cause  the  plant  to  asymptotically  track  Aym  when  Ax(0)  and 
Ax"(0)  are  not  equal.  For  the  AV-8A  aircraft,  however,  many 
flight  conditions  are  either  unstable  or  have  poor  perfor¬ 
mance.  The  problem  can  be  rectified  by  using  additional 
control  compensation . 

One  simple  compensation  is  state  feedback  about  the 
star  trajectory 

-J-  JU 

Au^  =  Au^  +  K(CAX^  -  CAx^)  (A-18) 

where  CAx^  indicates  that  full  state  feedback  is  not  nec¬ 
essary.  If  the  error  states, 

AXk  =  a*k  -  ax* 

aak  =  Auk  -  au* 

are  defined,  substituting  Eq.  A-18  into  the  plant  dynamics 
produces 

A-k+l  =  "  rKC)  A^k 

If  K  is  choosen  so  that  the  closed  loop  system  is  asymptoti¬ 
cally  stable  then  the  error  state  Ax,  goes  to  zero  as  k  in- 

*  *  ■"  i 

creases,  Ax^  goes  to  Ax^,  Au^  goes  to  Au ^ ,  and  Ay^  goes  to 
,  k ' 

The  control  law  shown  in  Eq .  A-18  can  be  rearranged 
to  indicate  the  feedback  and  feedforward  gains. 


A-k  '  KCA-k  +  <A2i'kcah  )A*ni,k  +  (A22"K''A12)A-m 


+  (A23-KCA13)Awk 


If  the  plant  dynamics  have  modeling  errors  (almost  always  the 
case)  but  K  keeps  the  closed-loop  plant  stable,  then  perfect 
asymptotic  tracking  will  probably  not  occur  because  the  feed¬ 
forward  matrices  are  incorrect.  Other  problems  in  imple¬ 
menting  the  feedback/feedforward  control  law  occur  because  Aw^ 
is  usually  not  measured  and  noise  disturbs  the  plant  and 
disturbances.  The  last  two  problems  mentioned  are  outside  the 
scope  of  this  work  but  are  discussed  in  Refs.  28  and  45.  The 
disturbance  is  dropped  from  the  rest  of  this  appendix. 


Modeling  error  mismatch  and  feedforward  matrix  in¬ 
accuracies  can  be  handled  using  additional  compensator  states. 
Examples  of  some  current  control  theory  work  in  this  area  for 
continuous-time  problems  are  given  in  Refs.  46  and  47.  Unfor¬ 
tunately  none  of  these  references  demonstrate  the  connection 
between  the  feedforward  matrices  and  compensators  and  no 
optimal  control  designs  are  presented.  The  discrete- time 
feedforward  matrices  presented  in  this  appendix  and  the 
discrete- time  star  trajectories  needed  in  the  optimal  control 
approach  are  new  developments. 


As  an  example  of  the  connection  between  the  feed¬ 
forward  matrices  and  compensators,  consider  the  problem  ad¬ 
dressed  in  Section  A.l  of  tracking  a  constant  command.  Com 
paring  Eqs.  A-3,  A-4,  and  A-17,  the  feedback/feedforward 
control  law  is 


AHk  =  KCA^k  +  (*22  ’  KCx12^A-m 


lf  x22 
to  Au 


and  X|2  are  incorrect  Ayk  will  not  necessarily  be  equal 
in  steady-state.  Consider  next  including  integrator 


■ 
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states  to  form  a  proportional- integral  control  law  as  shown  in 
Ref.  26.  The  control  law  in  incremental  form  is 


A-k  ‘  A-k-l  +  Klc(A^k'A-k-l)  +  K2(A^k-l'AHm) 

If  there  are  modeling  errors,  but  the  proportional  gain  K-^  and 
the  integral  gain  I<2  stabilize  the  plant,  Ay^,  is  guaranteed 
to  approach  Ai^  in  steady-state.  Although  X22  anc*  xi2  are 
needed  in  the  derivation  of  the  optimal  PI  control  law  as 
shown  in  Ref.  26,  they  are  not  needed  in  implementation. 

In  the  AV-8A  command  models  used  in  Chapter  4,  most 

of  the  Ay  commands  are  not  constants  hence  there  is  an  in- 
•Hn 

dication  more  advanced  compensation  is  needed.  In  the  fol¬ 
lowing  paragraphs,  the  discrete- time  PI  control  law  is  re¬ 
derived  including  the  generalized  command  model.  After  the 
derivation  is  complete,  a  discussion  is  presented  indicating 
why  the  PI  control  law  is  usually  adequate  when  using  the 
complex  command  models  presented  in  Chapter  4. 

We  begin  by  defining  the  error  states  and  controls  as 
before  and  include 


avk  =  avk  -  Av" 

where  Av^  is  the  control  difference, 

Au^  =  Au^_  ^  +  AtAv^_i  (A-19) 

JL 

The  ideal  trajectory,  Avj^,  can  be  constructed  as 

avk  =  k  (Au^-Au*^)  =  L  (Ai2[iS.,k-A*.,k.l 

Subtracting  the  ideal  model  response  shown  in  Eq .  A-l  from  the 
plant,  and  adding  and  subtracting  Atv^  from  Eq .  A-19  produces 


]) 


<t>  r 


o 


— k+l 


Hk+1 


0  I 


-k 


Hk 


Atl 


^k 


The  assumption  is  introduced  that  previous  to  k=0 ,  the  command 

Au  is  constant  and  after  k=0 ,  Au  is  constant,  but  at  k=0  the 
— m  — m 

command  instantaneously  changes, 


Au  ,  = 
— m ,  k 

k  =  0, 

1,  2,  .. 

4Vk  = 

4V-> 

k  =  .  .  , 

-2,  -1 

4Vt) 

-  Au.(-) 

?  o 

(A-20) 


In  implementation,  any  time  Au^  changes  the  linear  time- 
invariant  plant  is  considered  to  be  at  k=0  (superposition). 
Equation  A-20  allows  us  to  continue  with  the  derivation  as  if 
Aum  is  constant  but  will  indicate  the  correct  procedure  to 
undergo  if  Aum  ,,  is  different  from  Au  ,  ,  .  If  Au  .is 

slowly  varying  with  respect  to  the  closed-loop  system  time 
constants,  the  control  sequence  will  be  almost  always  optimal. 


The  actual  objective  of  the  optimal  control  is  to 
transfer  the  plant  from  perfectly  tracking  the  model  output 
for  Aum(-)  to  perfectly  tracking  the  model  output  for  Aum(+) 
while  minimizing  the  cost  function 


J 


Qj  Mj  M2' 

1 - 

««l 

<1 

■ _ 

Auk  Avk ] T 

q2  h3 

iuk 

R 

45k 
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The  cost  function  starts  at  k=-l  because  Av_^  is  not  speci 
f ied ,  i  .  e  .  , 

2-1  -  At  <2o  •  M> 


Before  minimizing  the  cost  function,  a  contradiction  in  the 
value  of  AXq  must  first  be  resolved. 

At  k=l  and  k=0,  Eq.  A-A  produces 


"  JL  ’ 

Ax^ 

A11 

A12 

Ax  , 

— m ,  -  i 

* 

Au_  ^ 

. 

A21 

- 

■ 

CN 

CN 
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I - 

> 

£ 

1 

s-/ 

1 _ 

(A-21) 
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and  the  ideal  plant  transient  is 


4>Ax”^  +  T  Au_  ^ 


(A-23 ) 


a’ 

AXq  specified  by  Eq.  A-22  is  not  equal  to  AXq  specified  by 
Eqs.  A-21  and  A-23.  To  eliminate  this  discrepancy  the  re¬ 
quirement  that 


Au__  , 

-m  j  R 


=  0 


(A-24) 
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is  invoked  only  at  index  k.  As  the  system  progresses  with  a 
constant  Au  ,  whenever  Au  ,  is  different  from  Au  ,  , ,  we  go 

back  to  k-1  and  restart  all  variables  we  have  control  over  in 

j,  j. 

time.  For  example,  Axm  ^-1’  ^-m  k-1*  ^-k-1’  anc*  ^-k-1  can 
changed  but  Ax^_^  and  Au^.^  are  immutable  for  casual  systems. 


All  that  happens  then  is  that  the  model  is  simulated 


as 


Ax  ,  =  4>  Ax  .  ,  +  r  Au  . 

— m , k  m  — m , k-1  m  — m,k 


and  Au  ,  can  replace  Au  ,  ,  in  implementation  as  long  as 

casuality  is  not  violated.  For  example,  Ax_^  and  Au_^  have 
already  occurred  when  the  cost  function  starts  to  be  mini¬ 
mized,  and  Aum(+)  becomes  known  at  k=0.  To  stress  the  im¬ 
portance  of  Eq.  A-24,  the  index  k  is  not  dropped  for  Aum  in  the 
following  derivation. 


The  control,  Av^,  which  minimizes  the  cost  function 
under  general  conditions  is  simply  the  solution  to  the  linear 
optimal  regulator, 


4vk  =  -Kj  Axk  -  K2  Auk 

where  and  are  gains  determined  from  the  linear  optimal 
regulator  algebraic  Riccati  equation  solution.  The  control 
Au^  becomes 

Au^  =  Au^_  ^  -  AtK^  Ax^_  ^  -  AtK2  Au^_  ^  (A-25) 

and  does  not  have  the  integral  property. 

The  incremental  form  of  the  discrete  control  which 
has  the  Type  1  property  is 
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Auk  =  AHk-i  “  Cl^A-k"A-k-l^  "  AtC2^H 


Ax 


k-1 


4Sk-i 


(A-26) 


The  question  to  be  answered  is  this:  What  transformation 
allows  the  Type  1  gain  matrices,  and  C2 ,  to  be  expressed  in 
terms  of  the  Type  0  gains  K-^  and  K2  •  The  solution  begins  by 
noting  that 

4Sk  -  4*k.l  *  <♦-!)  4*k-l  *  4“k-l 
and  substituting  this  into  Eq.  A-26  producing, 


4Sk  =  4Hk-! 


[  A  tC2 ] 


(0-1)  r 

4-k-l 

H  D 

4-k-l 

( A-27 ) 


Equation  A-27  is  equivalent  to  Eq .  A-25  if 


j 


4 


(0-1)  r 


[  A tC2 ] 


[AtK1  AtK2l 


(A-28 ) 


Equation  A-28  is  solved  only  for  £<^m,  because  for  £>m  the  cost 
function  is  almost  always  not  finite. 


The  left  weighted  pseudoinverse  of  the  quad  parti¬ 
tion  transpose  is  used  in  Eq.  A-28  so  that  the  unknown  gain 
magnitudes  have  the  best  possible  alignment.  If  the  same 
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V*  ~  * 


pseudoinverse  weighting  matrix  is  used  in  Eq .  A-28  as  in  Eq. 
A-17 ,  the  gains  are  determined  from 


[C1  AtC2] 


(AtK1  AtK2) 


X 


21 


X 


22 


If  £<m,  Eq.  A-28  is  not  necessarily  an  equality  after  and 
C2  are  calculated,  and  care  should  be  taken  when  using  the 
substitution . 


With  and  C2  determined,  the  PI  command  generator 
tracker  becomes , 


Au^  —  Au^ _ ^  ^ ^  k ” ^-k -*  1  ^  *  A tC2  \  (H  D) 


a£k-l 


AHk-l 


<H  D  ) 
m  m 


A^m,k-1 


AU  b 
— m ,  k 


(C1AU  +  A21]  [Ax^  -  , k- 1 1 


( A-29 ) 


where  Ax  and  Au  have  been  eliminated  using  Eq.  A-4  and  the 
following  substitution 
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A11  A12 

Ax  , -Ax  ,  , 

—ro ,  k  — m ,  k- 1 

I] 

*  * 

aHk-^U-l 

=  [C:  I] 

A21  A22 

0 

Equation  A-29  and  A-25  have  the  same  transient  if  £= m  but 
markedly  different  responses  to  sensor  noise  and  plant  varia¬ 
tion  . 


The  control  law  shown  in  Eq .  A-29  depends  on  A^  and 
A22  but  not  on  A^  and  ^22'  m°deling  errors  that  impact 

on  A^^  and  degrade  the  performance  of  the  control  law. 

There  are  basically  two  types  of  command  models  used  in  the 
AV-8A  control  law,  models  which  generate  ramps 

Ax  , ,  =  AXL,  h  +  A  tAu  , 

—m ,  kl  — m ,  k  — m ,  k 

and  model  states  that  are  dynamically  stable.  For  the  stable 
model  states,  the  term  Ax  , -Ax  .  ,  in  Eq.  A-29  has  an  ini- 

ID  )  K  — m  y  K  “  JL 

tial  contribution  to  the  control,  Au^ ,  when  Au^  ^  steps  but  as 

k  increases  Ax  ,  -Ax  ,  ,  approaches  zero  and  the  proportional 

— m,k  — m , k - 1 

integral  property  of  the  control  law  shown  in  Eq .  A-29  forces 
Ayk  to  Aym  The  CGT-PI  performs  well  for  sufficiently 
asymptotically  stable  dynamic  command  models. 

The  same  argument  applies  for  the  ramp  command  model. 
For  example,  when  the  pilot  inputs  a  constant  acceleration 
command,  Aa,  the  controls  move  to  gain  speed.  Eventually,  the 
pilot  will  stop  accelerating  the  vehicle,  Ax^  ^ be¬ 
comes  zero  and  the  control  law  will  force  Ay^  to  Aym  ^  re¬ 
moving  any  errors  caused  by  modeling  errors  in  A^  and  A2^.  1° 
many  instances  with  the  AV-8A  command  models,  A^  and  A2J  for 
ramp/step  models  are  very  sparse  and  any  numbers  at  all  are 
either  1  or  can  be  investigated  to  determine  their  origin  in 
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the  original  plant  dynamics.  Explicit  expressions  for 
and  A22  can  then  be  programmed  in  implementation  with  no 
numerical  inaccuracy.  An  example  of  sparse  A^  and  A21  are 
shown  in  Section  3.2.2. 

The  control  law  which  removes  the  knowledge  require¬ 
ment  for  Aj^  and  A2^  when  using  ramp  models  is,  of  course,  a 
Type  2  control  law  of  the  form 

A-k  =  A-k-l  +  C3(A-k'A^k-l)  +  AtC4A£-k-l 
+  AtC5(A^k-l"A^n,k-l) 

A^k  =  A^k-1  +  At^A^k-l”A^m,k-l } 

In  the  Type  2  control  law  AtC^A|k_^  estimates  the  constant 

^C1A11+A21^  ^A-m  k“A^m  k-1^’  A  area  of  research 

would  be  to  combine  the  star  trajectory,  feedforward  matrices, 
and  Davison's  servocompensator ,  (Ref.  46)  into  a  unified,  de¬ 
sign  approach  using  optimal  control  for  discrete  systems. 


A. 4  EIGENVALUE/EIGENVECTOR  PLACEMENT 

A  new  result  in  direct  digital  control  design  theory 

is  presented  in  this  section.  In  the  CGT  control  law,  using 

simple  state  feedback,  consider  the  case  where  Ax  has  the 

— m 

same  number  of  states  as  Ax,  has  desirable  eigenvalues  and 
eigenvectors  and  l~m  is  a  desirable  control  input  matrix.  The 
control  law  using  full  state  feedback  and  dropping  the  dis¬ 
turbance  term  is 

A— k  =  ^-k  +  (A21'KA11  )A^m,k  +  ( A22_KA12  )AHm 
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Assume  the  feedback  gain  is  chosen  so  that 

K  =  A21  Al\  (A-30) 

The  control  law  becomes 

Au^  =  KAxk  +  (A22  -  A22  A^J  Ai2>Au m  (A-31) 

Command  generator  tracking  is  obtained  without  feeding  forward 
Ax  ,  .  The  feedback  gain  shown  in  Eq.  A-30  structurally 

TO  9  K 

changes  the  plant  to  behave  like  the  model.  Equation  A-30  can 
be  shown  to  be  the  discrete  analog  of  an  eigenvalue/eigenvector 
placement  procedure  discussed  in  Ref.  48. 

The  structural  change  caused  by  using  K  in  Eq .  A-30 
for  feedback  is  clearly  evident  by  comparing  the  transfer 
function  of  the  model  with  the  transfer  function  of  the  plant, 
using  Eq.  A-31.  Substituting  Eq.  A-31  into  the  plant  produces 

A-k+l  =  A11  ^m  A11  A^k  +  r(A22  "  A21  A11  A12)A-id 


The  first  structural  change  noted  is  that  the  eigenvalues  of 
the  plant  closed-loop  matrix. 


<t>  +  TK 


4> 


m 


have  been  changed  to  the  eigenvalues  of  the  model,  4>m<  Noting 
that  (D  is  assumed  to  be  zero), 


r  a 


21 


(4>-I)Au 


r  A22  =  A11  Fm  '  (4>_I)A12 

the  transfer  function  of  the  closed  plant  becomes  after  some 
manipulation , 
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4Zk  =  H(zl  -  AU  *m  A’})'1  [A1;irm 


+  (zl  -  An  4>m  Aj^Aj^j  +  A12<  l-zl )  liUj! 


Using  the  following  information, 


(zl  -  An  4>m  A1J)  1  -  A11(zl-«t>m)  1  A]_J 


( 1  -  z  1 )  u  m  =  0 


H  A 


11 


=  H 


m 


the  transfer  function  reduces  to 

rn,  +  DJ4Hm  (A-32) 

Comparing  the  transfer  function  in  Eq.  A-32  with  the  transfer 
function  of  the  model 

Av  ,  =  [Hm(zl-4>m)'1  Tm  +  D  ]Au 
■4n ,  k  m  m  m  m  — m 

it  follows  that  the  gain  in  Eq .  A-30  structurally  changes  the 

eigenvectors  of  the  plant  so  that  Ay^  is  equal  to  A£m  ^  for 

constant  Au  when  the  release  condition  for  Ax  (0)  satisfies. 

— m  — m 

Ax  =  A, }  Ax  -  A, }  A, 0  Au 
-m,o  11  -o  11  12  — m ,  o 

The  eigenvalue/eigenvector  placement  procedure  is  demonstrated 
with  the  V/STOL  aircraft  using  a  PI  design  in  Chapter  3. 
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APPENDIX  B 

COMMAND  AUGMENTATION  MODELS 


Six  command  models  shown  in  Table  8  are  used  to 
control  the  AV-8A  at  various  flight  conditions.  Each  command 
model  has  advantages  and  disadvantages  depending  on  the  appli¬ 
cation.  This  appendix  presents  the  nonlinear  relationships 
needed  to  construct  the  command  vector  y  for  each  of  the 
command  models  from  the  aircraft  body  states.  The  nonlinear 
relationships  are  linearized  so  that  the  feedback  gains  and 
feedforward  matrices  can  be  determined.  Also  shown  are  the 
linear  and  nonlinear  command  model  dynamics  used  for  model 
following  and  eigenvalue/eigenvector  placement. 


B.l  LINEAR  AND  NONLINEAR  RELATIONSHIPS 

The  command  augmentation  system  design  method  used  in 
Chapter  4  requires  the  nonlinear  version  of  the  command  model 
and  the  aircraft  command  vector,  y,  for  implementation  and 
requires  the  linearized  model  and  command  vector  for  design. 
The  first  four  command  models  in  Table  8  are  linear  and  do  not 
have  to  be  linearized.  These  models  consist  of  simple  inte¬ 
grators  for  the  ramp  commands  and  direct  feedthrough  of  the 

model  control,  u  ,  for  the  step  commands. 

— m 

For  example,  the  enroute  ramp/step  velocity  command 
system  is  represented  as 
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k+1 


The  control  vector  [a 


no  m>V7  m>4<  J  is  the  Pilot's 

DR,m  Z,ro  m  CR,m  m 

input  vector,  Aum<  The  other  ramp/step  models  in  Table  8  have 
structures  similar  to  Eqs.  B-l  and  B-2  for  both  implementation 
and  design,  except  for  the  air  combat  maneuvering  (ACM)  model. 
The  nonlinear  ACM  model  for  the  AV-8A  at  87.5  m/s  (170  kt )  in 
continuous  time  is 
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(B-3 ) 
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(B-4 ) 
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The  complete  ACM  linearized  command  model  becomes 
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The  discrete  equivalent  of  the  linearized  ACM  command  vector  is 


formed  using  results  shown  in  Table  5.  Note  that  a^  and  am  are 


not  the  conventional  specific  contact  force  components,  but  do 
include  the  effects  of  gravity.  Hence  V,  a^T ,  and  an  completely 


describe  the  earth-relative  acceleration  of  the  aircraft. 


The  aircraft  states  that  are  commanded  by  the  CGT-PI 
control  to  follow  ym  form  the  aircraft  command  hold  vector,  y 
(Table  8).  Many  of  the  elements  in  the  y  vectors  have  non¬ 
linear  relationships  with  the  aircraft  body  states  used  in 


T  _ 


x  =[xz0uqwvrp<j>»|iy] 
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design.  In  order  to  present  the  nonlinear  and  linearized 
relationships  we  find  it  convenient  to  place  the  different 
hold  vectors  into  four  groups. 

The  aircraft  velocity  in  wind  axis  is  the  first  group 
and  is  related  to  the  body-axis  velocity  by 


V 

7u2  + 

2  2 
v  +  w 

p 
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tan  ^ 

(v/7 u2  +  w2 ) 
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„  -1 
tan 

(w/u) 

_  _ 

_ 

_ 

V  is  the  velocity  magnitude,  p  is  the  sideslip  angle  and  a  is 
the  angle  of  attack.  The  velocity  of  the  vehicle  expressed  in 
inertial  axis  (a  nonrotating,  flat  earth  is  assumed)  helps 
form  the  second  group  and  is  expressed  as 
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The  body  to  inertial  transformation,  Hg  ,  is  given  by 
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COS  9 

where  roll ,  $  , 

pitch ,  0  , 

and 

yaw ,  <l> , 

are 

the 

body 

Euler 

angles.  The  second  group  is  composed  of  the  vehicle  velocity 
expressed  in  local-level  axes. 
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In  the  models  in  Chapter  4,  the  vehicle  is  always  assumed  to 
be  heading  north  which  makes  the  nominal  yaw  angle  zero. 

The  third  group  is  the  inertial  velocity  magnitude 
and  inertial  velocity  Euler  angles  that  are  related  to  the 
inertial  velocities  by 


V 

h2  +  v2  +  v2 

x  y  z 

£ 

* 

arc  cos  <Vxa/v^  +  V^) 

Y 

arc  sin  (V z/Jv'2  +  V2  +  V*) 

. 

.  _ 

£  is  the  velocity  heading  angle  and  y  is  the  flight  path 
angle.  The  last  group  consists  of  the  wind-axis  roll  angle, 
4>y,  used  in  the  air  combat  maneuvering  command  vector.  The 
angle,  <$>y,  completes  the  transformations  between  the  different 
axis  systems  as  shown  in  Fig.  B-l.  The  angle,  4>y,  is  calcu¬ 
lated  using  one  of  the  elements  in  the  transformation 


H  v  ( ^  V  ^ 


1  0  0 
0  cos  4>y  sin 
0  -sin  cos  4>y 


=  Ho  ( a  ,  p  )  H?(<t>,6,40  Hj(y,U 


( B-5 ) 
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Figure  B-l  Reference  Axis  Transformations 

(Arrows  Indicate  Right-Hand  Rotation) 


The  body  to  wind  axis  transformation  is  the  transpose  of 


cos  a  0  -sin  a  cos  p  sin  p  0 


Hjj(a,p)  =010 

sin  a  0  cos  a 


-sin  p  cos  p  0 
0  0  1 


while  the  velocity  to  inertial  axis  transformation  is  the 
transpose  of 


cos 

Y 

0 

-sin 

■ 

Y 

cos 

i 

sin 

t 

0 

0 

1 

0 

-  sin 

i 

cos 

t 

0 

sin 

Y 

0 

cos 

Y 

0 

0 

1 

Using  available  aircraft  sensors  (IMU,  barometric  altimeter, 
gyros,  angle  of  attack  and  sideslip  sensors)  and  ground  radar 
when  the  aircraft  is  in  the  terminal  area,  Kalman  filtering 
can  be  used  to  provide  estimates  of  the  commands  used  in  all 
the  command  vectors,  y.  An  example  of  a  Kalman  filter  design 
for  a  simple  velocity  command  system  is  shown  in  Ref.  26. 


The  linearized  version  of  each  of  the  four  command 

vector  groups  are  shown  next.  Each  of  the  linearized  vectors 

are  to  be  related  to  the  linearized  body-axis  state  vector, 

T 

Ax  .  The  linearized  wind-axis  velocity  vector  is 


AV 

Au 

A0 

=  Jw1(Vo'Po)  HB<(W 

Av 

Aor 

Aw 

„  m 

Jv  is  a  diagonal  matrix  which  has  elements  1,  V  ,  and  VQ  cos 
0Q.  The  inertial  velocity  vector  varies  with  the  Euler 
angles,  hence  the  linearized  inertial  velocity  vector  includes 
the  sensitivity  to  A0  ,  A0  ,  and  At|> : 


AVx' 

Au 

A<t> 

AV 

=  Hp 

Av 

-  vD  Ld 

A6 

y 

AV 

Z 

B  ,o 

Aw 

B  ,  o  B  ,o  B  ,  o 

At)) 

»  a 

.  . 

»  « 
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v_  is  the  cross-product  equivalent  matrix  of  the  body  axis 
15,0 

velocities  (no  wind) 


0  -w^  v 

o  o 


Ld  transforms  Euler  angle  perturbations  to  body-axis  co¬ 
ts  ,  o 

ordinates  and  is  given  by 


sin  $  tan  8 
o  o 

cos  <b 


cos  <t>  tan 

o  o 

-sin  6 


sin  <b  sec  8^ 
o  o 


cos  <J>  sec 

o  o 


Since  the  nominal  yaw  angle  is  always  zero  in  the  design,  the 
perturbation  local-level  velocity  is  the  same  as  the  pertur¬ 
bation  inertial  velocity.  The  linearized  version  of  the 
inertial  velocity  magnitude  and  inertial  Euler  angles  is 
accomplished  by  first  relating  them  to  the  linearized  inertial 
velocity  vector, 


AC  = 


cos  yo  cos  iQ 
~sin 

Vo  cos  *o 


cos  Yo  sin  C0 

005  ^o 

V  cos  y 
o  1  o 


-sin  yo  AVx 


-sin  y„  cos  £ 


sin  yo  sin  £, 


-cos  y. 


Equation  B-6  is  substituted  into  Eq.  B-7  to  relate  AV,  a£ ,  and 
Ay  to  the  perturbation  body-axis  states.  The  final  and  rather 
involved  linearization  is  for  <p  .  Using  Eq .  B-5,  with  A<J>  the 

v  v  v 

first  element  in  the  three  state  vector,  AfcLy.  the  lineariza¬ 
tion  is 


AuX  '  Hp  (I  +  H?  „  Lv  j"*  V,.  ^  hY  „  Hp  1LR  „ 

*-W  B ,  o  I,o  V,o  V,o  V,o  I,o  B,o  B,o 


A  0 
A6 
At|) 


II  1  ~ ^  I  T  ~  ^  1 

B,olLW,o  jW,o  "1,0  V,o  v ,o  w , o  “B , o 


Au 

Av 

Aw 


The  matrices  which  have  not  already  been  identified  in  the 
expression  for  A^y  are 


V. 


V  ,o 


0  0 
0  0 


0 

-v. 


0  V  0 


L  J  ”  ^  = 

v ,oJV  ,o 


sin  £ 


V  COS  Y 

o  1  o 


-cos 

vo  cos  70 


V 


and 
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sm  a 


0 


o 


LWJW 


Using  the  equations  given  in  this  appendix,  each  of 
the  six  command  hold  vectors  can  be  placed  in  the  form 

Ay  =  H(Xq )  Ax 

The  matrix  D  is  not  required  for  any  of  the  command  vectors. 
For  the  command  vectors  in  Table  8  and  the  AV-8A  dynamics  at 
the  different  flight  conditions,  the  matrix 


v . 


0  0 


•cos  a 


V 


V  cos  p 
o  o 


rank 


F(x0) 


G(xq) 


H(xq  )  0 


n  +  £ 


always  has  full  rank.  When  the  quad  partition  matrix  has  full 
rank,  the  PI  feedback  gains  can  always  be  computed  and  the 
first  of  two  conditions  needed  to  compute  the  feedforward 
gains  is  satisfied. 


B . 2  MODEL  FOLLOWING  DESIGN 

While  most  of  command  models  used  in  the  AV-8A  con¬ 
trol  designs  are  simple  integrators,  a  dynamic  command  model 
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“t#  in  |i  MUM 


m 


is  also  used  for  comparative  purposes  from  both  a  design  and 
implementation  point  of  view.  This  section  describes  the 
dynamic  command  model  and  discusses  reasons  for  the  choice  of 
the  numbers  used  in  the  model. 


The  longitudinal  dynamics  in  the  command  model  are 
desirably  separated  from  the  lateral -directional  dynamics  and 
have  the  form 


Az 

m 

A0 

m 

A^DR,m 

= 
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^m 

AV  m 
z  ,  m 
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(B-8) 


The  down-range  velocity  of  the  vehicle,  height,  and  pitch 
angle  command  model  dynamics  are  decoupled.  The  f^  elements 
are  chosen  to  achieve  acceptable  closed-loop  properties  while 
the  g^  elements  are  chosen  to  achieve  response- to- input  sen¬ 
sitivities  compatible  with  the  appropriate  MIL-F-83300  re¬ 
quirements  (Ref.  15).  The  elements  f^  and  g^  are  -0.5  (2 
second  time  constant)  and  5.0  respectively.  The  numbers 
are  taken  from  Ref.  19  where  a  model  following  system  is 
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constructed  similiar  to  Eq .  B-8  for  the  X-22A  V/STOL  aircraft. 
The  pitch  angle  and  pitch  rate  time  constants  are  chosen  to  be 
0.625  seconds  and  0.4  seconds  respectively  using  Ref.  29  as  a 
guide.  The  control  element  g2  is  set  at  25.06  in  order  to 
meet  the  pitch  attitude  response  in  MIL-F-83300, 

A6 

1.18(3)  <_  -  <  7.87(20)  deg/cm  (deg/inch)  in  one  second 

°q,m 


The  g2  chosen  gives  a  A6/A6^  m  of  approximately  1.57  deg/cm 
(4.0  deg/inch)  for  the  command  model.  The  MIL-F-83300  re¬ 
quirement  on  velocity  response  for  the  vertical  axis  is 

AV 

12(100)  £  2-m ■  £  90(750)  mpm/cm  (fpm/inch)  in  one  second 

°Vz  ,m 

Setting  f^  to  -1.5  (0.667  time  constant)  for  good  response 

characteristics,  a  g~  of  11.26  provides  a  AV  /A6..  of 

J  z  j  in  v  )  z  ^  m 

42  mpm/cm  (350  fpm/inch)  in  one  second  for  the  command  model. 

In  the  lateral -directional  command  model,  a  turn- 
coordinated  system  is  considered  desirable.  Only  two  of  the 
states,  Av  and  A<J) ,  are  to  be  commanded  and  lateral-directional 
coupling  is  not  particularly  desirable  hence  the  model  becomes 
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g4  g5 

Ar m 
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(B-9 ) 


If  roll  angle  and  roll  rate,  A®  and  Ap  ,  are  included  in  the 

m  m 

model  without  coupling  to  the  above  states  and  controls,  the 
appropriate  columns  in  the  feedforward  matrices  (see  Appendix 


134 


A)  are  zero.  The  ideal  dynamic  model  does  not  have  to  have 
the  same  number  of  states  as  the  aircraft  model.  The  element 
f^  in  Fm  in  Eq .  B-9  is  set  approximately  equal  to  its  counter¬ 
part  in  the  plant  dynamics  (-2.9066  in  F,  -2.9  in  F  ). 

A  desirable  complex  pair  for  Av  and  Ar  has  a  nat- 

mm 

ural  frequency  of  2.07  rad/sec  and  a  damping  ratio  of  0.6. 

These  choices  place  the  models  lateral  model  approximately  in 

the  middle  of  the  optimum  frequency  and  damping  window  shown 

in  Fig.  43.  Assuming  f^  and  fg  are  equal  the  values  become 

-1.242  and  0.9456  for  f^  and  fy  ,  respectively.  Two  conditions 

are  used  to  determine  g^  and  gg .  First,  when  Afi^  c  is 

stepped,  the  steady-state  value  of  Av  should  be  zero  and 

m 

second,  the  yaw  response  should  satisfy  the  requirement  in 
MIL-F-83300  given  by 

2.36(6.0)  <_  —  <_  9.06(23)  deg/cm  (deg/inch)  in  one  second 

<l»  ,c 

Using  a  A*|»/A6  •  of  2.53  deg/cm  (6.43  deg/inch)  the  values  for 
4*  9  c 

g^  and  g^  are  23.2  and  9.936  respectively.  When  A6v  c  is 
commanded  to  2.54  cm  (1  inch),  the  steady-state  value  of  Av^ 
is  to  be  3.05  m/s  (10  fps)  with  Arm  in  steady-state  at  zero. 
The  values  for  g^  and  gy  become  12.42  and  -9.456  respectively. 
The  pilot  command,  A6v  ,  would  normally  be  considered  a  pedal 
control  and  is  designed  not  to  be  used  by  the  pilot  to  execute 
a  coordinated  turn. 


B . 3  EIGENVALUE/EIGENVECTOR  PLACEMENT  MODEL 

This  section  discusses  how  the  longitudinal 
continuous- time  eigenvalue/eigenvector  placement  (EEP)  model 
used  in  Section  4.3.4  is  constructed  and  indicates  the  nu¬ 
merical  values  in  the  model.  The  eigenvalue/eigenvector 
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placement  model  requires  a  structure  completely  different  from 
the  feedforward  gain  model  presented  in  the  previous  section. 


The  EEP  model  must  have  all  eigenvalues  in  the  open 
left-half  complex  plane.  The  plant  eigenvalues  will  be  changed 
by  feedback  to  the  EEP  model  eigenvalues  and  asymptotic  track¬ 
ing  requires  all  closed-loop  eigenvalues  to  be  stable.  The 
EEP  model  modes  should  be  structurally  similar  to  the  modes  in 
the  original  plant  in  order  to  avoid  large  gains  and  sensi¬ 
tivity  to  modeling  errors.  For  example,  complex  roots  in  the 
aircraft  open-loop  dynamics  such  as  Dutch  roll  and  short 
period  generally  should  also  appear  in  the  EEP  model  but  at 
more  desirable  locations.  The  EEP  model  to  be  used  in  Chapter 
4  is  compared  with  an  optimal  control  proportional - integral 
control  design,  hence,  additional  compensator  states  are 
adjoined  to  the  usual  aircraft  body-axis  states  in  the  EEP 
model  to  obtained  the  Type  1  property.  The  compensator  states 
are  the  same  as  the  control  states  adjoined  to  the  plant  which 
allow  the  control  rate  to  be  weighted  in  the  optimal  control 
cost  function  (Table  4). 

The  EEP  model  in  continuous-time  has  the  form 


Ax 
— m 

Au 
— m 


F  G 
m  m 


Ax 

m 

Au 
— m 


No  control  input  to  the  model  is  required  since  the  EEP  model 
is  not  to  be  used  for  control  interconnect  design.  The  EEP 
model  is  used  to  design  a  PI  control  law  for  low  speed  flight 
using  the  longitudinal  dynamics.  The  three  available  AV-8A 
longitudinal  controls  at  hover  can  be  used  to  force  three 
outputs  of  the  aircraft  to  have  transients  exactly  like  the 
EEP  model  and  at  the  same  time  force  all  the  eigenvalues  of 
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the  plant  to  be  those  of  the  EEP  model.  The  three  outputs 

chosen  are  Ax  ,  Az  ,  and  A 6  of  the  following  ideal  model, 
m  m  m 
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The  feedback  gains  determined  using  the  discrete  version  of 
Eq .  B-10  will  have  the  same  form  as  Eq .  A-25,  and  will  require 
the  same  gain  transformation  using  the  command  vector  H(xq) 
matrix  as  shown  in  Eq .  A-28.  After  the  transformation,  the 
EEP  gains  will  be  and  C2  in  Eq .  A-29  just  as  in  the  optimal 
control  design.  The  control  modes  Ty,  t  ,  and  tw  in  Eq .  B-10 
will  be  the  integrator  eigenvalues  after  the  transformation 
shown  in  Eq .  A-28.  The  EEP  model  dynamics  have  the  velocity 
and  position  states  Ax  and  Au  decoupled  from  the  rest  of  the 
states.  Decoupled  velocit3T  control  is  known  to  obtain  favor¬ 
able  pilot  ratings  (Ref.  3).  The  element  f  .  is  also  set  to 

CjW 

zero  so  that  a  Az  command  does  not  change  pitch  angle  in  the 
transient  response. 

The  rest  of  the  f  and  x  elements  are  chosen  so  that 
the  eigenvalues  of  the  EEP  model  are  approximately  the  same  as 
the  corresponding  optimal  control  design  eigenvalues  subject 
to  the  eigenvector  constraints  dictated  in  Eq .  B-10.  This 
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choice  satisfies  the  EEP  model  requirements  previously  men¬ 
tioned  and  makes  the  comparison  between  EEP  model  design  and 
optimal  control  design  easier.  Table  B-l  shows  the  closed- 
loop  optimal  control  eigenvalues,  the  EEP  model  eigenvalues, 
and  the  EEP  model  parameters  that  generated  the  results.  The 
parameters  g  ,  g  ,  and  g  have  little  effect  on  C,  and  C0  and 
are  chosen  to  be  2.5,  30.03,  and  26.86  respectively  using 
MIL-F-83300  and  the  results  in  the  previous  section.  More 
discussion  comparing  EEP  and  optimal  control  is  given  in 
Chapter  A. 

TABLE  B-l 

COMPARISON  BETWEEN  OPTIMAL  CONTROL  AND 
EIGENVALUE/EIGENVECTOR  PLACEMENT  CLOSED- LOOP  EIGENVALUES 


T-2912 


DYNAMIC 

MODE 

OPTIMAL  CONTROL 

EIGENVALUE/EIGENVECTOR 

PLACEMENT 

EIGENVALUE/EIGENVECTOR 
MODEL  ELEMENTS 

NATURAL  damp INC  TIME 

FREQUENCY,  CONSTANT, 

rad/eec  RATIO  gec 

NATURAL  r,»MPTMr  TIME 

FREQUENCY ,  CONSTANT , 

rad/sec  kahu  sec 

Aq-Aw 

1.68  0.569 

1.77  0.614 

/fq  =  -2 . 181  \ 

f  =  4.96  / 

AX-jAx 

0.694  0.687 

/1.43\ 

/fx  =  *0 . 35\ 

l 1 . 25/ 

U  -  *«■•  I 

Az- JAz 

0.635  0.689 

li.iJ 

Au 

2.02 

2.00 

fu  =  -1.2 

A6 

0.555 

0.682 

fe  =  -3.147 

JA6 

1 . 34 

1.33 

,q  =  -0.75 

APPENDIX  C 

AIRCRAFT  AERODYNAMIC  AND  CONTROL  SYSTEM  MODEL 


The  aerodynamic  and  mass  data  used  in  this  report  are 
discussed  in  this  appendix.  A  simple  ground-effect  model  is 
also  described,  and  a  description  of  the  AV-8A  control  system 
is  included. 


C.l  AERODYNAMIC  AND  MASS  DATA 

The  basic  source  of  aerodynamic  data  for  the  AV-8A 
V/STOL  aircraft  is  Ref.  6.  The  aircraft  size  and  mass  data 
used  in  this  report  are  given  in  Table  C-l.  The  aircraft 
model  used  in  this  report  included  independent  effects  of  ten 
control  variables,  as  listed  in  Table  C-2.  Note  that  on  the 
actual  AV-8A,  some  of  these  are  mechanically  interconnected 
as  discussed  later  in  this  appendix. 

TABLE  C-l 

CHARACTERISTICS  OF  THE  AV-8A  HARRIER 


T-3068 


Mass  .  m 

410.0 

slugs 

5983. 

kg 

Reference  Area,  S 

201.0 

ft2 

18.67 

i 

Mean  Aerodynamic  Chord,  c 

7.95 

ft 

2.42 

m 

Wing  Span,  b 

25.25 

ft 

7.70 

m 

Center  of  Gravity  Location,  x  0 . 09  c 

eg 


TABLE  C-2 

CONTROL  VARIABLE  RANGES 


T-3069 


Throttle  ,  6^-, 

105% 

to 

0% 

Stabilator, 

11.3 

to 

-10.0 

deg 

Nozzle  Angle ,  6.  . 

8  J 

98.5 

to 

0.0 

deg 

Front  Pitch  Jet,  <5pp 

7.23 

to 

0.0 

9 

in" 

Rear  Pitch  Jet,  6 

8.68 

to 

0.0 

9 

in" 

Aileron,  S 

12.0 

to 

-12.0 

deg 

Right  Roll  Jet, 

4.2 

to 

-4.2 

in2 

Left  Roll  Jet,  6^ 

4.2 

to 

-4.2 

in2 

Rudder,  6r 

15.0 

to 

-15.0 

deg 

Vav  Jet ,  6 

3.5 

to 

-3.5 

9 

in" 

yj 

The  linearized  aircraft  model  used  in  this  study  is 
based  on  the  complete  six-degree-of- freedom  nonlinear  aircraft 
equations  of  motion. 


hb  -B 

(C-l) 

^B 

=  lb2  — b 

(C-2) 

=  (£B+V/m  + 

UB  .1 

Hi*i  -  “B^B 

(C-3) 

•  I 
— B 

= 

T‘1~I,  1 

XB  U,B1B-B 

(C-A) 

The  following  paragraphs  describe  the  functional  relations 
used  to  calculate  the  aerodynamic  and  thrust  forces  and 
moments.  The  data  which  is  used  in  these  relationships  is 
given  in  Ref.  6. 


i 


The  thrust  forces  and  moments  and  those  due  to  the 
reaction  jets  are  not  direct  functions  of  dynamic  pressure, 
and  they  are  calculated,  for  the  purposes  of  this  report,  for 
sea-level  standard  atmosphere  conditions.  Figure  C-l  illus¬ 
trates  the  point  of  application  of  all  these  forces.  The 
reaction  jet  forces  are  calculated  first,  and  they  are  a  func¬ 
tion  of  the  pressure  ratio  in  the  reaction  ducting 

£rcs  =  £rcs(6ni’  6ej’  6atot) 

where 

6atot  =  6fp  +  6rp  +  I6lr|  +  |  6rr 

The  absolute  value  signs  are  necessary  because  a  "negative  jet 
area"  is  used  to  indicate  a  thrust  in  the  opposite  direction 
for  the  bi-directional  jets.  The  reaction  jet  duct  is  shut 


YAW 


(C-6) 


n  >i?M4 


Main  Nozzle  RCS  Jet  Location 
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Figure  C-l 


¥ 


off  for  nozzle  angles  below  20  deg,  which  is  the  only  effect 
of  6q j ,  and  the  duct  pressure  increases  with  engine  speed. 

The  effect  of  S^-pQ-p  is  to  decrease  the  duct  pressure  as  the 
reaction  jets  open.  The  magnitude  of  the  thrust  at  each 
reaction  jet  is  a  function  of  the  duct  pressure  ratio  and  the 
reaction  jet  opening 


TFP  =  TFP(6FP’ 

eRCS) 

(C-7 ) 

trp  =  TRP(6RP’ 

ERCS} 

(C-8) 

tlr  =  TLR(6LR’ 

eRCS) 

(C-9) 

trr  =  TRR(6RR’ 

eRCS) 

(C-10) 

YAW  =  TYAW(6YAW’  eRCS) 

(C-ll) 

The  total  reaction  control  system  mass  flow  can  now  be  cal¬ 
culated 


mRCS  ~  mRCS(TRCS’  6 

Nl} 

(012 ) 

where 

trcs  =  tfp  +  trp  +  |tlr|  +  |TRR  | 

+  |tyaw| 

(013) 

The  absolute  values  are  necessary  for 
jets  which  are  bi-dimentional . 

those  reaction 

control 

The  individual  reaction  jet  force  vectors  in  body 
axes  can  now  be  calculated. 


-FP 


Tpp  sin ( 0 pp ) 


"Tpp  cos(0pp) 


(014) 
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TrP  sln<eRp) 


-trp  cos(eRp) 


(C-15 ) 


where  0pp  and  0Rp  are  constants.  The  actual  value  of  0pp  is 
negative . 


(C-16) 


(C-17 ) 


T 

-YAW 


(C-18) 


The  main  nozzle  thrust  depends  primarily  on  throttle, 
reaction  control  mass  flow  and  forward  speed,  with  mR£g  causing 
a  decrease  in  thrust.  The  effect  of  nozzle  angle  on  the  thrust 
magnitude  is  small,  since  it  only  affects  the  rear  nozzle  splay 
angle . 


Tp  =  >  mRCS ’ 


(C-19 ) 


L- 


TR  "  TR(6N1*  60j’  mRCS ’  M)  (C” 2°) 

The  projection  of  the  thrust  magnitude  on  the  body  axes  is  a 
function  of  nozzle  angle. 


Tp  cos(6t) 
0 

“Tp  sin(8  t ) 


(C- 21) 


tr  cos(et) 

o 

-Tr  sin(6 1 ) 


where,  as  shown  in  Fig.  C-l 


(C-22 ) 


et  =  1.5  deg  +  60j  (C-23 ) 

The  last  engine- induced  force  is  ram  drag.  This  force  is 
aligned  with  the  velocity  vector  and  is  applied  to  the  vehicle 
at  the  engine  intakes.  It  is  a  function  of  engine  speed  and 
forward  velocity. 


D 


RAM 


DRAM(6N1 ’ 


M) 


(C-24) 


The  projection  of  this  force  on  body  axes  depends  on  the  aero¬ 
dynamic  angles 


-RAM  “  HW(° ) 


-D 


RAM 

0 

0 


(C-25 ) 
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The  total  thrust- induced  force  is  the  sum  of  the  reaction 
control  forces,  the  main  nozzle  forces  and  the  ram  drag. 


-B  — FP  +  -RP  +  -LR  +  -RR  +  -YAW  +  -F  +  -R  +  -RAM  (C'26) 

The  moment  vector  applied  to  the  vehicle  by  each  of  these 
forces  is  the  cross  product  of  the  location  of  the  point  of 
application  of  the  force  and  the  force  vector.  Using  Ax  to 
represent  the  location  of  the  point  of  application  relative  to 
the  aircraft  eg,  and  using  Ax  to  represent  the  cross-product 
equivalent  matrix,  the  thrust  moments  appear  as  follows: 


-B  '  AxFP  -FP  +  AxRP  ^RP  +  AXLR  -LR 


+  AxRR  -LR  +  AxYAW  ^YAW  +  AxF  -F 


+  AxR  -R  +  AxRAM  -RAM 


(C- 27) 


The  aerodynamic  forces  and  moments  depend  on  dynamic 
pressure  (this  model  is  not  designed  to  operate  at  hover),  but 
are  modified  by  the  engine- induced  flow  field.  This  flow 


field  is  defined  by  the  nozzle  angle,  6^^  ,  and  the  jet-to- 


forward  velocity  ratio,  V^/V. 


Vj/V  =  pv2  SN; 


(C-28 ) 


where  SXI  is  the  nozzle  reference  area  and  T  is  the  total  front 


The  aerodynamic  forces  in  body  axes  are  calculated 
from  wind-axis  force  coefficients  as  follows 


cos  a 

0 

-sin  a 

-% 

pV2 

s 

CD 

4  = 

0 

1 

0 

% 

pv2 

s 

CY 

sin  a 

0 

cos  a 

pV2 

? 

CL 

_ 

( 029 ) 


where  the  coefficients  are  calculated  f^om  tabular  values 
according  to  the  equations 


CD  =  cd  6s*  Vj//V>  60j^  +  CD  2V 
o  J  J  q 


(030 ) 


CY  =  CY  (p,  a,  V./V,  6  )  +  C  6  +  C  6 

Y  Y~  J  y6r  r  y-'a 


oa 


+  <C  Ps  +  C  rg)  jv 


(C-31) 


C,  =  C  (a,  6S,  V./V,  6  )  +  (CL.o  +  CL  q)  ^  (C-32) 

L  Lo  ^  ^  a  q 


The  stability  axis  roll  and  yaw  rates  are  found  as  follows 


cos  a  0  sin  a 

0  1  0 

■sin  a  0  cos  a 


P 

q 

r 


(033 ) 


The  aerodynamic  moments  are  calculated  in  a  similar  manner. 
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i 

1 


•im 


m 


cos  or  0 


-B  ' 


h  pv" 

h  pv: 


(C-3A ) 


sin  or  0  cos  or  Jg  pV* 


where 


O  =  C,  (p,  or,  V  A,  6,.)  +  C,  ( V  )6  +  C,  6a 

XS  Xo  J  6j  *6r  r  ‘6a  a 


+  (C^  ps  +  Cj^(a)rs) 


(C- 35) 


Cm  =  Cm  <«•  V  V*’  60j>  +  <Cm.“  +  Cm 

o  or  q 


m  q'  2V 

q 


-  ct-S*  -  -?2)  ♦  c  <-£8  -  -X£) 


where  x  =  0.09  c  ,  z  =0.0c 


(C-36) 


C  =  C  <p,  a,  V./V,  6fl.)  +  C  6  +  C  (o)6fl 

nS  no  J  0J  n6r  r  n6a  a 


+  ^Cnp(o,)Ps  +  CnrrS)  IV 


(037 ) 


For  use  in  the  study  reported  here,  the  nonlinear 
equations  described  above  are  used  to  calculate  vehicle  trim 
conditions.  Static  trim  is  straightforward,  since  all  of  the 
"inner  eight"  states  (6 ,u ,q ,w,v ,r ,p ,6 )  are  constant.  Quasi - 
static  trim  involves  a  force  and  moment  balance  which  produces 
local ly-constant  earth-relative  acceleration.  Specifically, 
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quasi-static  trim  involves  constant  V,  a,  £,  p,  q,  and  r  in 
this  application.  About  each  of  these  trim  conditions,  a 
linearized  aircraft  model  is  found  by  taking  into  account  all 
interactions  between  states  and  controls.  Analysis  of  this 
linear  system  provides  the  handling  quality  data  presented  in 
Chapter  2.  These  linear  systems  also  provide  the  basis  for 
flight  control  design  presented  in  Chapter  3. 


C . 2  GROUND  EFFECTS  MODEL 

The  ground  effects  included  in  this  study  are  due  to 
suck  down  (loss  of  lift  as  altitude  decreases)  and  fountain 
flow  effects  on  roll  and  pitch  moments. 

Reference  6  discusses  a  6  per  cent  lift  loss  during  a 
decent  from  20  ft  altitude  to  zero.  Assuming  the  thrust  at  20 
ft  is  about  equal  to  weight  and  the  suck  down  is  linear  be¬ 
tween  zero  and  20  ft  produces  the  following 


0.06  •  m* 
20  ft 


=  40  lbf/ft 


(C-38) 


=580  N/m 


Reference  17  includes  a  thorough  analysis  of  the 
effects  of  the  ground  on  pitch  and  roll  moments  for  the 
VAK-191B.  Since  the  configurations  of  the  VAK-191B  and  the 
AV-8A  differ  significantly,  the  following  adjustments  were 
made  to  Ref.  17  data  for  use  with  the  AV-8A.  The  AV-8A  has 
four  splayed  nozzles  on  its  lift-cruise  engine,  while  the 
VAK-191B  has  four  nozzles  on  its  lift-cruise  engine  and  two 
vertical  lift  engines.  The  lift-cruise  nozzles  of  the  VAK-191B 
do  not  appear  to  be  as  widely  splayed  as  those  of  the  Harrier. 
All  of  this  leads  to  the  subjective  judgement  that  roll 
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instability  will  not  be  as  severe  for  the  AV-8A  as  it  is  for 

the  VAK-191B.  However,  the  Harrier  roll  inertia  is  1.37  times 

the  VAK-191B  roll  inertia.  Hence,  the  same  value  of  L,  mea- 

$ 

sured  for  the  VAK-191B  (450  ft-lb^/deg)  is  used  for  the  Harrier, 
and  the  larger  roll  inertia  of  the  AV-8A  will  reduce  its  dynamic 
effect . 

The  stabilizing  pitch  stiffness  and  damping  terms 
observed  for  the  VAK-191B  are  probably  due  to  the  long  moment 
arm  of  the  lift  engines.  It  was  judged  advisable  to  reduce 
the  VAK-191B  values  by  50  per  cent  for  use  with  the  AV-8A. 

This  produces  a  pitch  stiffness  (Mg)  due  to  ground  effect  of 
-125  ft-lb^/deg  and  a  pitch  damping  (M^)  of  -250  ft-lb^/ 
(deg/sec).  Note  that  since  the  pitch  inertia  of  the  Harrier 
is  about  50%  smaller  than  that  of  the  VAK-191B,  the  actual 
dynamic  effects  on  the  two  aircraft  will  be  similar.  Table 
C-3  summarizes  the  aerodynamic  changes  due  to  ground  effect, 
as  modeled  in  this  study. 

TABLE  C-3 

GROUND  EFFECT  INCREMENTS 

T-3070 


zh 

40 

lbf/f t 

580 

N/m 

Me 

-125 

ft-lbf/deg 

-170 

N  • 

m/deg 

Me 

-250 

f t-lbf/(deg/sec ) 

-340 

N  • 

m/( deg/s) 

L* 

450 

ft-lbf/deg 

+610 

N  • 

m/deg 

C .  3 


AV-8A  CONTROL  SYSTEM  DESCRIPTION 


The  AV-8A  control  system  consists  of  the  pilot  con¬ 
trols,  the  control  linkage,  the  main  actuators,  the  control 
surfaces,  the  propulsion  system,  the  reaction  control  system 
(res),  and  the  stability  augmentation  system.  Except  for  the 
last,  these  constituents  of  the  control  system  are  discussed 
in  the  following  paragraphs. 

The  Pilot  Controls  -  The  pilot  controls  are  the  lat¬ 
eral  and  longitudinal  stick,  the  rudder  pedals,  the  throttle 
and  the  nozzle  angle  lever. 

The  lateral  stick  controls  the  ailerons  and  roll  res 
jets,  as  well  as  rudder  through  an  aileron-to-rudder  intercon¬ 
nect.  Lateral  stick  deflection  limits  are  -4.0  in  to  +  4.0  in 
with  positive  stick  giving  positive  aileron.  For  the  AV - 8 A , 
positive  aileron  produces  positive  roll  moments.  The  lateral 
stick- to-aileron  gearing  is  shown  in  Fig.  C-2. 


n-4?090 


Lateral  Stick-to-Ai leron  Gearing 
150 


Figure  C-2 


The  longitudinal  stick  controls  the  stabilator  and 
the  pitch  res  jets.  Longitudinal  stick  deflection  limits  are 
+3.7  in  to  -7.5  in,  with  positive  stick  giving  positive  sta¬ 
bilator.  For  the  AV - 8 A ,  positive  stabilator  produces  negative 
(nose  down)  pitch.  The  longitudinal  stick  to  stabilator  gear¬ 
ing  is  shown  in  Fig.  C-3. 


STABILATOR.  6S  Wegl 


1.54  deg /in 


Figure  C-3  Longitudinal  Stick-to-Stabilator  Gearing 


The  rudder  pedals  control  the  rudder  deflection  and 
the  yaw  res  jets,  with  pedal  limits  of  ±2.1  in.  Positive 
pedal  deflection  corresponds  to  positive  rudder  deflection 
(which  produces  nose  left  yaw  and  right  wing  down  roll. 
Rudder  pedal  to  rudder  gearing  is  shown  in  Fig.  C-4. 


R-4  2092 


Figure  C-4  Rudder  Pedal - to-Rudder  Gearing 

The  throttle  controls  engine  thrust  level,  and  is 
calibrated  directly  in  percent  of  fan  rpm  at  normal  rated 
thrust.  In  the  actual  aircraft  the  idle  limit  is  somewhat 
above  0%,  while  the  absolute  limit  is  105%.  There  are  further 
limits  (not  modeled  in  this  work)  of  102%  if  water  injection 
is  off  or  nozzle  angle  is  less  than  10  deg.  There  are  also 
temperature  limits  which  restrict  the  time  that  can  be  spent 
above  102%. 

The  nozzle  angle  control  (physically  located  next  to 
the  throttle)  is  calibrated  in  degrees  of  nozzle  angle.  The 
limits  are  zero  degrees  (nozzles  fully  rearward)  and  98.5  deg, 
(nozzles  fully  down  and  somewhat  forward).  These  angles  are 
measured  relative  to  a  line  which  is  1.5  deg  down  from  the 
body  x-axis.  Hence,  even  at  zero  nozzle  angle,  there  is  some 
upward  body  z-axis  force  due  to  thrust. 
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The  Control  System  Linkage  and  Main  Actuators  -  The 
aileron  main  actuator  is  driven  by  the  aileron  command  from 
lateral  stick  and  b}’  the  roll  sas  command.  The  aileron  com¬ 
mand  from  lateral  stick  is  limited  to  ±12  deg,  while  the  roll 
sas  limit  is  ±2  deg.  The  main  actuator  is  essentially  a  first 
order  lag  with  a  time  constant  of  0 . 04  sec  (1/25  sec).  No 
rate  limit  is  given. 

The  stabilator  main  actuator  is  driven  by  the  stabi- 
lator  command  from  longitudinal  stick  (limited  to  -10  to  +11.3 
deg)  and  by  the  pitch  sas  command  (±1.5  deg  limit).  The  main 
actuator  time  constant  is  0.0833  sec  (1/12  sec).  No  rate 
limit  is  given. 

The  rudder  has  no  actuator,  and  is  driven  by  rudder 
command  from  rudder  pedals  (±15  deg)  and  by  a  mechanical 
ai leron- rudder  interconnect  (0.58  deg  rudder/deg  aileron) 
driven  by  lateral  stick.  There  is  a  discrepancy  in  Ref.  6 
between  the  AV-8A  computer  simulation  and  the  rudder  control 
system  description;  the  above  description  corresponds  to  the 
AV-8A  computer  simulation. 

The  engine  is  modeled  as  a  first-order  system  with 
time  constant  and  rate  limits  which  depend  on  engine  rpm.  The 
time  constant  is  plotted  in  Fig.  C-5,  while  the  limit  on  posi¬ 
tive  rpm  rate  is  shown  in  Fig.  C-6.  There  is  no  limit  on 
negative  rpm  rate. 

The  nozzle  angle  drive  is  a  constant  rate  drive  that 
operates  at  150  deg/sec.  The  front  nozzle  splay  angle  is  con¬ 
stant  at  5.13  deg,  while  the  rear  nozzle  splay  angle  varies 

with  6.  as  shown  in  Fig.  C-7. 

J 
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ENGINE  TIME  CONSTANT  (sec) 


t 

I 


H-4209'j 
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0  20  40  60  80  100 

NOZZLE  ANGLE.  6i  (deg) 

Figure  C-7  Rear  Nozzle  Splay  Angle  as  a  Function 
of  Nozzle  Angle 

The  Reaction  Control  System  -  Jet  Effectiveness  -  The 
reaction  control  system  (res)  consists  of  front  and  rear  pitch 
jets,  left  and  right  roll  jets  and  a  rear  mounted  yaw  jet. 

The  locations  of  these  jets  are  shov/n  in  Fig.  C-l.  Note  that 
the  pitch  jets  only  thrust  upwards,  while  the  roll  jets  exhaust 
both  up  and  down.  The  yaw  jets  are  also  bi-directional.  The 
jets  are  supplied  with  air  bled  from  the  fan  plenum;  this  bleed 
decreases  the  available  thrust  from  the  engine  nozzles.  The 
bleed  air  connection  is  fully  closed  at  a  nozzle  angle  of  0 
deg,  and  gradually  opens  as  0^  increases  to  20  deg,  as  shown 
in  Fig.  C-8.  The  maximum  thrust  from  each  nozzle  at  100%  N^ 
and  no  duct  pressure  loss  is  given  in  Table  C-4 ,  along  with 
the  body  axis  moments  that  these  thrust  values  produce. 


RCS  THRUST  GAIN 


R-42096 


NOZZLE  ANGLE,  6i  (deg) 


gure  C-8  RCS  Thrust  Gain  at  Low  Nozzle  Angles 

TABLE  C-4 

RCS  THRUSTER  POWER 

T-3071 


THRUSTER  JET 

NOMINAL  THRUST 

MOMENT 

Front  Pitch 

1 

4230  N 

21,960  N  •  m 

Jet 

950  lb. 

16,198  lbf- ft  pitch 

Rear  Pitch 

4890  S 

-32,740  N  ■  m 

Jet 

1100  lbf 

-24.145  lbf-ft  pitch 

Y  av.  Jet 

2890  N 

-18.850  N  ■  m 

650  lbf 

-13.90u  lbf-ft  yav 

-1,730  N  •  m 
-1,310  lbf-ft  roll 

Roll  Jet 

2720  N 

9,190  N  •  m 

( Up  Thrust ) 

612  lbf 

6,781  lbf-ft  roll 

-5,320  N  ■  ra 
-3.926  lbf-ft  pitch 

Roll  Jet 

1060  S 

3.570  S  •  m 

'Down  Thrust) 

238  lbf 

2,632  lbc-ft  roll 

i. 

2.C70  N  •  m 

To  summarize,  the  nominal  thrust  is  changed  by  the 
plenum  pressure  ratio  (which  depends  on  N^),  by  main  bleed 
value  position  (which  depends  on  6^  and  is  fully  open  for  0^. 
greater  than  20  deg) ,  by  duct  pressure  losses  (which  depend 
on  the  open  exhaust  area  of  the  thrusters)  and  each  individual 
thruster  exhaust  area. 

These  last  thrust- to-area  relationships  are  detailed 
in  Figs.  C-9  and  C-10. 

R-42097 


Figure  C-9  Pitch  Jet  Thrust  Ratio  vs  Jet  Area 

The  Reaction  Control  System  -  Jet  Gearing  -  The  res 
jets  are  opened  and  closed  by  the  pilots  controls  and  the  sta¬ 
bility  augmentation  system.  The  front  pitch  jet  is  connected 


R-42099 


Figure  C-10  Roll  and  Yaw  Jet  Thrust  Ratio  vs  Jet  Area 

directly  to  the  longitudinal  stick,  as  described  by  Fig.  C-ll. 
The  rear  pitch  jet  is  driven  by  the  actual  stabilator  deflec¬ 
tion,  which  itself  is  a  function  both  of  pilot  stick  and  pitch 
sas.  The  rear  pitch  jet  to  stabilator  gearing  is  shown  in 
Fig.  C- 12. 


The  wing-mounted  roll  jets  are  connected  to  their 
respective  ailerons.  Figure  C-13  shows  the  form  of  this  gear¬ 
ing.  Note  that  the  ailerons  deflect  differentially,  hence 
the  net  effect  is  that  due  to  both  roll  jets.  The  yaw  jets 
are  driven  by  both  the  rudder  deflection  and  the  yaw  sas  out¬ 
put.  Interestingly,  the  yaw  sas  does  not  drive  rudder  at  all. 
The  yaw  sas  output,  expressed  in  equivalent  degrees  of  rudder, 
is  limited  to  ±5  deg.  Figure  C-14  gives  the  yaw  res  jet  gear¬ 
ing  as  a  function  of  the  rudder  and  yaw  sas  commands. 


Figure  C-13  Roll  Jet  vs  Aileron  Gearing 
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Figure  C- 14 


Yaw  Jet  vs  Equivalent  Rudder  Gearing 


REACTION  JET  INTERCONNECTIONS 


I 
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As  listed  in  Table  C-2,  the  AV-8A  model  analyzed  in 
this  report  has  ten  independent  controls.  For  the  purposes  of 
digital  flight  control  system  design,  however,  the  number  of 
controls  was  reduced  to  eight.  This  was  done  by  combining 
front  and  rear  pitch  jets  and  left  and  right  roll  jets  in  much 
the  same  way  as  they  are  combined  in  the  actual  aircraft. 


For  the  pitch  jets,  a  single  control  variable,  6 


PJ  ’ 


is  created  which  directly  commands  front  pitch  jet  with  unity 
gain  when  positive.  The  gearing  to  the  rear  pitch  jet  when 
6pj  is  negative  is  chosen  so  that  pitch  moment  due  to  6pj  is 
linear.  From  Table  C-4,  this  gain  is  as  follows: 


K 


-0.81 


(C-39) 


The  lower  limit  of  6 

PJ 

pitch  jet  as  follows: 


is  found  from  the  limit  of  the  rear 


6  .  .  =  ( 1/K )  6 

pj ,min  /  rp.max 

=  -10.78  in2  (C-40 ) 

The  roll  jet  gearing  is  chosen  to  be  linear  from  the 
minimum  to  the  maximum  deflection.  As  one  roll  jet  thrusts 
upwards,  the  other  thrusts  downwards.  Note  from  Fig.  C-13 
that  this  is  somewhat  different  from  the  AV-8A  roll  jet  gear¬ 
ing.  The  down- thrusting  roll  jet  is  much  less  effective  than 
the  up-thrusting  one,  so  the  roll  moment  excited  by  both  jets 
is  only  1.388  times  that  excited  by  the  up-thrusting  jet. 
Additionally,  the  vertical  forces  are  not  balanced  by  this 
gearing.  The  combined  roll  jet  control  is  referred  to  as  6rj 
in  Chapter  3. 


APPENDIX  D 
LIST  OF  SYMBOLS 


In  general,  matrices  are  represented  by  capital 
letters  and  vectors  are  underscored;  exceptions  to  these 
rules  are  only  made  when  they  are  contradicted  by  standard 
aerodynamic  notation. 


Variable 


A 

a 


a 


y 


B 

b 

C 


c 

D 

E 

F 


Description 

Discrete  time  feedforward  matrix 
Acceleration 

Component  of  the  earth-relative  accelera¬ 
tion  in  the  aircraft  x-z  plane  normal  to 
the  velocity  vector 

Component  of  the  earth-relative  accelera¬ 
tion  normal  to  the  aircraft  x-z  plane 

Disturbance  observation  matrix 

Wing  span 

State  observation  matrix 

Type  1  control  law  feedback  gains 

Partial  derivative  of  the  nondimensional 
coefficient  of  force  or  moment  1  with 
respect  to  the  nondimensional  variable  2 
( scalar ) 

Mean  aerodynamic  chord 

Control  observation  matrix 
Drag 


F 


Disturbance  input  allocation  matrix 
System  dynamics  matrix 
Aerodynamic  contact  force  vector 


LIST  OF  SYMBOLS  (Continued) 


Variable 


F 


CL 


f 


f 


G 


g 


H 

H 

h 

I 


i 

J 

j 

K 

k 

L 


In 

SL 

M 


m 


N 


Description 

Closed-loop  aircraft  system  matrix 

Vector-valued  nonlinear  function 

Scalar  element  in  system  dynamics  matrix 

Control  input  allocation  matrix 

Magnitude  of  gravitational  acceleration 
vector 

Command  variable  transformation  matrix 

Euler  angle  transformation  from 
Frame  1  axes  to  Frame  2  axes 

Altitude 

Identity  matrix 
Moment  of  inertia 

Index  integer 

Cost  functional  matrix 

v^T 

Type  0  gain  matrix 
Index  integer 

Aerodynamic  moment  about  the  x-axis 
( scalar) 

Natural  logarithm 
Number  of  commands 

Aerodynamic  moment  about  the  y-axis  (scalar) 
Cross  weighting  matrix  between  states  and 
controls 

Mass  of  the  vehicle 
Number  of  controls 
Meters 

Aerodynamic  moment  about  the  z-axis  (scalar) 

■  2 

Newtons  (kg  m  sec  ) 
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LIST  OF  SYMBOLS  (Continued) 


Variable 

N1 

n 

P 

P 

Py 

Q 

q 

^00 

R 

Ra 

r 

S 

T 

t 

u 

u 

V 

v 

w 

w 

X 


Description 

Engine  speed  expressed  as  a  percentage  of 
full  speed 

Number  of  states 

Riccati  matrix  in  the  optimal  regulator 
problem 

Rotational  rate  about  the  body  x-axis 

Stability-axis  roll  rate 

State  weighting  matrix 

Rotational  rate  about  the  body  y-axis 

2 

Free  stream  dynamic  pressure  (=^pVQ) 

Control  or  control-rate  weighting  matrix 

A  subspace  with  dimension  a 

Rotational  rate  about  the  bod}'  z-axis 

Reference  area  (usually  wing  area) 

State  command  matrix  used  in  eigenvalue/ 
eigenvector  placement 
(Diagonalizing  matrix) 

Thrust 

Time 

Body  x-axis  velocity  component 

Control  vector 

Inertial  velocity  magnitude 

Body  y-axis  velocity  component 
Control  difference 

Body  z-axis  velocity  component 

Aircraft  disturbance  vector 

Aerodynamic  force  along  the  x-axis  (scalar) 
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LIST  OF  SYMBOLS  (Continued) 


x 


Position  along  the  x-axis 


x 


State  vector 


x 


eg 


Y 


y 


z 


Normalized  longitudinal  distance  between 
actual  c.g.  location  and  point  used  for 
aerodynamic  moment  measurements  expressed 
in  body  axes 

Aerodynamic  force  along  the  y-axis  (scalar) 
Position  along  the  y-axis 

Aerodynamic  force  along  the  z-axis  (scalar) 


z 


Position  along  the  z-axis 
Laplace  operator  (discrete  time) 


Variable 

(Greek) 


a 


3 


r 

Y 

6 

6 

6 

6 

6 

6 


a 

PJ 

rj 

r 

s 

yj 


£ 


c 


Description 

Wind-body  pitch  Euler  angle 
(angle  of  attack) 

Negative  of  wind -body  yaw  Euler  angle 
(sideslip  angle) 

Discrete  time  control  effect  matrix 

Inertial-velocity  axis  pitch  Euler  angle 
(flight-path  angle) 

Aileron  deflection 

Pitch  jets  nozzle  area  opening 

Roll  jets  nozzle  area  opening 

Rudder  deflection 

Stabilator  deflection 

Yaw  jet  area  opening 

Reaction  jet  force 

Damping  ratio 
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laiail 


Inertial-body  pitch  Euler  angle 
Thrust  nozzle  angle 
Eigenvalue 

Inert ial -veloci ty  axis  yaw  Euler  angle 
(heading  angle) 

Eigenvector 
Integrator  state 

Matrix  in  quad  partition  matrix  inverse 

Air  density 

Summation 

Real  part  of  an  eigenvalue  in  radians/sec 
Time  constant 

Discrete- time  system  matrix 

Inertial -body  axis  roll  Euler  angle 

Matrix  in  adjusted  quad  partition 
matrix  inverse 

Inertial-body  axis  yaw  Euler  angle 

Frequency  in  radians/sec 
Imaginary  part  of  an  eigenvalue 


Variable 
(Subscript  or 
Superscript )  Description 


Eigenvalue/eigenvector  placement  models 
Body  axes 
Closed  loop 


Discrete  time 
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LIST  OF  SYMBOLS  (Continued) 


Variable 
(Subscript  or 
Superscript ) 

DR 

f 

I 

£ 

m 

N 

n 

V 

W 


w 

X 

y 


z 


X 


Y 


Z 


Description 

Downrange 

force 

Inertial  axes 

Aerodynamic  moment  about  the  x-axis 

Aerodynamic  moment  about  the  y-axis 
Model  to  be  followed 

Nozzle 

Aerodynamic  moment  about  the  z-axis 
Velocity  axes 

Wind  axes  (same  as  stability  axes  for 

^0  =  a0  = 

Disturbance 

Component  along  the  x-axis 
Component  along  the  y-axis 
Component  along  the  z-axis 
Aerodynamic  force  along  the  x-axis 
Aerodynamic  force  along  the  y-axis 
Aerodynamic  force  along  the  z-axis 


Operator  Definition 

(  )  Time  derivative 


(~)  Matrix  equivalent  to  vector  cross  product. 

Specifically,  if  x  is  the  three-dimensional 
vector 
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LIST  OF  SYMBOLS  (Continued) 


Operator  Definition 


X  = 

x 

y 

,  then  x  = 

'  0 

z 

-z 

0 

y“ 

-X 

z 

-y 

X 

0 

and  the  cross  product  of  x  and  f  is  equal 
to  the  product  of  the  matrix  x  and  the 
vector  f , 

xx  f  =  xf 


<:> 

Difference  between  a  vector  and  its 
expected  value 

(  )i 

The  i^  column  of  the  matrix 

(  )T 

Transpose  of  a  vector  or  matrix 

(  I'1 

Inverse  of  a  matrix 

<  >o 

Reference  or  nominal  value  of  a  variable 

A(  ) 

Perturbation  about  the  nominal  value  of 
a  variable 

^  ^max 

Maximum  value,  usually  due  to  displacement 
limit  of  an  actuator 

Ac  ronym 

Corresponding  Phrase 

AV-8A 

Harrier  jet  vertical  or  short  take-off 
and  landing  aircraft  in  operation 

AV-8B 
( YAV-8B ) 

Harrier  jet  vertical  or  short  take-off 
and  landing  aircraft  in  development 

CCV 

Control  configured  vehicle 

eg 

Center  of  gravity 

CGT-PI 

Command  generator  tracker-proportional 
integral  controller 

DFCS 

Digital  flight  control  system 

i 
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LIST  OF  SYMBOLS  (Continued) 


Acronym 

Corresponding  Phrase 

EEP 

Eigenvalue/eigenvector  placement 

IGE 

In  ground  effect 

kt 

Knots 

RPM 

Revolutions  per  minute 

res 

Reaction  control  system 

V/STOL 

Vertical  or  short  take-off  and  landing 
aircraft 

VAK-191B 

Jet  vertical  or  short  take-off  and 
landing  aircraft 

X-22A 

Fixed-wing  variable  fan  vertical  or  short 
take-off  and  landing  aircraft 
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